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Abstract 

The nearly Minimal Supersymmetric Standard Model (nMSSM) is one of the 
promising models of the new physics, since this model can avoid hierarchy problem, 
H problem, cosmological domain wall problem, and tadpole problem simultaneously. 
In this thesis, we consider the phenomenology of the nMSSM. Especially, we focus on 
the phenomenology of the dark matter and the baryon asymmetry in the universe 
generated by the electroweak baryogenesis mechanism. We find that with high- 
scale supersymmetry breaking the singlino can obtain a sizable radiative correction 
to the singlino mass, which opens a window for the singlet dark matter scenario 
with resonant annihilation via the exchange of the Higgs boson. We also propose 
a new electroweak baryogenesis scenario in the nMSSM with additional vector-like 
multiplets. If the soft supersymmetry breaking scale is 0(10) TeV, these scenarios 
are compatible with each other and an observed mass of the Higgs boson, constraints 
by the electric dipole moments measurements and the flavor experiments. As a result 
of these two studies, we conclude that the nMSSM with a high-scale supersymmetry 
breaking is valid and can be probed by the direct direction of the singlino dark 
matter. 
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Chapter 


1 

Introduction 


1.1 Overview 

The standard model (SM) of the particle physics has worked very well for a long time. Amazingly, 
ten observations and eight predictions or theories, which are related with the SM, had received 
the Novel Prize. In 2012 the Higgs boson, which is a last missing piece of the SM, was observed 
by the Large Hadron Collider (LHC) experiments at CERN [1,2], and F. Englert and P. Higgs 
have won the Nobel Prize for the discovery of the Higgs mechanism. This is a triumph of the SM, 
and the observation of the SM Higgs boson has given an important step towards understanding 
the electroweak symmetry breaking. 

However, there are many unsolved problems within the SM, for example, the observed dark 
matter particles and baryon asymmetry of the universe. From theoretical viewpoint, the gauge 
hierarchy problem is still in question. Hence, there have been many attempts to solve such 
problems in framework beyond the SM. 

The supersymmetric (SUSY) models are good candidates as the physics beyond the standard 
model [3-7]. It is because that they can solve the hierarchy problem naturally and ensure the 
unification of the gauge couplings. In addition, the lightest SUSY particle can be a natural 
candidate of the WIMP dark matter if the R parity is conserved. 

The minimal SUSY extension of the SM (MSSM) contains a supersymmetric dimensional 
parameter fi, which is the mass term of the superpartner of the Higgs boson. However, this 
parameter causes problem”, which is also one of the hierarchy problem [8]. Although /r has 
to be a size of the SUSY breaking scale to realize the electroweak symmetry breaking properly, 
there is no reason for /r to be small compared to the Planck scale. 

One of the simplest ways to solve the problem is introducing a gauge-singlet superfield [9]. 
There are several models of singlet extension of the MSSM depending on the imposed addi¬ 
tional symmetry. However, the additional singlet superfield causes a cosmological domain wall 
problem [10,11] and tadpole problem [12]. The nearly Minimal (or new Minimal) Supersym- 
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metric Standard Model (nMSSM) [13-15] is based on a discrete R-symmetry. Actually this 
model can avoid naturally the cosmological domain wall problem and tadpole problem, unlike 
Z 3 symmetric models [11,16]. Therefore, the nMSSM is one of the promising models of the 
new physics: this model can avoid the /r problem, the domain wall problem, and the tadpole 
problem simultaneously. In addition this model has natural candidate of the dark matter and 
can generate the baryon asymmetry of the universe. 

In this thesis, we consider the phenomenology of the nMSSM, which are the based on the 
works by the author [17,18]. Especially, we focus on a phenomenology of the dark matter [17] in 
Chapter 4 and the baryon asymmetry in the universe generated by the electroweak baryogenesis 
mechanism [18] in Chapter 5. 

In Chapter 4, we consider a singlino dark matter scenario in the nMSSM. The singlino is 
a fermion component of the additional singlet superfield. We find that with high-scale SUSY 
breaking the singlino can obtain a sizable radiative correction to the mass, which opens a window 
for the dark matter scenario with resonant annihilation via the exchange of the Higgs boson. We 
show that the current dark matter relic abundance and the Higgs boson mass can be explained 
simultaneously. This scenario can be probed by the search of the Higgs invisible decay and the 
direct direction of the dark matter. 

In Chapter 5, we propose a new electroweak baryogenesis scenario in high-scale SUSY mod¬ 
els, and consider the nMSSM introducing additional vector-like multiplets. We show that the 
strongly first-order phase transition can occur at a high temperature comparable to the soft 
SUSY breaking scale. In addition, the proper amount of the baryon asymmetry of the universe 
can be generated via the lepton number violating process in the vector-like multiple! sector. 
The typical scale of our scenario, the soft SUSY breaking scale, can be any value. Thus our 
new electroweak baryogenesis scenario can be realized at arbitrary scales and we call this sce¬ 
nario as a scale free electroweak baryogenesis. This soft SUSY breaking scale is determined by 
other requirements. If the soft SUSY breaking scale is 0(10) TeV, our scenario is compatible 
with the observed mass of the Higgs boson and the constraints by the electric dipole moments 
measurements and the flavor experiments. Furthermore, the singlino can be a good candidate 
of the dark matter. 

As a result of these two studies, we will conclude that the nMSSM with a high-scale SUSY 
breaking is valid and can be probed by the direct direction of the singlino dark matter. 

1.2 Organization of this thesis 

This thesis is organized as following. 

In Chapter 2, we review the SM, the supersymmetry and current status of the supersym¬ 
metric minimal model. In Section 2.1 and 2.2, in order to solve the hierarchy problem of the 
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standard model, we first introduce the supersymmetry and the MSSM. In Section 2.3, we review 
the current situation of the MSSM. In fact, an observed mass of the SM Higgs boson is 125 GeV, 
and it gives a meaningful constraint on the parameter space of supersymmetric models. Thus, we 
focus on the one-loop, two-loop and higher-loop radiative corrections to the Higgs boson mass. 
These calculations for the Higgs boson mass are reused in the study of the nMSSM (Chapter 4). 
Furthermore, we also discuss the constraints from the flavor violation and CP violation process 
in supersymmetric model. In Section 2.4, we summarize the current status of the MSSM. 

In Chapter 3, we review singlet extension models of the MSSM and the nMSSM. In Sec¬ 
tion 3.1, we first explain the g. problem in the MSSM, and in order to solve the g. problem we 
introduce the additional gauge singlet superfield. Next, we show that when one imposes extra 
symmetries to forbid unwanted terms of singlet superfield, these symmetries lead to the domain 
wall problem and the tadpole problem in Section 3.2. The nMSSM is the one of the models 
which can solve the /r problem, the domain wall problem and the tadpole problem, and so we 
review the nMSSM in Section 3.3 and 3.4. Anyhow, we need an extra symmetry to solve these 
problems. 

In Chapter 4, we study the phenomenology of the singlino dark matter in the nMSSM. In 
Section 4.1, we briefly review a situation of the singlino dark matter in the nMSSM, and we also 
explain why we have considered it. In Section 4.2, using the low energy effective Lagrangian 
we calculate thermal relic abundance of the singlino dark matter which annihilate via the SM 
Higgs boson. We point out that one-loop corrections to the singlino mass can raise its mass with 
relatively high-scale SUSY breaking, in Section 4.5. In Section 4.4, we numerically investigate the 
singlino resonant dark matter scenario with high-scale SUSY breaking, and show this scenario 
is compatible with the observed SM Higgs boson mass. Section 4.5 and 4.6 are devoted to the 
conclusion and discussions in this chapter. 

In Chapter 5, we study the baryon asymmetry in the universe generated by the electroweak 
baryogenesis mechanism in the nMSSM, and we propose a new electroweak baryogenesis scenario 
in the nMSSM with high-scale SUSY breaking. In Section 5.1, 5.2 and 5.3, we introduce the 
model: the nMSSM with vector-like multiplets, and also present the overview of our scenario. In 
Section 5.4, we discuss about the strongly first-order phase transition and this section is divided 
into three parts. In subsection 5.4.1, we show the full thermal potential at high temperatures. 
In subsection 5.4.2, we provide an intuitive understanding for the behavior of the potential at 
high temperatures. In subsection 5.4.3, we analyze the full potential and show that the strongly 
first-order phase transition actually occurs at a temperature comparable to Msusy- We also 
show that the region with low tan /3 and a light charged Higgs boson is favored in our scenario. 
In Section 5.5, we demonstrate the generation of the BAU with the lepton number violating 
process. In Section 5.6, we discuss the singlino dark matter scenario paying particular attention 
to the lifetime. Section 5.7 and 5.8 are devoted to the conclusion and discussions in this chapter. 

Chapter 6 is devoted to the conclusion of this thesis. 
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Chapter 



An Introduction to Supersymmetry 

In this chapter, in order to solve the hierarchy problem of the standard 
model, we first introduce the supersymmetry and the minimal supersymmetric 
standard model. In fact, an observed mass of the SM Higgs boson is 125 GeV, 
and it gives a meaningful constraint on the parameter space of supersymmetric 
models. Therefore, we first review the one-loop, two-loop and higher-loop ra¬ 
diative corrections to the Higgs boson mass. Furthermore, we also discuss the 
constraints from the flavor violation and CP violation process in supersymmet¬ 
ric model. These facts imply that the naive low scale minimal supersymmetric 
standard model, which is in spite of being favored in terms of the natural¬ 
ness, is disfavored. One of the solutions of these problems is the high-scale 
supersymmetry. 


2.1 The Standard Model 

In nature, there are four fundamental forces: the electromagnetic, weak, and strong nuclear and 
gravitational interactions. These interactions can be understood by interactions of elementary 
particles. The standard model (SM) can describe the electromagnetic, weak, and strong nuclear 
interactions as a quantum field dynamics. 

The standard model is one of the gauge theories. In the standard model, the imposed 
gauge symmetry is Gsm =U(l)yX S\J{2)lX SU(3)c [19-22]. Gauge bosons are introduced at 
every gauge symmetry: a B boson for U(l)y gauge, isospin triplet W bosons for SU(2)x, gauge 
and color gluons octet for SU(3)c gauge. However, in order to describe the electromagnetic 
and weak nuclear interactions, U(l)y x SU(2)i gauge symmetry (electroweak symmetry) should 
be spontaneous broken to gauge symmetry (electromagnetic symmetry) by the Higgs 

mechanism [23-26]. This mechanism predicts an existence of the Higgs boson h (SM Higgs 
boson). The unitarity requirement for the high-energy scattering of the longitudinal W boson 
leads to the upper bound on the mass of the SM Higgs boson nih < 700 GeV [27]. 


The standard model has been worked very well for a long time, and its last missing piece, the 
Higgs boson, was finally discovered by the LHC experiment at CERN [1,2]. This is a triumph 
of the SM and a great step to understand physics at the electroweak scale. 
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The Lagrangian of the standard model is given as follows, 


C = C 


gauge 


+ 


fermion T ■^Yukawa “ 1 “ ■^scalar; 


+ 


with 


C 


gauge 


3 8 

-1 E - 7 E 


( 2 . 1 ) 


( 2 . 2 ) 


a=l 


a=l 


'^fermion — 


iQ^r (^5^ - j Q, 

+iUiY (^9^ - U, + Wa^ + i^-B^ - i|A“G“^ A 

+iLir (^9^ + i^B^ - U + iEa^^ (9^ + ig’B^) A, (2.3) 


-^Yukawa 

-^scalar 


— Ui{yu)ijHQj Qj Lj H- H.c., 

J 




On - i^Bmu - ) H 


- V{H), 


(2.4) 

(2.5) 

( 2 . 6 ) 


where the field strength is defined as = 9^^“ — 9 j, 74 “ + gAf°'^’^A^^A'^, cj“ (a = 1, 2, 3) is 
the Pauli matrix, A“ (a = 1, 2,... , 8 ) is the Gell-Mann matrix, index i represents the generation 
(i = 1,2,3), y is the Yukawa couplings, the gauge couplings for U(l)y, SU(2)i and SU(3)c are 
denoted as 5 ', g and gs respectively. The contraction of the two SU(2) doublet is AB = A'^ia'^B. 
A more detail definition is written in Appendix A. H is the Higgs doublet, 


H = 


H+ 


The Higgs potential can be written as follows, 

V = -fi^\H^H\ + 


Aquartic 




Then, the Higgs doublet obtains the vacuum expectation value vew (VEV), 

2 _ 

''Jew — 


(2.7) 


( 2 . 8 ) 


(2.9) 


Aquartic 

It breaks the electroweak symmetry to the electromagnetic symmetry, SU( 2 ) 2 ,x U(l)y —>■ 
U(1 )em- Using the freedom of SU(2) rotations, one can always align the VEV with the neutral 
Higgs direction, 

™ ) + (;),('.+ iG») )• 

where G is the Numbu-Goldstone bosons and h is the SM Higgs boson. Then, the mass of the 
SM Higgs boson is 

— 2/i = 2 Aquartic(2-11) 


H = 
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The Higgs VEV vew generates the Dirac mass to the all fermion via the Yukawa interaction, 
nif = VfVEW- The gauge bosons also obtain the mass, 


f7'2 -I- 


'^EW^ 


71 ^2 9^ 2 


Tfl^ — ^gluon — 0, 


( 2 . 12 ) 


where ^( 7 ) = cos + sin Z = —smOwB + cos 6wW^, = (VT^ =F ilY^)/\/ 2 , and 

9w is the Weinberg angle. 

Therefore, the Higgs vev can determine the electroweak scale. In fact, vew = 174.10363 ± 
0.00004 GeV, which has been given by the measurement of the Fermi constant from the muon 
decay. In the other words, the parameter in the Higgs potential determines the electroweak 
scale and the SM Higgs boson mass. 


Hierarchy problem 

However, the radiative corrections to the parameter p? contain the quadratic divergence. 
For example, the top quark loop gives the following quadratic divergence, 

= + ■■■). ( 2 . 13 ) 

where A is the ultraviolet momentum cutoff. All diagrams which contain the quadratic diver¬ 
gence are shown in Figure 2.1. If one assume that the standard model is valid up to the Planck 
scale, the ultraviolet cutoff is naively the Planck scale, A ~ 10^® GeV. Although these radia¬ 
tive corrections can be renormalized by the bare parameter /Tq, it requires incredible fine-tuning 
cancellation between the bare mass Pq and the quadratic radiative corrections Ap'^. In fact, 
in order for the 0(100) GeV electroweak scale to be realized, the 10“^^ fine-tuning is required. 
This difficulty in the standard model is called the hierarchy problem. 

2.2 Supersymmetry and Minimal Model 

The supersymmetry (SUSY) is the one of the symmetry which can solve the hierarchy problem. 
In this section, we first introduce the supersymmetry and its minimal model, and show that the 
hierarchy problem is solved actually. Next, we briefly review the Lagrangian and the mass of 
the Higgs boson in the supersymmetric minimal model. 

2.2.1 Motivations of Supersymmetry 

First, let us briefly introduce the supersymmetry, which can solve the hierarchy problem [3-7]. 
The supersymmetry is the extension symmetry of the Poincare group through the introduction 
of anticommuting spinor generators Qa, Qa, where a, d are spinor index [28]. It is equal to the 
symmetry between the boson and fermion. 
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Standard model 


Supersymmetric minimal model 



Figure 2.1: The quadratic divergence to the mass of the Higgs boson in the standard model and 
the supersymmetric model. Since the supersymmetry assures that the Higgs quartic coupling is 
related to the gauge boson coupling like Eq. (2.20), all quadratic divergence are always canceled 
out in the supersymmetric model. 


Q|boson) = |fermion) (2-14) 

Qlfermion) = |boson). (2-15) 

Then, ordinary space-time is extended to the superspace {x^,0a,0a), where 0a,0a are the 
anticommuting Grassman coordinates. The supersymmetric theories are described by chiral 
superfield, vector superfield and spinor chiral superfield. The chiral superheld contains the 
scalar and fermion field on the superspace, 

l>(x, 9, 9) = cj>{y) + V29^l;{y) + 99F{y), (2.16) 

with y^ = x^ + i9a^9, and the hat (') represents the superheld. 0 is a complex scalar held 
(sfermion), ijj is two-components Weyl spinor held and F is an auxiliary held. The vector 
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Table 2.1: All chiral and vector superfields with their components for spin 0, 1/2 and 1, and 
their representations for SU(3)cX SU(2)£,x U(l)y gauge group in the MSSM. 

Chiral Supermultiplet 

SpinO Spin! SU{‘i)c SU{2 )l U{1)y 

Quark-Squark Q 

U 

D 

Lepton-Slepton L 

E 

Higgs-Higgsino Hi 

H2 

Q = {UL, dhf Q = [uL^dif 

Ur 

dn 

L={Va,eLf L={ve,eLY 

V e^R 

Hi = {Hi Hif Hi = {Hi Hif 
H2 = {Hl,Hlf H2 = {Hl,Hlf 

3 2 1 

3 1 -1 

3 1 i 

1 2 -i 

111 

1 2 -i 

1 2 i 

Vector Supermultiplet 

Spin i Spinl SU{3)c SU{2 )l U{1)y 

Gluon-Gluimo Vq 

W boson-Wino IV 

B boson-Bino Vb 

~g 

1 U±, lUO W±, 

8 10 

13 0 

110 


superfield V°’ contains the vector boson and fermion field on the superspace, 

e, 9) = + m6X‘^{y) - me\%y) + {D\y) - id^V^{y)) , (2.17) 

where a is the index of the generator of the gauge group, is the gauge boson, Xa is two- 
components Weyl spinor field (gaugino) and Da is an auxiliary field. The spinor chiral superfield 
1 T“ also contains the vector boson and fermion field on the superspace. 


w^{x,e,e) 


-iKiv) + 



ey + eeard^M\y), 


(2.18) 


where - d^V^. 

The supersymmetric minimal model is called Minimal Supersymmetric Standard Model 
(MSSM). Table 2.1 shows all chiral and vector superfields with their components for spin 0, 
1/2 and 1, and their representations for SU(3)cX SU(2)iX U(l)y gauge group in the MSSM. 
Here the tilde Q represents the SUSY partner of the SM particle. Note that, although the 
number of Higgs doublet is one in the SM, we must introduce two Higgs doublets in the SUSY 
model. 


Hi = 


Hf 

Hi 




(2.19) 


where H is called Higgsino. It is because that an existence of the Yukawa interaction with 
quark/lepton and the gauge anomaly cancelation require two kinds of the Higgs doublets that 
the hypercharge is opposite. 
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As we will discuss in detail later, the Higgs quartic coupling Eq. (2.8) and the gauge coupling 
are related. In the MSSM, this relation is given as 

Aquartic = 5"(i“)". (2.20) 

Figure 2.1 shows the quadratic divergence to the mass of the Higgs boson in the MSSM. Obvi¬ 
ously, the quadratic divergence by the top quark loop Eq. (2.13) is canceled out by the stop loop 
which are SUSY partners of the top quark. Furthermore, the relationship Eq. (2.20) can cancel 
out the quadratic divergence by the gauge boson, Higgs boson, gaugino and Higgsino. Actually, 
all quadratic divergence are always canceled out in the supersymmetric model [6,7]. Therefore, 
the supersymmetry can solve the hierarchy problem. 


Other motivation of the supersymmetry is dark matter. The SM does not include dark matter, 
which is stable and does not interact with the electromagnetic force. In the supersymmetric 
model, in order to forbid all harmful terms, which break baryon number (B) and lepton number 
(L) and thus cause the proton decay, one should add a new symmetry. This symmetry is called 
the R parity [29]. It is defined as 

Pj^ = (_i)3(B-L)+2s^ ^2.21) 


where s is the spin and B = 1/3 for Q, B = —1/3 for U, D, B = 0 for all others, L = 1 for 
L, L = — 1 for E and L = 0 for all others. Thus, the SM fermions, Higgs bosons and gauge 
bosons have even R parity (Pr = -|-1), while the squarks, sleptons, Higgsinos and gauginos 
have odd R parity (Pr = — 1). This symmetry is equal to discrete Z 2 R-symmetry. If the R 
parity is exact symmetry, the lightest supersymmetric particle (ESP) becomes stable. Here the 
supersymmetric particles are defined as the R parity odd ones. Therefore, when the LSP does 
not have the electromagnetic charge, it can be a natural candidate of the dark matter. 

Another motivation of the supersymmetry is the gauge unification. In fact, the electromag¬ 
netic and weak nuclear interactions are unified by the electroweak theory above the unification 
energy 0(100) GeV, it is so-called Glashow-Weinberg-Salam theory [19,20]. In this sense, the 
electroweak interaction and strong interaction may be unified by grand unification theory (GUT). 
In fact, although the gauge couplings can not unify in the SM, in the SUSY model the unifica¬ 
tion of the gauge couplings can occur [30]. Let us briefly observe this fact. The one-loop level 
renormalization group equations (RGEs) for the gauge couplings in the SM are 


dg' ^ 1 41 ,3 dg ^ 1 19 3 dgs ^ _^_7 3 

dlnQ ( 47 r )2 6 ^ ’ dlnQ 6 ^ ’ dXnQ (dvr)^ ’ 


( 2 . 22 ) 


where Q is the renormalization scale. While, the corresponding equations in the MSSM are 

^ = ^11P3 = _^33 ( 2 . 

dlnQ (AttY ’ dlnQ (dTr)^ ’ dlnQ (dTr)^ ^ 
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Logio(Q/GeV) 


Logio(Q/GeV) 


Figure 2.2: Two-loop level running gauge couplings as a function of the renormalization scale 
Q. The black lines represent y^5/3g', the blue lines represent g and the red lines represent gs- 
Left : The dashed lines correspond the running gauge couplings in the SM, while the solid lines 
correspond the running gauge couplings in the MSSM with Msusy = lOTeV. Right : The 
running gauge couplings in the MSSM with Msusy = 1 TeV (upper) and Msusy = 100 TeV 
(lower). 




Logio(Q/GeV) 


Logio(Q/GeV) 


Figure 2.3: Two-loop level running Yukawa couplings {yt- red, yf,- blue, yr- green) and the 
Higgs quartic coupling (Aquartic^ black) as a function of the renormalization scale Q. The dashed 
lines correspond the running couplings in the SM, while the solid lines correspond the running 
couplings in the MSSM with Msusy = 10TeV. We take tan/3 = 2 in the left panel, and 
tan/3 = 50 in the right panel. 
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Two-loop level RGEs for all couplings are summarized in Appendix B.l. 

Figure 2.2 and 2.3 show the two-loop level running couplings as a function of the renormal¬ 
ization scale Q. In Figure 2.2, the black lines represent ^/5/3g', the blue lines represent 
g{Q) and the red lines represent gsiQ)- In the left panel, the dashed lines correspond the run¬ 
ning gauge couplings in the SM, while the solid lines correspond the running gauge couplings 
in the MSSM with Msusy = lOTeV, here Msusy represents typical SUSY particles mass 
scale. As one can see, in the MSSM three gauge couplings are actually unified at GUT sale: 
A/gut — 2.0 x 10^® GeV, while in the SM it can not achieve. This fact gives one of the main 
motivations for the SUSY GUT models [31,32]. In the right panel, we show the running gauge 
couplings in the MSSM with Msusy = ITeV (upper) and Msusy = 100TeV (lower). We find 
that when A/susy ^ 0(100) TeV, the gauge coupling unification is no longer achieved if one 
takes universal SUSY particle mass in the MSSM. Thus, this figure implies that when mass 
scales of all SUSY particles are the same order, the GUT suggests that typical SUSY particles 
mass scale is A/susy ^ 0(100) TeV. 


In Figure 2.3, the red lines represent the running top Yukawa coupling yt{Q), the blue lines 
represent the running bottom Yukawa coupling yh{Q), the green lines represent the running 
tau Yukawa coupling yriQ) and the black lines represent the running Higgs quartic coupling 
Aquartic(f5)- The dashed lines correspond the running couplings in the SM, while the solid lines 
correspond the running couplings in the MSSM with Msusy = 10 TeV. We take tan/I = 2 in 
the left panel, and tan /3 = 50 in the right panel. The definition of tan [3 is given in next section. 
The important point in these figures is that the Higgs quartic coupling becomes negative at the 
high energy scale in the SM. If one take my ~ 125 GeV, then the Higgs quartic coupling becomes 
negative at Q A 10^^“^^ GeV [33]. This situation is avoided in the supersymmetric theory. It 
is because that the relation between the quartic coupling and the gauge couplings Fq. (2.20) 
assures the stability of the Aquartic- 


2.2.2 Lagrangian and the Higgs Boson Mass in the MSSM 


Next, we brief review the Lagrangian and the Higgs boson mass of the MSSM [34]. 

Using the chiral superfield Fq. (2.16), the vector superheld Fq. (2.17) and the spinor chiral 
superfield Fq. (2.18), the SUSY invariant Lagrangian of the MSSM is given as follows. 


/' 


£ = / de^S^K 




+ 


dO'^W 


$ 


+ H.C. + 


i/ 


dO'^W'^W^ + H.c. 


(2.24) 


where K is the following Kiihler potential 


K 


^,¥,V 


^$^Fxp(5t“U“)4>, 

!> 


factor ^5/3 is a normalization factor for the GUT. 


(2.25) 
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and W is the following superpotential 

W = PH 2 H 1 + V^Yukawa, (2.26) 

WVukawa = UyuQH 2 - DydQHi - EyeLHi, (2.27) 

where p is the supersymmetric mass of the Higgs multiplets. When the supersymmetry is 
broken, p becomes the mass of the Higgsinos. While the superpotential WVukawa gives the 
supersymmetric Yukawa couplings of the standard model. In fact, the Yukawa matrices yu, yu 
and ye are given approximately as follows 

/0 00 \ /0 00 \ /0 00 \ 

yu~ 0 0 0 , yd~ 0 0 0 , ye~ 0 0 0 . (2.28) 

VO 0 yt) Vo 0 Vb) Vo 0 Vr) 

Here and the following, we use these approximations of Yukawa matrices. Then, the superpo¬ 
tential becomes 


WYukawa = ytiRQH2 “ UbbRQHl - yrTRLHi 

= ytitRiiH^ - IrIlH^ ) + yb^^RbiHi - hRiiH^) 

+ yT{TRTLHl - TRVrHi). (2.29) 

Finally, as the spinor chiral superfields IF“ we take the U(l)y gauge multiplets (the B boson 
and the bino), the SU(2)i gauge multiplets (the W bosons and the winos) and the SU(3)c gauge 
multiplets (the gluons and the gluinos). 

Because the supersymmetry is broken in nature, one has to also introduce SUSY breaking 
terms in Lagrangian. The supersymmetry must be broken softly in order to ensure the cancela¬ 
tion of the quadratic divergences. Generally, the soft SUSY breaking terms are gaugino masses, 
scalar masse, and trilinear coupling terms for scalars. In the MSSM, the soft SUSY breaking 
Lagrangian is given as 


•^soft - I^oft gaugino - I4oft Yukawa I^oft Higgs j 


(2.30) 


with 


^oft gaugino 

= ^ (^MiBB + M 2 WW + Msgg^ + H.c., 

(2.31) 

^oft Yukawa 

= + + m])b*b + mlVL + m\E*E 



+ {by^A^Q ■ H 2 - by^A^Q ■ Hi - by^AeL ■ Hi + H.c.^ , 

(2.32) 

^oft Higgs 

= ml\Hi\^ + ml\H 2 \^ + ^{BpH 2 Hi + H.c.), 

(2.33) 


where Mj are the gaugino masses, jDiLjE squark/slepton masses, m ^/2 

SUSY breaking Higgs mass, Ai are called A terms and B is called B terms. Note that m? 
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are flavor generation mixing 3x3 matrix generally. Without additional assumptions, these off- 
diagonal squark/slepton masses are not suppressed. Adding the soft SUSY breaking Lagrangian 
Eq. (2.30) into the SUSY invariant Lagrangian Eq. (2.24), one can obtain the full Lagrangian 
of the MSSM. 


Now, let us see the Higgs sector and the mass of the neutral Higgs boron of the MSSM. The 
tree-level Higgs scalar potential is given directly from Eqs. (2.30, 2.24) as 


Uni, 


2 


HIH2 


+ 


g" + g' 


- \H2?f + (l/rp + ml) + (j^j^ + mi) 




(2.34) 


Using the freedom of SU(2) rotations, one can always choose {H 2 ) = 0. Then, a potential 
minimization condition dV/dH^ = 0 leads to {H^) = 0. Thus, scalar potential of the neutral 
Higgs boson is given as 


Uh 


iggs 


g^ + 


1 


(l^iT - + {\g\^ + ml) + (|/r|2 + mi) \hI 


0|2 


The potential minimization conditions dV/dHf = 0 and dV/dH^ = 0 leads to 

(mi + ^ -\v2\'^)v*i-^BfiV2 = 0, 


1 


(mi + |/ip)u2 - ^ -\ v 2 \'^) v 2 -= 0 , 


(2.35) 

(2.36) 

(2.37) 


where Ui /2 is the dev of If th® determinant of d'^V/{dH^dH^) is negative, one linear 

combination of and Lfg a negative square mass, then and H 2 obtain nonzero vevs. 
This condition at the vicinity of = H 2 = 0 is 


> 4(mi \g\%m^2 + \g'\^)- 


(2.38) 


In this case, using a redefinition of the phase of and H 2 , one can always choose that Bfj,, vi 
and V 2 are real and positive. Let us defined the ratio of the Higgs vev as 


tan P = 


ui’ 


vew = — 174.1 GeV. 


Then, the potential minimization conditions Eqs. (2.36, 2.37) become 

ml = — I/ip-|-— - (M^-|-M|) cos 2/3, 


(2.39) 

(2.40) 

(2.41) 
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^2 ~ d- -M\ + - (^M\ + M^') cos2f3, (2.42) 

where 


M\ = 2|//|^ + mfy ml- 


(2.43) 


The vacuum fluctuations of the Higgs field are dehned as 


= 


vi + + iHii), 

^2 = V2 ~^^{H2R + iH2l)- 

Now, one can obtain the following mass matrix for the CP even Higgs bosons 
1 


Ml sin^ ,5 + M| cos^ ,5 -i(M| + M|) sin 2/3 

— ^(M^ + M|) sin 2/3 cos^/3 + M| sin^/3 


- ^CP -even — 2 {HirH2r) 

Diagonalizing this mass matrix, the following mass eigenvalues are given^^ 
mjj = 

ml = 


(Hir 

\H2R 


2 + \/ i^A ~ + 4M^M| sin^ 2/3 

i (mI + M| - ^(Ml - M|)2 + 4MlM| sin^ 2/3 


(2.44) 

(2.45) 

(2.48) 

(2.49) 

here we call a eigenstate, which has the heaviest/lightest mass eigenvalue, H/h. In the decou¬ 
pling limit Ml ;§> Ml, the eigenstate h becomes the SM Higgs boson perfectly. It is because 
that the couplings of h with SM particles are equal to the one of the SM Higgs boson with SM 
particles in this limit. Here and the following, we call this eigenstate h SM-like Higgs boson. 
Thus, one can recognize Eq. (2.49) as the tree-level mass of the SM-like Higgs boson in the 
MSSM. 

In fact, the tree-level mass of the SM-like Higgs boson Eq. (2.49) has an upper bound, 

m\ < MI cos^ 2/3. (2.50) 


In the SUSY model, in order to cancel out the quadratic divergence in the radiative correc¬ 
tions, the Higgs quartic couplings must be related with the gauge couplings (see Figure 2.1 and 
Eq. (2.20)). In addition, the Higgs quartic couplings decide the mass of the SM-like Higgs boson 
(see Eq. (2.11)). These are the reason why there are the upper bound on the Higgs boson mass 
in the SUSY model. Thus, the SM-like Higgs boson mass must be lighter than the Z boson mass 
at tree-level in the MSSM. However, this mass bound is too severe to explain an observed SM 
Higgs boson mass, m/j ~ 125 GeV [1,2]. 

’^^Similarly to the CP even Higgs boson mass, the tree-level mass of the CP-odd Higgs boson A and charged 
Higgs boson are given as 

Ml = + = (2.46) 

(2.47) 
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Figure 2.4: The Feynman diagrams for the dominant one-loop radiative corrections to the Higgs 

boson mass. 

Here, the bullet • represents At term interaction. 





Radiative corrections to the Higgs boson mass 

On the contrary, in 1990, several groups found that radiative corrections can contribute to 
the mass of the Higgs boson significantly [35-40] The dominant radiative correction comes 
from the top/stop loop diagrams, which are depicted in Figure 2.4. 

The top/stop one-loop effective potential (the Coleman-Weinberg potential [42]) is given 
as. 



(2.51) 


where the top quark mass is = y||R 2 P ^.nd Q is the renormalization scale. The stop masses 
mi^ (lighter) and (heavier) are given as the eigenvalues of the following stop mass matrix. 


with 


Mf tr' 

r — /j.* j.* \ I til ti2 n 

-C-stop mass — I ]^2 

\ i21 l22. 



(2.52) 


Mf 

til 


Mi = 


t22 


M/ 

t\2 


- J#) (|Hf- IffSP) 

_L ffuiz , I 2 2 M% 


3 fEVI/ 3 v\^r 


l^iT - \Hl?) , 


f‘2 \ * 


yt {AKHlY - ^H«). 


(2.53) 


This effective potential gives the following one-loop corrections to the CP even Higgs mass 
matrix Eq. (2.45), 

that time, the Large Electron-Positron Collider (LEP) had not excluded the Higgs boson mass up to the 
Z boson mass region [41]. 
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AMji 


4 2 2 d2/ 2 2 \2 / 2 2 \ 


+ 


AM^ 


22 


+ 


IOtt^ 

3 „, 4^,2 
Stt^ 

3 


VtAtp tan, m 




(2.54) 


2 2 

yt '^2 ( In —+ ‘^AtRt In 


m~ 


IGtt^ 


mj 

^ 12 
y'^Atp cot j3f{Q^,ml, 


+ A Rt {mp - ) 5(m^- , ) 


t2' 


(2.55) 


AM^ 


12 


2 

mf 


= --^ytV29Rt(ln-A- + AtRtim^^ - 


2 
m~ 

C 2 


^2 ^ 


t2 ■ 


lOvr^ 




(2.56) 


with i?t = (^4 — /icot l3)/{m‘^ — ), here we assume and y to be real for simplisity. The 

tl t2 

loop functions / and g are given in Appendix B.3. These corrections give the following upper 
bound on the SM-like Higgs boson mass, 


< M| cos^ 2/3 + sin^/3{ In 


2 2 
m~ m~ 

tl C 2 


mr 


+ 


{At - p cot PY 
(m? — m? )2 

tl t2 


2 - 


+ 2 (A->.cot/i) 

Cl C2 C2 

2 2 2 \ 

TTlr + m~ m~ \ \ 

~J - T ~T ) • 

m~ — m~ m~ / / 

Cl C2 t2 / ' 


If we take the universal soft SUSY breaking stop masses, 




the stop masses are written down easily as follows, 


mr 


* 1/2 


m? + =F mt{At - g cot /3). 


Then, the upper bound Eq. (2.57) becomes 


mf, < M| cos'" 2/3 + sin^ /3 ( In ^ + 


dvr 2 


rriq (At -/icot/3)2 1 (At -/r cot/3)^ 


mr 


2 

m| 


12 


4 

m| 


(2.57) 


(2.58) 


(2.59) 


(2.60) 


Therefore, the radiative corrections can give significant contributions to the mass of the SM-like 
Higgs boson. Note that the logarithmic term is given by the first and second diagram from the 
left in the Figure 2.4. As one can see, when the stop mass is heavy, this term generates large 
correction. While, the last two terms are given by the first and second diagram from the right 
in the Figure 2.4, where the bullet represents the At term interaction. A notable point is that 
these contributions are maximized by 


At — g cot /3 = ±\/ 6 mg. 

Then these contributions can be comparable to the logarithmic one. 


(2.61) 
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2.3 Current status of the MSSM 

In this section, we will review the current situation of the MSSM in terms of an observed 125 GeV 
Higgs boson and the constraints from the flavor violation and CP violation process. 

2.3.1 125 GeV 

On Wednesday, July 4 2012, the ATLAS and CMS collaborations at the LHC experiment de¬ 
clared astonishing announcements that they had observed a new particle which is consistent 
with the SM Higgs boson [1,2]. A mass of the new particle had been around 126 GeV. The 
latest measured value of the Higgs boson mass is 

my = 125.36 ± 0.37 (stat.) ±0.18 (syst.) GeV (ATLAS) [43], (2.62) 

mh = 125.03 (stat.) (syst.) GeV (CMS) [44], (2.63) 

where correspond to integrated luminosities are 4.5 fb~^ at 7 TeV and 20.3 fb~^ at 8 TeV (AT¬ 
LAS) and 5.1 fb“^ at 7 TeV and 19.7 fb“^ at 8 TeV (CMS). A naive average of the ATLAS and 
CMS results is 125.15 ± 0.25 GeV [45]. 

As we discussed in the previous section, such a Higgs boson mass can not be realized in 
the tree-level estimation of the MSSM. However, the one-loop radiative corrections can actually 
raise the mass of the SM-like Higgs boson, and so it can be realized. Therefore, considering the 
radiative corrections to the Higgs boson mass is important and essential in SUSY models. 

Although the one-loop order radiative corrections can contribute significantly to the Higgs 
boson mass, in fact two-loop order radiative corrections give a negative contribution to the one 
and it is not a negligible contribution [46-51]. It is because that the QCD corrections first appear 
in two-loop order diagrams, and they give opposite contributions to the SM-like Higgs boson 
(cf. see the last term of second line of the two-loop RGE for Aquartic Eq. (B.16)). Therefore, 
in order to predict a reliable Higgs boson mass in SUSY models, one should take the two-loop 
order radiative corrections to the Higgs boson mass into account. 

It is known that there are three ways to achieve the 125 GeV SM-like Higgs boson considering 
the two-loop order radiative corrections. First way is the heavy stop scenario. When stop masses 
are about 10 TeV, the logarithmic corrections which are generated by the stop loop become an 
appropriate magnitude. Second way is the large stop mixing scenario, namely the large At term 
scenario Eq. (2.61). In this way, even when the stop masses are about 1 TeV (and At ~ 2.5 TeV), 
the radiative corrections become an appropriate magnitude [52]. Third way is an extended 
models of the MSSM [9]. The appropriate extension models which can explain the observed 
Higgs boson mass are singlet extension models, vector-like matter extended models [53, 54], 
U(l) gauge extended models [55], etc. In the singlet extension models, an additional E-terms 
can raise the tree-level mass of the Higgs boson^^. In the vector-like matter extended models, 
the additional CP-even singlet scalar is lighter than the Higgs boson, the singlet-doublet mixing can raise 
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similarly to the stop loop, the vector-like particle loop gives the sizable radiative corrections to 
the Higgs boson mass. In the U(l) gauge extended models, an additional D-terms can raise the 
tree-level mass of the Higgs boson. Especially, we will focus on a singlet extension model in this 
thesis. 

In the following two sub sections, we will review the two-loop and higher-loop order analysis 
of the Higgs boson mass in the MSSM. It is because that we have applied this two-loop order 
analysis of the Higgs boson mass to a singlet extension model. 


2.3.1.1 Two-loop level analysis of the Higgs boson mass 


Inclusion of the two-loop level radiative corrections is important and indispensable in the calcu¬ 
lation of the mass of the Higgs boson in supersymmetric models. In this section, we hrst review 
the two-loop level calculation of the mass of the Higgs boson using the RGE [45,57,58]. Then, 
we will show a behavior of the mass of the Higgs boson as a function of SUSY breaking scale 
and tan /3 in the MSSM. 

Let us assume uew ^ -^SUSY and Mgaugino ~ /U ~ = 0{Ms\jsy) for simplicity, where 

ttiq represents dimension two soft SUSY breaking mass term of Higgs and sfermion, namely all 
the sfermions, heavy Higgs doublet A, Higgsinos and gauginos are integrated out at the scale 
Msusy- This assumption of the mass spectrum represents the effective theory below SUSY 
breaking scale Msusy to be the SM. In this section, the SM-like (surviving) Higgs doublet 
and heavy Higgs doublet are defined as 


[ \ ^ [ cos 13 sin/3 \ / -eH^ \ 

\ J V “ cos/3 J \ Hu y ’ 


(2.64) 


where e is the antisymmetric tensor ei 2 = 1. In the component representation, this equation is 
the same as follows, 

= cos/3^ ^ -hsin/3 ^ , (2.65) 

<hH = -sin^^ ^-hcos/3^^|y (2.66) 

The potential of the SM-like Higgs below SUSY breaking scale can be given by 

= (2.67) 

This potential becomes as follows when it is expanded by vacuum fluctuation of the Higgs boson, 
‘ b ® = vew + 

U(h) = 24^^Aquartic + 2UEwAquartic ' (2-68) 

the mass of the Higgs boson [56]. 
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Namely the tree-level mass of the SM Higgs scalar h is 

— 21;^^ Aquartic • (2.69) 

In order to derive the physical mass of the SM Higgs scalar h from SUSY breaking scale, we 
connect the two scale, that is SUSY breaking and electroweak scale, using the two-loop RGE for 
the Higgs quartic coupling. The full set of two-loop RGEs for the coupling constants of the SM 
using MS regularization are presented in appendix B.1.1. Then, we impose mating conditions 
of the couplings for the RGE at the SUSY breaking (high) and electroweak (weak) scale. Note 
that we estimate the mass of the Higgs boson at two-loop level (next-to-leading order), we need 
to include one-loop threshold corrections in these matching conditions. 


Matching at high scale 

Because supersymmetry ensures the relationship among the dimensionless coupling con¬ 
stants, the Higgs quartic coupling must satisfy the following matching condition at SUSY break¬ 
ing scale, 

Aquartic( Y/susy) = Alo( ^susy) = ^ (fl'^( ATsusy) + 5'^( ATsusy)) cos^ 2/3. (2.70) 

Now, in order to accurately calculate the mass of the Higgs boson, we need the matching 
condition including next-to-leading order corrections. The matching condition including one- 
loop level threshold corrections is given as follows [45,57] 

Aquartic( MsvSy) = Alo( ^SUSy) + (' 4 y,-) 2 '^NLo( YIsuSy), (2-71) 

where 


Anlo 


\reg I \<P 

^NLO ^NLO 


+ -^NLO 


+ -^NLO- 


(2.72) 


The First term in the right hand A^^o ^ convention factor from MS to DR regularization 
scheme, which gives a correction to the tree-level relation of Eq. (2.70) even in SUSY limit. 


^reg _ _ 

^NLO ~ 


1 




cos^ 2/3 
6 




(2.73) 


The other terms, A^^o’-^ nlo> "^nlo computed using the DR regularization scheme. 

The second term A^^o obtained when one integrate out the heavy Higgs multiplet and 
sfermion at the matching scale. Neglecting all Yukawa coupling except the top quark yt, this 
term is given as follows. 


'^NLO 




1 


1 


yt + ^{9^ - ^ 5 ^) cos2/3 


In- 


m. 


Q 


m.c-\ 




MsvSY 
, 3 o . 




2,^/2 0/3 

Vt + -g cos 2/3 


In- 


m 


u 


Msusy^ 






m/n 
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-ylx^cos^2^{g'^+g^)H{'^) - ^{g'^ + g^) sinH^ 

4 rriTf lb 


1 


+ Y^ [29ff'^ + 42g'^g^ + 53g^ - 4 cos 4/3(ff'^ + Qg'^g^ + 7g^) 


—9cos8/3(g('^ +5^)^] In 


m. 


Msusy^ 


+ 


cos2 2/3 


2 g'Un 


ml 


Msvsy 


m 


u 


+ -ff In , 

3 Msusy^ 


D 


,2 ,4 m 

H —q In-^ „ 

3 Msusy^ 


1 


m, 


+ 3(5"^ + 95"^) In 2 + "" + 5 '"') 2 

J iWsuSY IWSUSY 


ml 


(2.74) 


where loop functions F, G, 77, Hi,H 2 are defined in appendix B.3, and they are normalized such 
that F{1) = G(l) = .^(1) = .^ 1 ( 1 ) = .^ 2 ( 1 ) = 1- The stop mixing parameter Xt is defined by 

X, = (2.75) 

y/mqmu 

When the all masses of the heavy Higgs multiplet and sfermons are the same as SUSY breaking 
scale, A^lo becomes as follows 


a: 


NLO 


= Qyt 


X?- 


12 


+ ( 5 '^ + g^) (3 - cos 2/3) cos 2/3 


-^{g'^ + g^)smHp. 
lb 


(2.76) 


Here dominant term is the hrst one, and it reproduces the last two terms of the one-loop 
corrections Eq. (2.60). As we discussed before, this yf correction is maximized when Xt — \/6 
Eq. (2.61). Therefore, the largest threshold correction comes from the mixing of the stops, 

<LO, max ^ 18l/i • (2-77) 

Note that the contribution to the mass of the Higgs boson from this threshold correction (2.77) 
becomes smaller when SUSY breaking scale is higher. There are two reasons why the contri¬ 
bution becomes small. Eirst, the value of the top Yukawa coupling becomes smaller at high 
scale by the RGE corrections. Second, the renormalization flow of Aquartic has a focusing effect. 
Namely, in order to obtain the 125 GeV Higgs boson, Aquartic of the matching condition should 
be small when SUSY breaking scale is high. Then the dominant contribution comes from not 
threshold corrections but RGE corrections from SUSY breaking scale to weak scale. 

The third one Aj^^q is the modification term to the tree-level relation of Eq. (2.70) through 
the Higgsino-gaugino loop. 


1 


^NLO ~ 0 


25 ^ In 






-9 + (y + 9) In o 

Msusy' Msusy" 


(2.78) 


The last term A(^lo obtained as follows when one integrate out the Higgsinos and gauginos 
at the matching scale. 


\X^ 

^NLO 


1 


-7 


-^h{ri){g\d + 9iu) - ^f2{r2){9td + sL) 
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with 


3 ~ 7 ~ 8 ~ 

-2Mn)gM9lu - 2^4r2)gh9L - :^Mri,r2)hdhuhdhu 

-'^Mri,r2){9ld92d + 9iu92n) “ ^/7(n,r2)(5?d5L + 9iu92d) 

4 ~ 

— -^fs{'l"l,r2){hd92u + huhd){hdhd + 5ln52u) 

+-firi)gidgiu [Alo ~ “^{gid + 9iu)] + ‘^f{'<'2)92d92u [Alo ~ 2(^1^ + 
+ 25(7'l)ALo(ffld + glu) + 5(7’2)ALo(ffM + 92u) I 


(2.79) 


9id = 9 sin 13, g 2 d = 9 sin 13, 

giu = g'cos 13, hu = 9 cos/3, 


(2.80) 


Px — 2ALo(3ld + 9lu + ^92d + ^92u) ~ 9ld ~ 9lu ~ ^92d ~ ^92u 

— 4gidgiu92dg2u — ‘^{g'id + 92u){9lu + 92d)i (2-81) 


and ri = Mi/fi, r 2 = M 2 /iJi^^. The loop functions fi,f,g {i = 1,2, ..., 8 ) are defined in 
appendix B.3, and they are normalized such that /(I) = ^(1) = /i/ 2 / 3 / 4 (l) = / 5 / 6 / 7 / 8 (l! 1) = 1- 

y2 

When the all masses of the Higgsino and gauginos are the same as SUSY breaking scale, A(^lo 
becomes as follows 


-^NLO ~ o(l + sin2/3) 


1 


-f s'* - - 179* 


+( 39 '^ + 5^) {{ 9 '^ + 9^) cos 4/3 + 2{-g''^ + g^) sin 2/3) 


(2.82) 


Matching at weak scale 

In order to calculate the mass of the Higgs boson at the next-to-leading order, we should 
include the one-loop corrections to the tree-level mass (2.69). The pole mass of the Higgs boson 
and the top quark are related to Aquartic(<3) and yt{Q) at the MS scale Q as 


^h,pole 


TTJ-tjpole 


27’Ew(Aquartic(Q) T ^xiQ))} 
yt{Q)vEw 
1 + 6t{Q) ’ 


(2.83) 

(2.84) 


Here 6\ is the full one-loop radiative corrections via SM particle loop derived by Sirlin and 
Zucchini [59] 7 ^®. 6t is the three-loop corrections 0{a^) and the full one-loop radiative corrections 
via the SM particle loop [61,62]. 

MQ) = -Y:|ji^|^(er.(e) + r„((3) + l^), (2.85) 

7^® In the split case, vew ~ Mgaugino ~ /r <C Msusy, one must include this threshold correction in matching 
condition at not the high scale but weak scale, as we discuss later. 

7^®Recently, the full two-loop corrections to the Higgs quartic coupling are evaluated [60]. 
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HQ) = + ( 2 . 86 ) 

where Gp is the Fermi constant from the muon decay, = 246.21971 ± 0.00006 GeV, 

and ^ Mz = 91.1876 ± 0.0021 GeV. The loop function Fq, Fi and T_i are defined in 

appendix B.3. 

We have used the matching condition for the Higgs quartic coupling Eq. (2.83) at Q = mt 
with 

mt = (174.34 ± 0.37(stat) ±0.52(syst)) GeV [63], (2.87) 

asiMz) = 0.1184 ±0.0007 [64]. (2.88) 

The values of the SM couplings at Q = mt are computed by Ref. [58], using two-loop 5-flavor 
MS RGE (as initial values they use the MS mass of the bottom quark mb{mh) = 4.18 GeV and 
the pole mass of the tau lepton mr = 1.777GeV.). We use these values in matching condition 
at the weak scale, 

gs{Q = mt) = 1.1666 ±0.00314^^-^T__-o.00046(-^ - 173.35), 

(2.89) 



g{Q = mt) = 0.6483, 

(2.90) 


g’{Q = mt) = 0.3587, 

(2.91) 


yb{Q = mt) = 0.0156, 

(2.92) 


yr{Q = mt) = 0.0100. 

(2.93) 

Then, the numerical values of the next-to-leading order threshold correction Eqs. 

(2.85, 2.86) 

are 

S\{mt) ~ 0.0075Aquartic("if), 

(2.94) 


6^^^{mt) = -^as{mt) - 0.92al{mt) - 2.64a|(mt), 

OTT 

(2.95) 


dt{mt) = 5?^°(mt) ± 5f^(mt) ~ -0.0600 ± 0.0013. 

(2.96) 


We have numerically calculated the mass of the Higgs boson at two-loop level, using two-loop 
SM RGEs (appendix B.l) from SUSY breaking scale to top quark mass scale and including one- 
loop threshold corrections (Eqs. (3.62, 2.83)). In the Figure 2.5, we show the predicted mass of 
the Higgs boson as a function of SUSY breaking scale and tan j3. The green regions represent the 
appropriate Higgs mass 125 GeV < mh < 126 GeV. We take Mgaugino = h = \/~ ^SVSY for 
simplicity. The stop mixing parameter is fixed at At^Msusv) = 0 {upper figure), ^^(Msusy) = 
Msusy {middle figure) and Vt(MsusY) = VG {lower figure). Our results are consistent with the 
Figure 2 in Ref [45]. 
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We find that at large tan/3 region, an appropriate Higgs mass is achieved when SUSY 
breaking scale is 0(10) TeV. On the other hand, at small tan /3 region, it is achieved when SUSY 
breaking scale is 0(10^"^®) GeV. We also show that the stop mixing contribution Eq. (2.77), that 
is the dominant threshold correction, is effective at large tan/3 region. For example, even when 
the soft SUSY breaking scale is 0(1) TeV, the maximal stop mixing can raise the Higgs boson 
mass to be observed value (see the lower figure). However, we discussed before, this contribution 
is not effective at small tan /3 and high-scale SUSY breaking region. 

Note that in split case, that is uew rsj Mgaugino ~ At < Msusy ~ Higgsino and gaugino 

also affect the RGEs of the quartic coupling of the Higgs boson at uew < Q < Msusy in addition 
to the SM particles. In other words, the following terms is still active at uew < Q < Msusy^ 

C = -^‘h^zu2(-52dtT“lU“ + ~9idB)Hd - + huB)Hu + H.c., (2.97) 

where gid/u s-iid g 2 d/u Yukawa-like gaugino couplings. Then, since supersymmetry is no 
longer ensued at Q < Msusy j the Yukawa-like gaugino couplings are different from the corre¬ 
sponding gauge couplings'^Namely, the relations Eq. (2.80) are not satisfied at Q < Msusy- 
Therefore, in the split case, we should take into account the RGE of not only the SM couplings 
but also the Yukawa-like gaugino couplings at vew < Q < Msusy- The study of the split mass 
spectrum case is written in Refs. [45,57] in detail. 

2.3.1.2 Higher-loop radiative corrections to the Higgs boson mass 

Recently, the contributions of higher-loop correction to the mass of the Higgs mass have been 
studied. As a result, one find that these new contributions are important when SUSY breaking 
scale is not low. In this section, we briefly review the higher-loop corrections. 

Figure 2.6 shows the mass of the Higgs boson as a function of SUSY breaking scale including 
the higher-loop radiative corrections. We have numerically analyzed the mass of the Higgs boson 
using the public code FeynHiggs2.10.0 [67-74]. We take mg = mu = M 3 = g. = Ma = Msusy, 
tan/3 = 50 and At = 0. In this analysis, the mass of the top quark is 174.34GeV [63], which 
is the latest result. The green region is 125 GeV < my < 126 GeV. The dotted line represents 
the full one-loop result and the dashed line represents leading 0{y^g‘l) plus subleading 0{yf) 
two-loop result. Note that full one-loop, leading two-loop and subleading two-loop corrections 
are calculated by Feynman-diagrammatic approach, and these contributions include not only the 
logarithmic corrections but also finite corrections. The thick line represents the result of two-loop 
plus higher-loop corrections, which is evaluated by a resummation of the leading and subleasing 

We have also considered the corrections to the gaugino couplings from the corresponding gauge couplings 
[65,66] in the split mass spectrum case, in which sleptons are also light. We found that this corrections can be 
as large as ~ 10% in the parameter region in which the muon g — 2 anomaly can be solved, and that the gaugino 
couplings can be measured from the production cross section of the right-handed selectrons at 1% accuracy at 
ILC with y/s = 500 GeV. 
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Figure 2.5: Two-loop level mass of the Higgs boson as a function of SUSY breaking scale and 
tan/3. The result is evaluated by using two-loop SM RGEs from SUSY breaking scale to top 
quark mass scale and including one-loop threshold corrections. The green regions represent the 
appropriate Higgs mass 125 GeV < rrih < 126 GeV. We take Mgaugino = Ai = = VfsuSY- 

The stop mixing parameter is fixed at Hi(MsusY) = 0 {upper figure), Hi(MsusY) = YAsusy 
{middle figure) and ^^(Msusy) = {lower figure). 
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^susY [TeV] 


Figure 2.6: The mass of the Higgs boson as a function of SUSY breaking scale including the 
higher-loop radiative corrections. We take mg = mu = M 3 = // = Ma = Msusy, tan/3 = 50, all 
H = 0, and mt = 174.34 GeV [63]. The dotted (dashed) line represents the full one-loop (leading 
plus subleading two-loop) result. The thick line represents the result of two-loop plus higher- 
loop corrections, which is evaluated by a resummation of the leading and subleasing logarithmic 
corrections from the scalar top sector. The blue and red regions represent one sigma bands from 
the theoretical uncertainty. The green region is 125 GeV < my < 126 GeV. 


logarithmic corrections from the scalar top sector. The resummation have been obtained from 
an analysis of the RGE at two-loop level. The blue and red regions represent one sigma bands 
from the theoretical uncertainty. The theoretical uncertainty is dominated by two contribution, 
the experimental error on the mass of the top quark and unknown higher-order corrections. The 
theoretical uncertainty from the experimental error on the mass of top quark is numerically 
estimated [75], and it is obtained as follows, 

6m^^ ~ P ~ 1 GeV. (2.98) 

On the other hand, the theoretical uncertainty from the unknown higher-order corrections is 
estimated at ~ 3 - 5 GeV in two-loop level calculations. However, as you can see 

Figure 2.6, including the resummation of the leading and subleading logarithmic corrections, 
this uncertainty dramatically decrease [74,76] at 

^^higher <^QeV. (2.99) 
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2.3 Current status of the MSSM 


In this Figure, we show that the higher-order corrections by the resummation of the leading 
and subleading logarithmic corrections can raise the mass of the Higgs boson. We also find that 
when there is no help of At term, SUSY breaking scale is predicted Msusy ~ 5TeV in order to 
obtain an appropriate mass of the Higgs boson (tan/3 = 50 case). 


2.3.2 SUSY FCNC / CP Problem 


In general, the supersymmetry introduces the new flavor and CP violating sources through 
the SUSY breaking sector. These new contributions seem to cause the flavor changing neutral 
currents (FCNC) and electric dipole mordents (EDM) of quarks and leptons. On the other hand, 
in the SM, the FCNC is suppressed by the GIM mechanism [77], and the EDM is suppressed 
by the three (four) loop suppression and by the smallness of component of the CP violating 
phase in the CKM matrix. Hence, the FCNC and EDM which are generated from the SUSY 
breaking sector can give dominant contributions, and they can be probed by the low energy 
flavor experiments and the EDM measurements. Nevertheless, the corresponding signals have 
not been observed yet. The current experimental results set the constraint to the new flavor 
and CP violating sources in the supersymmetric models. 

In the MSSM, the new flavor violations are generated by the off-diagonal sfermion soft 
SUSY breaking terms (i ^ j) and the off-diagonal A terms Atj (i ^ j). If the A terms 

are proportional to the Yukawa matrix, the main sources of the flavor violation are given by 
the sfermion masses. Besides these off-diagonal components are expected the same order as the 
diagonal components if there are no additional symmetries'^®. Thus, the supersymmetry would 
cause relatively large FCNC process. 

One of the severe constraints comes from the branching ratio of // —)• ey. In the MSSM, when 
we assume (m| )ij = and neglect the A term contributions, the branching ratio 

of /r —)• ey at the one-loop order is given as [65,78], 


T3 / \ ^ OtTfl^ 

Br(/i —ey) ~ 


tot 


16 


("lf)23("lf)31 m. 


m~ 


m. 


bino(neutralino) loop 


+ 




mt 


chargino loop T wino-Higgsino(neutralino) loop 


T bino-Higgsino(neutralino) loop 


( 2 . 100 ) 


with 


m: 


chargino loop — 


IOtt^ M2fJ, 


tan 13 fc 


M| 




2 ’ 2 


( 2 . 101 ) 


large off-diagonal components cause the negative eigenvalue of sfermion mass square, which leads to the 
charged/colored breaking vacuum. Thus, actually the off-diagonal components are expected 0(0.1 x (mpii)- 
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wino-Higgsino(neutralino) loop 


bino-Higgsino(neutralino) loop 


blno(neutralino) loop 


327r2 M 2 IL 


327r2 Min 


tan/3/Ar 


tan Pfjsf 


(^ 

[my ml) 
[my ml) 


g'"^ mlMin 


tan 13 fN 


(^ rnl \ 

[my Ml) ’ 


( 2 . 102 ) 

(2.103) 

(2.104) 


where (m| )22 = m^ the total decay width of muon Ftot = 2.99 x 10”^® GeV, and the loop 
functions fc and /w are defined in Appendix B.3. Note that the SUSY contributions are given 
by the chargino-sneutrino loop and the neutralino-smuon loop. As one can see, these effects are 
proportional to tan^ /3. Naturally, these SUSY contributions are decoupled in the Msusy 00. 

The current bound on the /i —)• ey had been given by MEG Collaboration, and the result 
is Br(// ^ ey) < 5.7 x 10“^^ (90 % CL) [79]. The Ref. [80] showed that when ■sjml = /r = 
Msusy, 3Mi/{5g'y = M 2 /g^ = Msusy/ gl, tan/3 = 50, (m|) 2 i/m? = 0.1 and (m|) 23 /m? = 
(^Tipsi/m? = 0 are taken, the current bound implies that the soft SUSY breaking scale should 
be heavier than 30TeV, Msusy ^ 30TeV. Thus, this bound severely constrains the low scale 
SUSY models if there are no additional symmetries or mechanisms. Note that this bound can 
be relaxed if we take small tan/3. 


One the other hand, in the MSSM the new CP violation sources are the off-diagonal sfermion 
soft SUSY breaking masses Im((mj)jj) {i / j), A terms Im(Ajj), the gaugino mass Im(Mgaugino) 
and the Higgsino mass Im(/i). The physical CP violation source is relative phases of those. 
Although the quark (lepton) EDM is generated by three (four) loop diagram in the SM, the EDM 
is actually generated by one loop diagram in the SUSY models. Therefore, the supersymmetry 
also would cause relatively large EDM. 

One of the severe constraints to the CP violation phenomena comes from the electron EDM. 
In the MSSM, when we take the same assumption as the above, the electron EDM at the one-loop 
order is given as [65], 


4 

e 




2mp 


arg[M2/x] (oe, chargino loop 3~ wino-Higgsino(neutralino) loop) 


-F arg[Mi//] (oe, 

bino-Higgsino(neutralino) loop T ®e, bino(neutralmo) loop) 


,(2.105) 


where we have replaced subscript /x with e in the Eqs. (2.101-2.104). Similarly to the /x —)• ey 
process, the SUSY contributions are given by the chargino-sneutrino loop and the neutralino- 
smuon loop. These contributions are proportional to the CP violating relative phase of the 
Higgsino and the gaugino masses and proportional to tan /3. When the all dimensional parameter 
are the same values Mi = M 2 = /x = ^ ^ m| = Msusy, and if we take the universal 
CP violating phase of gaugino mass arg(Mi/x) = arg(M 2 /x) = sin^, the one loop order electron 
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EDM is given as, 


\de 


+ g‘ 


n 


rfip 


■tan/3sin(/) x 197 x 10 [e cm] 


2 • 487r • 47r Msusy' 


= 2.9 X 10 


_27/10TeV\ /tan/ 3 \ . . , 

-^7 - sin 0 e cm] 

V Msvsy y V 50 y 


e cm . 


(2.106) 


The current bound on the electron EDM had been given by ACME Collaboration, and 
the result is |de| < 8.7 x 10“^® [e cm] (90 % CL) [81]. The Ref. [80] showed that when 
\f^e = h = Msvsy, 3Mi/(5c/'2) = M2/52 = Msusy/Gs, tan/3 = 50, {mj)ijlml = 0 (i / j) 
and the maximal CP violating phase are taken, the current bound implies that the soft SUSY 
breaking scale should be heavier than 100 TeV, Mgusv ^ 100 TeV. Although this bound can 
be relaxed if one takes small tan/3, it also severely constrains the low scale SUSY models. 

Another severe constraint comes from measurement of the CP violation in the kaon decay [82] . 
Under the CP transformation a neutral kaon (ds) becomes an anti-neutral kaon {sd), 
and it is represented by convention. 


CP|A°) = +\K^), 

CP|A°) = +\K°)- (2.107) 


Hence, (|iL°) -|- \K^))/V2 and (|Ar°) — \K^))/\/2 states are the CP even and odd eigenstates. 
Since the CP symmetry is slightly broken in nature, the mass eigenstates of the neutral kaons 
become the following forms. 



1 

71 

1 

71 


e-*<^x(l + eK)\K^) + - eK)\K^)\ , 

e-*'Ax(i + eK)\K^) - e'‘^^(l - eK)\K^)\ , 


(2.108) 


where ex is an indirect CP violation parameter (CP violation in K^-K^ mixing) and (pK is a 
direct CP violation parameter (CP violation in the neutral kaon decay), and actually eK ^ 
1) 4 ’k ^ 1- The indirect CP violation parameter ex is determined by an imaginary part of the 
K^-K^ mixing [83,84], 


ex = sin 


Im(A47^o_^o) 

AMx 



with 


(2.109) 


Mxo.xo = {K^\Ks 

2AMx 


tan d>f = 


AFx 


{dsds)\K'^ 

= (43.52 ±0.05)" 


( 2 . 110 ) 


where (j^xfi is the phase of the 0-isospin amplitude in kaon decay —>• 27r, and 


AMx = Ml - Ms = (3.484 ± 0.006) x 10"^^ MeV, 
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ATk = Ts-r^ = (7.2823 ±0.0098) X MeV. (2.111) 

The second term of the Eq. (2.109) contributes an 0(5%) correction to €k [84]. 

In the SM, the leading contribution of the four-fermion Hamiltonian 77gg^^™((Js(is) comes 
from one-loop box diagrams with the weak interactions and the imaginary component of the 
CKM matrix. While in the MSSM, the leading SUSY contribution comes from one-loop gluino- 
squark box diagrams with the strong interaction, and these contributions are easy to affect 
estimation of ck- The main CP-violating source of this diagram is the off-diagonal sfermion soft 
SUSY breaking masses Im[(m^ )i 2 ]. The detailed formulae of the SUSY contributions 

to ex are given in Ref. [85]. 

The latest experimental value is = (2.228 ± 0.011) x 10“^ [ 86 ] and the SM prediction 

(including next-to-next-to-leading-order) is = (1.81 ±0.28) x 10“^ [87]. These numerical 

values set a conservative upper bound of the SUSY contributions, < 0.98 x 10“^ [80]. 

The Ref. [80] showed that when = M 3 = MsusYj ~ 7^ 

j) and the maximal CP violating phase are taken, the current bound implies that the soft SUSY 

breaking scale should be heavier than 500 TeV, Msusy ^ 500TeV. If one imposes the SO(IO) 

relation at the GUT scale; (m^ )ij = (m^ )*j, the imaginary components of the soft SUSY 

breaking squark masses vanish at the GUT scale and they are generated only through the RG 
fSUSY) 

effects, and thus e]^ ' is suppressed. Then, the current bound implies Msusy ^ 40 TeV. 

Therefore, when the soft SUSY breaking scale is low scale which is favored in terms of the 
naturalness, too large FCNC and too large CP-violation are induced by the SUSY particles one- 
loop radiative corrections. These difficulties are called “the SUSY FCNC problem” and “the 
SUSY CP problem”. 

2.4 Discussions 

In previous section, we have reviewed the current situation of the MSSM in terms of an observed 
125 GeV Higgs boson and the constraints from the flavor violation and GP violation process. 
Actually, we have shown the following two points in the MSSM, 

• In order to realize the observed Higgs boson mass, the masses of the stops have to be heavier 
than 5 TeV or the stop mixing has to be maximized (or their compromised parameter 
region). 

• In order to avoid the severe constraints from the flavor violation and CP violation mea¬ 
surements, the masses of the sleptons should be heavy and the parameter tan (3 should be 
not so large if there are no additional symmetries or mechanisms. 

Therefore, in this chapter we naively conclude that one of the suitable solutions is the high-scale 
supersymmetry, which has 0(10-100) TeV soft SUSY breaking terms. 
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2.4 Discussions 


Once we give up thought on the naturalness, such a high-scale supersymmetric model is very 
attractive. First, as we discussed, the observed 125 GeV Higgs boson and the explanation of 
the current bound on the flavor violation and CP violation process are simultaneous realized 
naturally. Second, all the SUSY particles are out of the reach of the current LHC, thus it can 
explain that there are no signals of the SUSY particles at the LHC experiments. Third, it can 
naturally solve the cosmological gravitino problem [88]: although an unstable gravitino spoils 
the big-bang nucleosynthesis, the heavy gravitino (typically heavier than 5-8 TeV) can avoid the 
disaster"^®. Fourth, this heavy SUSY breaking scale remains the gauge couplings unification at 
the GUT scale (see Figure 2.2). 

Furthermore, the future prospects of the flavor violation and the CP violation measurements 
have the potential to probe the supersymmetry beyond the reach of the LHC by orders of 
magnitude [89-96]. Therefore, we hope that these precision measurements can indirectly probe 
the high-scale supersymmetry. 


#'^Here, we implicitly assumed that the soft SUSY breaking scale and the gravitino mass are the same scale. 
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Chapter 


3 

Singlet Extension 


In this chapter, in order to solve the /r problem of the supersymmetric min¬ 
imal model, we first introduce a gauge singlet superfield and singlet extension 
supersymmetric minimal model. Then, we should impose extra symmetries to 
forbid unwanted terms of singlet superfield, which spoil the solution of the /r 
problem. In general, these symmetries lead to the domain wall problem and 
the tadpole problem. The nearly minimal supersymmetric standard model is 
the one of the models which can solve the ^ problem, the domain wall problem 
and the tadpole problem simultaneously. Finally, we review the Lagrangian 
and the Higgs sector of the nearly minimal supersymmetric standard model, 
and we show that there is a sizable tree-level contribution to the Higgs boson 
mass due to an extra F-term contribution to the Higgs quartic coupling. 

3.1 ji Problem 


In the MSSM, the potential minimization conditions Eqs. (2.41, 2.42) lead to the following 
relationship, 


M| = 


ms — , ^9 mS — m 

2 1 - 2 


- — nil — 1712 — 


cos 2/3 cos2/3 

where /r is the supersymmetric mass of the Higgs multiplets. Especially, when tan/3 /:$> 1, this 
equation can be expanded by tan /3 and becomes 

2.9 9. 1 . 


— —2(m2 + |/x|^) -b 


tan^ /3 


{mf — m^) + 0{ 


tan^ /3 ^ 


(3.2) 


These relationships imply that since actually Mz ~ 0(100) GeV (electroweak scale), mf, sj 
and [I should be naively at the electroweak scale, or ^^/vn^ and should be the same scale 
to able to cancel out. In other words, the magnitude of /i have to be the soft SUSY breaking 
scale Msusy or be less than the scale. 

One the other hand, the /z term is stable under the all orders in perturbation theory of 
the effective Lagrangian due to the non-renormalization theorem [97], and the supersymmetry 
provides a valid description that the scale of the /z parameter is as large as GUT scale or Planck 
scale^^. If /z is at the Planck scale, it leads to Mz ~ Mpi (see Eq. (3.1)), and the appropriate 

Another elegant possibility is /z = 0, which is respected the some symmetries (e.g. Peccei-Quinn symmetry). 
However, nature have not chosen this values. It is because that the lightest chargino search by the LEP set lower 
bound on chargino mass, |zt| > 103.5 GeV [98]. 
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Table 3.1: The singlet superfield with their components for spin 0 and 1/2, and their represen¬ 
tations for SU(3)cX SU(2)iX U(l)y gauge group. 


Chiral Supermultiplet 

Spin 0 Spin 4 

SU{3)c SU{2 )l U{1)y 

Singlet scalars-Singlino S 

S = {SR + iSi)/V2 S 

110 


electroweak symmetry breaking can not occur. Therefore, pL has to know the soft SUSY breaking 
scale in order to realize the nature, namely the Z boson mass is the electroweak scale. This 
problem is called problem” [8]: 

Why pL <C MguX) .Tfpi ? 

Why does /i know Msusy ? (3-3) 


One of the natural solutions of the /i problem is the singlet extension models of the MSSM [9]. 
These models have the following superpotential, 

W = XSH 2 H 1 + /[5] + WYukawa, (3.4) 

where S is an additional gauge singlet superfield, A is a dimensionless coupling constant, and f[S] 
is the superpotential which does not depend on superfields of the MSSM at the renormalizable 
level. The singlet superheld S with their components for spin 0 and 1/2 are shown in Table 3.1. 
When supersymmetry is broken, singlet superfiled also receives the soft SUSY breaking mass or 
A term. At this time, singlet scalar boson can naturally obtain a vev which is the order of the 
soft SUSY breaking scale, (S) ~ 0( Msusy), and so its vev gives the effective /r term for Higgs 
multiplets, 


IJ-eS = HS). 


(3.5) 


Therefore, the fi problem can be solved by the singlet extension models of the MSSM 

One of the other solutions of the p, problem is provided by the Giudice and Masiero [99]. Let us consider the 
following Kahler potential, 

K = ci^X^H2Hi+C2H2Hi+'R.c., (3.6) 

M-pi 

where ci and C 2 are 0(1) couplings and X is an additional gauge singlet superheld. Note that the second term is 
permitted in the supergravity theory. This Kahler potential can be rewritten to the following superpotential, 


IV = 



^C 2 + Cl 


(Y*) A 
Mpi J 



H2H1, 


(3.7) 


where 1713/2 is the gravitino mass which is equivalent to a vev of the superpotential. Therefore one can obtain the 
appropriate effective p term, 


Peff = Cl 


Ft 

Mpi 


+ C 2 + Cl 


(Y*) 

Mpi 


7713/2 ~ 0 { Msusy). 


(3.8) 
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3.1 Problem 


In the case that the additional superfiled for the MSSM is only one singlet superfield, the 
superpotential f[S] can be written as 

f[S]=CFS+^fr'S^ + ^KS^. (3.9) 

Then, the soft SUSY breaking terms are 


Usoft = 771111^11^+77121^^21^+^1151^ 

+ (xA^SH 2 Hi + ^sS + + H.c. 


(3.10) 


here \/m\, y w?g, Ax, (Cs)^^^; y SUSY breaking terms, and their 

magnitudes are typically Msusy- The minimization equations of the singlet scalar potential^^ 
lead to the vev of the singlet scalar boson. In other words, the vev of singlet scalar (S) is a 
solution the following equation, 

- A7;i7;2(^a + +') 


(5) = - 


mg + m'g + + 2k^p + kAk{S) + 2k?{S)‘^ + 3Kfj,'{S) + + v^) — 2 XkviV2 


AIsvsy^ + 


2Msusy^ + + 2k^p + K Msusy (5) + 2k? [3)“^ + 3k^'{S) 


(3.11) 


here we have neglected the terms which includes the Higgs vev since vi,V 2 <C Msusy■ Because 
the solution of this equation can not be estimated intuitively, let us consider the case of k <C 1 . 
The equation becomes 


A'hlJSY^ + 
2 Msusy ^ 


(3.12) 


When ^p < Msusy^ and p' < Msusy j then the singlet scalar can obtain an appropriate vev, 
(5) 0( Msusy)- However, since ^p and /r' are the dimensional supersymmetric tadpole and 

mass term, they have no reason that ^p < Msusy^ and /r' < Msusy j and their magnitudes are 
typically the GUT scale or the Planck scale if there are no symmetry. Then, the singlet scalar 
obtains an inadequate vev for the solution of the // problem, (5) ~ O(Mpi). 

These facts imply that in order to solve the p problem we have to impose some symmetries, 
which suppress or forbid some terms in the superpotential /[.S']. For example, when some 
symmetries forbid the dimensional supersymmetric tadpole ^p and mass term p', the singlet 
scalar obtains the following vev [100], 


(S) ~ 


1 

Ik 




(3.13) 


Thus, the singlet scalar can obtain an appropriate vev, (5) ~ O(MgusY)- In this manner, 
there are various singlet extension models which are classified by the symmetries to control the 

Strictly speaking, we should also solve a condition for an absolute minimum of the scalar potential. 
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Table 3.2: The charge assignments under the Abelian symmetries which are considered in he 
text. The index i denotes the generations. 



Hi 

H2 

s 

Qi 

u^ 

Di 

Li 

E^ 

W 

U(l)y gauge 

-1/2 

1/2 

0 

1/6 

- 2/3 

1/3 

-1/2 

1 

0 

Z3 C U(l)pQ 

1 

1 

-2 

-1 

0 

0 

-1 

0 

0 


0 

0 

2 

1 

1 

1 

1 

1 

2 

C U(l)p, 

1 

1 

4 

2 

3 

3 

2 

3 

6 (1 (mod 5)) 


superpotential f[S]. Note that, when f[S] = 0, there are extra global U(l) symmetries in the 
Lagrangian. This symmetry leads to an unwanted visible Nambu-Goldstone boson or a visible 
axion, when the symmetry is spontaneously broken that is associated the electroweak symmetry 
breaking. 

What symmetry is useful for controlling the superpotential f[S]7 

The imposition of extra global symmetries is unreasonable. Because when the singlet scalar 
boson obtains vev, this symmetry is spontaneously broken and gives an unwanted visible Nambu- 
Goldstone boson or a visible axion. On the other hands, discrete symmetries are suitable in 
order to control the singlet superfield [101-103]. It is because that even if this symmetry is 
spontaneously broken, all extra scalar bosons can have a heavy mass which is phenomenologically 
acceptable. Furthermore, discrete R-symmetries could be more suitable. The R-symmetry can 
be imposed in the supersymmetric models, and the R-symmetry breaking is related to the 
supersymmetry breaking. Therefore, we are able to easily estimate and control terms in the 
superpotential which are generated via the R-symmetry breaking"^^. 

3.2 Domain Wall Problem / Tadpole Problem 

In order to control the singlet superfield, one can impose a discrete symmetry. However, the 
spontaneous broken of the discrete symmetry, which is caused by the electroweak symmetry 
breaking, leads to a disastrous cosmological domain walls. Furthermore, if we introduce the 
explicit breaking terms of the discrete symmetry to avoid the domain walls, there are cases 
where these terms lead to a disastrous tadpole. In this section, we briefly review these two 
problems with the desecrate Z3 symmetry as an example. 

The Next-to-Minimal Supersymmetric Standard Model (NMSSM) [9,105-111] is one of the 

Moreover, such a desecrate R-symmetry should be embedded in a gauge symmetry. In other wards, the 
discrete R-symmetry should be anomaly-free. Otherwise since the R-symmetry is anomalous broken, one can 
not estimate the terms which break the classical R-symmetry. Therefore, the anomaly-free discrete R-symmetry 
should be employed. Actually, the discrete Z 5 * R-symmetry can be anomaly-free if one introduces the extra 
charged singlet superfields to the nearly minimal supersymmetric standard model [104]. 
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singlet extension models, and its matter content is the MSSM matter and an additional gauge 
singlet superfield. In order to control the singlet superfield, the NMSSM is imposed the desecrate 
Z 3 symmetry. In this symmetry the superhelds are transformed as 

/V . 27rq^ y.. 

^ e*—(3.14) 

where the charge assignments qi is listed in Figure 3.2. Then, the NMSSM has the following 
superpotential, 

W = \SH2Hi + ^kS^ + Wyu^^^. (3.15) 

Thus, the superpotential of the NMSSM does not have the dimensional supersymmetric tadpole 
and mass term. The vev of the singlet scalar is given Eq. (3.13), and the p problem is solved. 
Note that without the singlet cubic term , this superpotential becomes invariant under an 
anomalous global U(1)pq symmetry, which includes the desecrate Z3 symmetry as a subgroup. 
When the scalar bosons obtain its vev, this symmetry is spontaneously broken and gives an 
unwanted visible axion [112]. Thus, one needs this singlet cubic term as discussed previous 
section. 


Domain wall problem 

This desecrate Z3 symmetry has to be spontaneously broken in order to give the appropriate 
electroweak scale. The discrete symmetry is broken in different ways in different domains, which 
are separated in a larger distance than the horizon size or correlation length. Eventually, these 
different domains are divided by the domain walls (domain boundaries). If the discrete symmetry 
is exact symmetry, these domain wall configurations are topologically stable. 

When the domain walls had not disappeared in the early universe, they would produce 
destruction of the observed homogeneity and isotropy of our universe. In addition, in the 
radiation dominated era, the energy density contributions to the universe by the domain walls 
can become comparable to the energy density of the universe, and so the existence of the domain 
walls would change the evolution of the universe significantly. Therefore, one should avoid the 
disastrous cosmological domain walls. This difficulty is called “domain wall problem”. 

One of the solutions of the domain wall problem is an addition of tiny explicit discrete 
symmetry breaking terms. These tiny explicit breaking terms remove the vacuum degeneracy 
of the different domains. It can be interpreted in the decay of the domain walls, and eventually 
the universe is covered by a unique vacuum. 

Let us consider the NMSSM case, that is desecrate Z3 symmetry and its explicit breaking 
terms. The additional dimension-5 operators, which are the Planck suppressed explicit Z3 
symmetry breaking terms, are given as, 

Tipi ’ T/pi ’ Tfpi 
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in the superpotential Eq. (3.15). Note that the additional dimension-5 operators in the Kahler 
potential, {S + S'^){HiH\)/M-p\ and + H.c.)/Mpi, can be absorb into the last two terms 

of the superpotential Eq. (3.16) by the redefinitions of the superfields. These explicit breaking 
terms give the pressure to the domain walls. The strongest constrains come from the big-bang 
nucleosynthesis. It is because that the decay of the domain walls produce the entropy, and it 
destroys the great success of the nucleosynthesis. Thus, we require that the domain walls have 
to decay before the onset of the nucleosynthesis, and it sets the following lower bound on the 
parameter A' [11], 

A' > 10-". (3.17) 


Tadpole problem 

However, such a model, which is no longer protected by the exact discrete symmetry, permits 
the renormalizable explicit discrete symmetry breaking terms in the superpotential at the loop 
level. Especially, there is the quadratic divergence in the tadpole, which is generated via the loop 
diagram including the explicit symmetry breaking non-renormalizable terms. The cancelation 
of the quadratic divergences of the two-point functions remains even when the supersymmetry 
is softly broken. However, in fact, when the supersymmetry is softly broken, the quadratic 
divergences of the tadpole are not canceled out [113]. It is because that the coefficients of the 
quadratic divergence have dimension one, and so they depend on the soft SUSY breaking masses 
and A terms. Thus, the supersymmetry can not keep the cancelation of the quadratic tadpole 
divergences due to the soft SUSY breaking terms. It means reintroduction of the hierarchy 
problem. These facts lead too large tadpole of the singlet superfield at the loop level, and it 
gives too large vev of the singlet scalar terms Eq. (3.11). Therefore, the singlet scalar no longer 
solve the // problem. This difficulty is called “tadpole problem”. 

Let us consider the NMSSM which includes the explicit symmetry breaking terms Eq. (3.16). 
The leading quadratic tadpole divergences appear from the two-loop diagrams. These diagrams 
are shown in the Figure 4 in Ref. [11]. These contributions to the Lagrangian are 




^ ^(5 + 5*)A2m| + 


A' 


(167r2)2 Mpi 3 
1 A' K 


(167r2)2 M-p\ 


\{S + S*)A^m% 


+ 


■-{Fs + F*s)K^A,+ 


A' 


-\{Fs + F*s)k^Ax, 


(167r2)2 Mpi 3' (167r2)2Mpi' 

-h A^ Mpi Msvsy'^S + f o + -^ ) ^P1 MsvsyF s 


(167r2)2 V3 
+H.C., 


(Ibvr 


(3.18) 


where we have taken the UV cut off A to be the Planck scale. The first term is regarded as 
the tadpole in the soft SUSY breaking scalar potential, and the second term is regarded as the 
tadpole in the (effective) superpotential. Thus the tadpole is generated by the loop diagram 
which has the soft SUSY breaking effect. 
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Including these tadpoles in the scalar potential, the vev of the singlet scalar becomes too 
large (see Eq. (3.11)). The demand that the scale of the effective p term should be the soft 
SUSY breaking scale in terms of the solution of the p problem sets the following upper bound 
on the parameter X' [11], 

A'<10-“. (3.19) 

In this estimation, A and k are assumed 0(1). 

As one can see, two demands for the solution of the domain wall problem Eq. (3.17) and 
for the solution of the tadpole problem Eq. (3.19) are inconsistent obviously. Therefore, the 
NMSSM with Z3 symmetry, which is explicitly broken by the Planck suppressed operators, can 
not solve the domain wall problem and tadpole problem simultaneously [11,114-116]^^. 

3.3 Solution of /i Problem, Domain Wall Problem and Tadpole 
Problem 

Refs. [13-16] have shown that the desecrate Z5 R-symmetry, which can control the superpotential 
f[S] and the singlet scalar vev, can avoid the the domain wall problem and tadpole problem 
simultaneously. So, we review this symmetry in the following. 

When the superpotential f[S] = 0, apart from ordinary Lepton and Baryon number symme¬ 
tries, the Lagrangian has two additional global continuous symmetries. That is the anomalous 
Peccei-Quinn symmetry U(l)pg and a non-anomalous R-symmetry U(l)p. The charge assign¬ 
ments for these symmetries is given in Table. 3.2. Therefore, the Lagrangian is also invariant 
under global U(l)p/ transformation, where charges of U(l)p' symmetry are defined as 

R' = 3R + PQ. (3.20) 

In order to avoid an unwanted visible Nambu-Goldstone boson related with the spontaneous 
global U(l)p' symmetry breaking, let us introduce the discrete symmetry. Actually, the maximal 
discrete sub-symmetries of U(l)p/ is discrete R-symmetry. In this symmetry the superfields 
are transformed as 

/V .27rqA yv 

W e^^W, (3.21) 

where the charge assignments qi is listed in Figure 3.2. 

When one imposes this discrete R-symmetry on both the Kahler potential and the superpo¬ 
tential, the Lagrangian becomes the desired form up to a possible singlet tadpole term, which 

Hamaguchi, Nakayama and Yokozaki have pointed out that the NMSSM in gauge mediation SUSY breaking 
with vector-like exotic matters, which are charged under the hidden QCD, can solve the domain wall and tadpole 
problem [117,118]. In this model, Z 3 symmetry is anomalous, so that the domain wall problem can solve [119]. 
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is generated from the non-renormalizable sector. The discrete symmetry is imposed up to the 
higher dimensional Kahler potential and superpotential completely. This is the different point 
to the previous section. Then, the desecrate R-symmetry is spontaneously broken by R- 
symmetry breaking sector whose charge is one. This R-symmetry breaking sector 

generates the vev of the superpotential Wo, which is required for small cosmological constant. 
In addition, Wq is related with the gravitino mass as follows, 

= ^0 = ™'3/2-^pi- (3.22) 

The tadpole is generated by the following higher dimensional interaction with the R-symmetry 
breaking sector, 

P'23) 

where ci and C2 are dimensionless 0(1) coupling constants. Because of an imposition of the 
discrete R-symmetry, a harmful quadratic tadpole divergences does not appear from loop dia¬ 
grams. Once the discrete R-symmetry is broken by Wq, the tadpoles of the singlet fields are 
induced, 

Aad ~ Msvsy^S + Msusy'^Fs + H.C., (3.24) 

where we assume mo /2 ~ -^SUSY- As we will show next section explicitly, the singlet scalar 
obtains the following appropriate vev (see Eq. (3.11)), 

(S) ~-^ ~ 0 {Msvsy), (3.25) 

with 

ts ~ Msusy^- (3.26) 

Therefore, the generated tadpole and effective fi term do not destabilize the hierarchy, and this 
symmetry can naturally solve the /i problem and tadpole problem. 

Furthermore, the existence of the tadpole term, which is effectively generated by the R- 
symmetry breaking, breaks the discrete symmetry explicitly. Thus, this symmetry can naturally 
avoid the domain wall problem [14]. 

Hence, the singlet extension model imposed the desecrate R-symmetry can naturally 
solve three problems simultaneously: the /r problem, the domain wall problem and the tadpole 
problem. This singlet extension model is called “nearZy minimal supersymmetric standard model 
(nMSSM) [13-15]. 

In addition, the desecrate Z^ R-symmetry prohibits the dangerous D < 5 Baryon or Lepton 
violating operators like QQQL and UUDE in the nMSSM matter contents [14]. Thus, the 
constraint of the proton decay is satisfied. 

Next section, we review the Lagrangian and the Higgs sector of the nMSSM. 
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3.4 Nearly MSSM 

In this section, we review the Lagrangian, the Higgs sector and the Landau pole constraint of 
the nearly minimal supersymmetric standard model [13-15]. 

3.4.1 Lagrangian 

In the nMSSM, the superpotential is given as 

W = ASB 2 B 1 + + HVukawa , (3.27) 

where HVukawa is defined as Eq. (2.29). The soft SUSY breaking terms are given as 

= m?|17i|2 + mi|H2|2 + m||5|2 + (AHAH2Hi5 + ts5 + H.c.) 

Th^oft gaugino T l^oft Yukawa j (3.28) 

where I4oft gaugino and I4oft Yukawa are defined as Eqs. (2.31, 2.32). As discussed previous section, 
although the terms m^2 and ts are forbidden by the discrete R-symmetry, when the R- 
symmetry is broken they are generated. Let us parameterize these tadpole terms as follows. 


m 


12 


= Acf Msusy^ , 


ts = cs Msusy^ , 


(3.29) 

(3.30) 


where Msusy denotes the SUSY breaking scale. Here, cp and cs are 0(1) complex constants 
and then mf 2 and ts become O(Msusy^) and O(Msusy^) respectively^®. With these values, S 
has a vacuum expectation value (S) ~ —ts/rng ~ ^( -^susy) as we will discuss next subsection. 
Thus the generated effective term is 0( Msusy) and the ^ problem can be solved. 

Note that this model is imposed the desecrate Z^ R-symmetry, which is broken by the R- 
symmetry breaking terms in the hidden sector. Then, as a low-scale (TeV scale) effective theory, 
the nMSSM Lagrangian Eqs. (3.27, 3.28) are generated. It is known that, however, there are 
other models which generates nMSSM Lagrangian Eqs. (3.27, 3.28) as a low-scale (TeV scale) 
effective theory. An example is Peccei-Quinn invariant NMSSM [120,121]. The Lagrangian of 
this model is given as. 


C = [ (feXSH2Hi + [ d^e^(X^^S + R.c.), 
J J Mpi 


(3.31) 


where X is an axion superfield. This Lagrangian is controlled by U(l) Peccei-Quinn symmetry 
for (Hi, H 2 , S, X) carrying the charges as (1, 1, -2, -1). The vev of the scalar component of 
X becomes the axion decay constant fa ~ 10^®-10^^ GeV, and the vev of the auxiliary field Fx 
is zero. Then, similarly to the Giudice Masiero mechanism, by the supergravity interaction the 


Although the trilinear kS^ term is also generated, it is highly suppressed by Planck scale. 
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holomorphic term in the Kahler potential can be rewritten to the superpotential. As a result, 
the (TeV scale) effective superpotential is given as 

Weff = \SH 2 H 1 + (3.32) 

Aipi 

As one can see, the tadpole is generated and its dimensional coupling is actuary C>(Msusy^)- 
Therefore, this model can solve not only /r problem, the tadpole problem but also the strong 
CP problem, and effective Lagrangian can be regarded as the nMSSM. 

Other example is secluded U(l)'-extended minimal supersymmetric standard model (sMSSM) 
[122]. The sMSSM contains a U(l)' gauge symmetry, Z' gauge boson, and four additional sin¬ 
glets. The superpotential is 

W = ASH2ffi + ^sSiS2S3 + fiiSSi +/ 12 SS 2 , (3.33) 

and this additional gauge symmetry is motivated by the GUT. The U(l)' charges satisfy qjj-^ + 
QH 2 + 95 = 0 and —qs = qsi = QS 2 = “1/2953- If the dimensional couplings /i* is controlled by 
some mechanisms, and three singlet Si masses are heavier than soft SUSY breaking scale, the 
effective Lagrangian can be regarded as the nMSSM, 

lUeff = XSH 2 H 1 + (m(«5i) + MS 2 ))S. (3.34) 

Other example is Fat Higgs model [123]. In the Fat Higgs model, the Higgs doublet fields are 
composite bound states of fundamental fields, which couple to a extra supersymmetric strong 
SU(2) gauge. This gauge theory is UV complete and calculable. Below a scale Ah, although 
the strong SU(2) gauge theory becomes non-perturbative, the theory can be described by the 
composite fields and their perturbative couplings. As a low energy effective Lagrangian, the 
following superpotential is dynamically generated, 

We{i = XN{H3H2-vi), (3.35) 

where is a SU(2)x, singlet composite held and Hi, H 2 are SU(2)i doublet composite helds. 
The coupling A is generated at X{Ah) ~ dvr, and it decreases towards the low energy scale due 
to the RGBs. Vq is given by vl ~ mAn/, where m is a fundamental parameter in the UV 
theory. The tadpole of N is generated and N obtains a its suitable vev. Therefore, Fat Higgs 
model also can solve /r problem, and the effective Lagrangian can be regarded as the nMSSM. 
One of the features of this model is that the coupling A is typically 0(1). 

The important thing is that the phenomenology of the nMSSM is almost independent of the 
symmetry of a UV theory (i.e. the desecrate symmetry, the U(l) Peccei-Quinn symmetry, the 
secluded U(l)' gauge symmetry, the strong SU(2) gauge symmetry, etc.). Therefore, we consider 
only the nMSSM Lagrangian Eqs. (3.27, 3.28) as a low-scale (TeV scale) effective theory in the 
next two chapters. We will discuss the phenomenology of the nMSSM, especially the dark matter 
and the baryon asymmetry of the universe. 
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3.4.2 Higgs Sector 

In the nMSSM, a Higgs potential is given as, 

Ho = m?|Hi|2 + m2|H2|2 + m||5|2 + A2|H2Hi|2 + A2|5|2(|Hi|2 + |H2|2) 

+ ^(|H2|2-|Hi|2)2 + ^|F|i?2p 

+{XAxSH 2 Hi + tsS + + H.c.), (3.36) 

where is defined as (f' = g'"^ + where g' {g) is the U(l)y (SU(2)) gauge coupling constant. 
In this potential, there are seven independent parameters, 

A, ml, ml, ms. Ax, mi2, ts- (3.37) 


Thanks to SU(2) rotation, we can take {H^ ) = 0 at elsewhere. Thus, when one takes = 0, 
the Higgs potential can be expanded as follows 


Vo 


= + ml{\H^\^ + + m||5p + X^\H^\‘^\H^\‘^ + 

+ + \H^\^ + \H+\^ - 2|H?|2|H^|2 _ 2\H^\‘^\H+\‘^ + 

8 

+ + {-\AxSH^H^ + tsS-ml^H^H^ + H.c.). (3.38) 


Next, the vacuum fluctuations of the scalar field are defined as 


1 


= vi + ^{HiR + iHu), 
^2 = V2 - 7={H2R 4-iH2l), 


S = s 


V2 

1 


V2 


{SR + iSi), 


(3.39) 


here we can choose (H^) = V 2 to be real and positive by phase redefinition of H 2 . The mini¬ 
mization conditions for electroweak symmetry breaking give the following conditions, 

m-i = ("ii2 + ^-4 as)*—- ^(|ui|^ - U 2 ) - |Ap(u2 + |sp), (3.40) 

vi 4 

"^2 = iml 2 +>^Axs)*—+ ^{\vi\‘^-vl)-\X\^{\vi\‘^+ \s\‘^), (3.41) 

V2 4 

ml = AM^^-^-|A|2(lui|2 + ui). (3.42) 

By the phase redefinition of H 2 H 1 , one can choose (mf 2 + AH^s) to be real and positive, and so 
it leads ui is real form Eq. (3.41). Similarly to the MSSM, let us dehne the parameter tan/3, 


tan/3 


Vl 


VEW 



+ vl 


174.1 GeV. 


(3.43) 
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In addition, the Eq. (3.42) leads to the vev of the singlet scalar, 

t*g- \*A\VlV2 


(3.44) 


Here, we can choose s and {t*g — \*A\viV 2 ) to be real by phase redefinition of S. At the this 
time, the Higgs potential Eq. (3.36) can obtain an e effective /r term and a effective H/r term 
(see Eq. (2.35)), 


Ateff = As, (3.45) 

(H/r)eff = 2(m^2 + AAAs). (3.46) 

Since m^2 — Cl( Afsusv^) and ts — C>(Msusy^)) Eq. (3.44) leads to s ~ O(Msusy)- Therefore, 
/Teff - C)(Msusy) and (H/r)efr ~ ©(Msusy^) are realized. 

Note that since this model has four complex parameter Eq. (3.37) and one can carry out 
three phase redefinition, one physical CP-violating phase remains in the Higgs sector. If one 
take Ax to be zero, there is no CP-violating phase in the Higgs sector. But relative phase 
with gaugino masses and other A terms, arg(m^2^5-A7gaugino) and a.Yg{m\ 2 t*s-^i)^ keep being the 
physical CP-violating phase [124]. 


Mass of Scalar bosons 

We assume that A, ts-, mf 2 and Ax are real in the following for simplicity. Then, a classifica¬ 
tion according to the CP-even or CP-odd for the scalar boson is justified. The singlet superfiled 
consists of singlet CP-even and CP-odd scalar components and singlet Majorana spinor compo¬ 
nent. Therefore, the nMSSM has 3x3 CP-even mass matrix, 3x3 CP-odd mass matrix and 
5x5 neutralino mass matrix, which are points of difference in the MSSM. 

Now, one can obtain the following mass matrix for the CP-even and CP-odd Higgs bosons 
from the Higgs potential Eq. (3.38), 


with 


-C 


CP—even/odd 


1 (Ml,, 

- {H,rH2rSr) MI,2 

' \mI„ 

1 (M\„ 

+- [ll,iR2iSi) 

' \Mj„ 


MI,2 Ml„\ ( H,r \ 
MI22 Ml,, 

Mhs MlJ \ Sr J 
M],2 M]„\ I H,J \ 
M]22 M]„ , 

-^?23 M]J \ Sj ) 


Ml,, = M\ sin^ /3 -|- M| cos^ /3, 

MI22 = M\cos^ fi + mIsIv? fi, 

MI„ = —— {ts — XAxVriy sin (3 cos P) = m'g + 


(3.47) 
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= ~2(Ml + Ml — 2\‘^vl^r)sm2f5, 

(3.51) 


= vew{2X‘^s cos f3 — XAxsin f3). 

(3.52) 


= VEW (2A^ s sin /3 — XA\ cos /3), 

(3.53) 

Mjii 

= Ml sin^ /3, 

(3.54) 

-Yf/22 

= Ml cos^ /3, 

(3.55) 

-Yf/33 

= — — (ts — XAxVew sin /3 cos /3) = ml X^Vew , 

(3.56) 

-Yf/12 

= ^MI sin 2/3, 

(3.57) 


= XAxVew sin l3. 

(3.58) 

-Yf/23 

= XAxVew cos l3. 

(3.59) 


where we have defined the following Ma, 

M\ = m\ + m2 + + X^Vew = ^ (3.60) 

olll 

The determinant of the CP-odd mass matrix is zero, and so the Nambu-Goldstone boson apperas. 
At the large Msusy case, these mass matrices give 

mn ~ Ml, mjjs ~ m|, 

m\ ~ M| cos^ 2/3 -I- ^ 

V ^s/ 

ml ~ Ml, mis ~ ( 3 - 61 ) 


where HS (A5) is the heavy CP-even (-odd) scalar boson, whose component is mainly S, and 
h is the SM-like Higgs boson. Note that the mass of the SM-like Higgs boson can become large 
at a low tan/3 region in comparison with the MSSM. It is given by an additional F-terms as we 
have discussed in Section 2.3.1. 

However, the radiative contributions to the Higgs boson mass are still important. In order 
to estimate the mass of the SM-like Higgs boson including two-loop radiative corrections, we 
have extended the two-loop level calculation using the RGE of the MSSM, which is given in 
Section 2.3.1.1. In the nMSSM, since the SM-like Higgs boson receives an extra F-term con¬ 
tribution to the Higgs quartic coupling Aquartic, there is a sizable tree-level contribution to the 
Higgs boson mass. When integrating out heavy SUSY particles, the matching condition at high 
scale from the SUSY to the SM is shifted by [57] 

Aquartic( Y/suSy) = Alo( MsvSy) + (3Aquartic( MsvSy) + (4y,-)2 ^NLo( APsUSy), (3.62) 


with 


(3A 


quartic 




A^ - Al 
2 


sin^ 2/3, 


(3.63) 
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where Alo is given by Eq. (2.70), and Anlo is given by Eq. (2.72). Therefore, at large A and 
small tan /3 can give an additional sizable contribution to the Higgs boson mass. Note that this 
extra contribution can be controlled by Ax. We have not considered extra loop corrections in 
Anlo term, which depend on the coupling A. When A is larger than 1, these extra radiative 
corrections become signihcantly contributions to the SM-like Higgs boson mass [125]. In this 
thesis, we have taken A < 0.8 due to consideration to a Landau pole constraint below the GUT 
scale^^, and thus these extra radiative corrections do not give the signihcantly contributions. 


Next, one can obtain the following mass matrix for the charged Higgs boson from the Higgs 
potential Eq. (3.38), 


-C 


charged 


where 



+ (Ma- - 






(3.64) 


= — sin/3i7]^ * — cos 

iVG+ = -cos/377f*+sin/3772^. (3.65) 

Thus, the mass of the charged Higgs boson is given as, 

= M\ + M'^ — (3.66) 

Note that the mass of the charged Higgs boson is lighter than the one of the MSSM by 
term (see Eq. (2.47)). Thus, the theoretical condition M\ + is needed to avoid a 

vev of the charged Higgs scalar held in the nMSSM. 


3.4.3 Landau Pole Constraint 

A large A coupling can raise the mass of the Higgs boson at the tree level. However, when 
the A coupling is too large at low energy, it causes a Landau pole (or Landau singularity) 
at a higher energy scale. The Landau pole means that some dimensionless running coupling 
constants become non-perturbative couplings at the hnite energy scale, and so a perturbative 
approximation is broken down by the strong (non-perturbative) couplings. In the nMSSM, a 
one-loop RGE for A is given as 

IMQ ^ (47r)2'^ + yl- 9'^ - • (3-67) 

This RGE implies that an absolute value of A monotonically increases with energy scale Q, and 
so A may eventually become non-perturbative. Let us demand that the theory is perturbative (A 

’^^The detail of the Landau pole is written in the next subsection. 
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tanfi 

Figure 3.1: The upper bound on the coupling A (Amax) under the condition of the no Landau 
pole up to the GUT scale using two-loop RGBs. The horizontal axis is tan f3. Amax is the values 
at Q = Msusy- We take Msusy = ITeV (dotted line), lOTeV (dashed line) and 100 TeV (solid 
line). 


and other dimensionless couplings do not blow up) below the GUT scale, Mqut — 2 x 10^® GeV. 
Then, one can obtain the upper bound on A at the soft SUSY breaking scale. 

Figure 3.1 shows the upper bound on the coupling A (Amax) under the condition of the no Lan¬ 
dau pole up to the GUT scale as a function of tan /3. Here, we have demanded gf {Mquj')/A ir < 1, 
where gt denotes all dimensionless couplings. Here we have used two-loop RGBs for calculations 
of the running couplings. The two-loop level RGBs for all couplings are summarized in Appendix 
B.1.2. Amax is the values at Q = Msusy, where we have used the SM RGBs below the soft SUSY 
breaking scale Msusy and the matching condition at electroweak scale (see Section 2.3.1.1). We 
take Msusy = ITeV (dotted line), 10TeV (dashed line) and 100TeV (solid line). 

We find that the A coupling should be smaller than 0.75-0.85, which depends on the soft 
SUSY breaking scale. Note that this bound can be alleviated by introducing an additional 
gauge symmetry or extra particles [126]. It is because that the contributions of additional gauge 
coupling to the RGB for A is negative or the ordinary gauge couplings g ', g become large at the 
high energy scale due to the extra matters, and so its contributions to the RGB for A is negative 
(see Bq. (3.67)). 

At low tan/3 region, tan/3 < 1.5, the upper bound on A decreases drastically. It is because 
that the top Yukawa coupling becomes large values at the low energy scale, yt = rnt/{vEW sin/3). 
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then it causes the Landau pole of top Yukawa coupling at the high energy scale. While, at large 
tan/3 region, tan/3 > 65, the upper bound on A also decreases drastically. It is caused by 
the bottom/tau Yukawa couplings, which become large values at the low energy scale, = 
m^irlivEW cos/3). 
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Chapter 


4 

Resonant Singlino Dark Matter 


This chapter is based on the work by the author [17]. We consider a singlino 
dark matter scenario in the nearly minimal supersymmetric standard model. 

We find that with high-scale supersymmetry breaking the singlino can obtain 
a sizable radiative correction to the mass, which opens a window for the dark 
matter scenario with resonant annihilation via the exchange of the Higgs boson. 

We show that the current dark matter relic abundance and the Higgs boson 
mass can be explained simultaneously. This scenario can be probed by the 
search of the Higgs invisible decay and the direct direction of the dark matter. 

We have shown that the nMSSM can solve the // problem, the domain wall problem and the 
tadpole problem simultaneously in the previous chapter. Next two chapter, we will consider the 
phenomenology of the nMSSM, which are the based on the works by the author [17,18]. As 
a results, both the two studies conclude that the nMSSM with a high-scale SUSY breaking is 
valid. 


4.1 Dark matter in the nMSSM 

In this chapter, we focus on the dark matter phenomenology in the nMSSM. This chapter is 
based on the work by the author [17]. First, we briefly review a situation of the dark matter in 
the nMSSM, and we also explain why we have considered it. 

Recent various cosmological observations have established the ACDM cosmological model 
and the relic abundance of the cold dark matter is measured accurately by WMAP and Planck [127, 
128]. In the nMSSM, as discuss later, the singlino, which is the fermionic component of the extra 
gauge singlet superfield, can be a candidate of the dark matter [15,129-133]. In fact, the singlino 
mass and its couplings with SM particles are suppressed by the soft SUSY breaking scale in the 
nMSSM, and such a dark matter leads to the overabundance in the universe. Therefore, the 
singlino dark matter scenario seems to be incompatible with relatively high-scale (TeV scale) 
supersymmetry breaking, which is inferred from the measured SM Higgs boson mass [134,135] 
and the null signals of the sparticle searches at the LHC [136,137]. 

However, if one-loop corrections to the singlino mass are taken into account, this situation 
will change. As will be shown later, similar to the mass of the SM Higgs boson in the SUSY 
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models, the singlino can obtain a sizable mass, which opens a window for a resonant dark matter 
scenario via the s-channel annihilation with the exchange of the SM Higgs boson. 

Let us consider the helicity of the annihilation process of the neutralino dark matter to 
fermions XX //■ Since neutralinos are Majorana fermions, the helicity of the two neutralinos 
are opposite to each other in the center-of-momentum frame when the case that neutralino pair 
annihilates via the CP-even scalar boson s-channel exchange. While the helicity of the final state 
two fermions are facing in the same direction. Therefore this process needs either the p-wave 
(L = 1) suppression or chirality suppression. Since the annihilation rate of the singlino via the 
SM Higgs boson s-channel exchange is also p-wave snppressed, one needs a relatively large value 
of SM Higgs boson-singlino coupling compared with the scalar dark matter in the Higgs portal 
model [138]. This fact implies that the singlino dark matter could be probed more easily than 
the scalar one. 


Now, let us focus on the phenomenology of the singlino in the nMSSM. The singlino partic¬ 
ipates in a member of usual neutralino as a new gauge eigenstate. At the tree level, the 5x5 
neutralino mass matrix in the basis = (H, H^, H^, S) is given by 


jr = -^X^M£,,X + R.c., 
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Ml 
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g ui 

V2 

g U 2 

V2 

0 \ 
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gri 

^/2 
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g'ri 
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V2 
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ReE 

-XV2 


g V 2 

V2 

gv2 

V2 

/^eff 

0 

—Xvi 

V 

0 

0 

-XV2 

—Xvi 

0 / 


(4.1) 


(4.2) 


where S is the fermionic component of S. Since a determinant of this mass matrix is nonzero, all 
(five) neutralinos obtain a nonzero mass. If one impose additional matter parity Pm, the lightest 
neutralino becomes stable and can be a candidate for the dark matter in the universe^^. The 
one of the origin of the matter parity Pm is a remnant discrete subgroup of the local U(1 )b-l 
symmetry. This U(l)s_i symmetry is broken above the electroweak scale. The neutralino mass 
matrix can be diagonalized by a unitary matrix N, 


= diag(m2 q , 0 , 0 , o, m Jo), 


(4.3) 


here we call Xi the lightest neutralino. 

The desecrate R-symmetry does not contain the R parity Eq. (2.21). In fact, although the R parity is 
conserved accidentally in the renormalizable teams, it is broken in the non-renormalizable terms: SSLH 2 , SLLE 
and SLQD, which can not lead to an observable proton decay. Considering the effects of such as R parity breaking 
non-renormalizable (Planck suppressed) operators, the life time of the lightest neutralino is longer than the age 
of the universe [14]. However, there are experimental bound on the lifetime of the dark matter from the cosmic 
ray searches, which is longer enough than the age of the universe. It implies that the life time of the dark matter 
which is comparable to the age of universe have been excluded. Therefore, we should impose the additional matter 
parity or some symmetries which do not forbid the tadpole. 
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From the tree-level calculations, in a very well approximation the mass of the lightest neu- 
tralino is given as [139], 


= Min 

' Xi ' 


B - - 4C B + \/52 - AC 


(4.4) 


where 


Ml Ms / z/2 

Ml + Ms ^ 


M| M \ 

+ 1/2 Ml + Ms J 


/ieff sin 2/3 


(Ml -|- Ms)(/igg -|- 


(4.5) 


C 




Ml Ms 
Ml + Ms 


/ieff sin 2/3 


M 

Ml +Ms 



(4.6) 


with and M = Micos^Ow + Mssin^^iy. When the case M^ <C //eff and 

Mz Mgaugino) a mass-eigenstate neutralino whose component is mainly S becomes the lightest 
neutralino. We denote s as the mass-eigenstate neutralino whose component is mainly S. Note 
that, here and in the following, we call s as a “singlino” in order to make understanding easy. 
Then, the mass of the singlino is evaluated by expansions of Mz/HeS and M^/Mgaugino, 




^tree ^ 


^EW 


Msusy 


sin 2/3, 


(4.7) 

(4.8) 


where we denotes the typical soft SUSY breaking scale by Msusy and we use the fact that a value 
of /igff become 0( Msusy) Eq. (3.45). As you can see, the mass of the singlino has suppressions 
by soft SUSY breaking scale and by sin 2/3. Therefore, when the soft SUSY breaking scale is 
relatively high (Msusy ^ ITeV) as suggested by the LHC experiments [134-137], the singlino 
becomes the LSP and it can be a candidate of the dark matter. Furthermore, by the tree-level 
analysis one can see < 50GeV. 

Since the singlet superfield S interacts only the Higgs multiplets, the singlino can be coupled 
with the SM particles only through the mixing with Higgsinos. Thus, the singlino-SM particle 
coupling are suppressed by the soft SUSY breaking (Higgsino mass) scale in the nMSSM. More¬ 
over, since the singlino is the LSP it can not decay. Generally, an annihilation cross section of 
such a stable particle is small and it would freeze-out at relatively early time in the universe. 
In other words, the relic density of such a stable particle would be overabundant. Therefore, 
the singlino dark matter typically leads to the overabundance in the universe, and one needs to 
dilute the relic density of the singlino by some mechanisms. 

In the literature, it is known that there are two solutions for avoiding overabundance of the 
singlino dark matter. First, when the lightest GP-odd Higgs boson Ai is dominantly singlet¬ 
like and its mass is mAi ~ 2ms, tbe singlino annihilation cross section can be resonant via 
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the s-channel Ai exchange [132,140-143]. This resonant annihilation cross section gives much 
dilution of singlino dark matter, then the relic density is suppressed. After discovering the SM 
Higgs boson /i, however, this scenario is severely constrainted from branching ratio of h —)■ ss 
(invisible) and h —>• AiAi —)■ All [142]. Strictly speaking, if is light, the mass spectrum 
becomes rris < w,/i/2 < rriAi ~ 2ms, then the SM Higgs boson decay to the two singlino is 
kinematically allowed. While, if //eff is heavy, the singlino mass becomes suppressed and the 
mass spectrum becomes < mAi ~ 2ms < rnh/2, then the SM Higgs boson decay to the two 
Ai is also kinematically allowed. Such an additional decay channel of the SM Higgs boson is 
currently constrained from the measurements of the Higgs coupling. 

Next, when mz ~ 2ms, tbe singlino annihilation cross section can be resonant via the s- 
channel Z boson exchange [131,144]. This scenario is constrained from Z ^ ss and h —)• ss. 
Furthermore, this scenario leads to light /ieff in order to obtain the sizable singlino mass mg ~ 
45GeV Eq. (4.8). Although such a light soft SUSY breaking scale is favored with naturalness, 
it is disfavored from the point of view of the SUSY flavor/CP problem. 

However, these arguments are incomplete. It is because they are given by the tree-level 
analyses, and in fact, we find at first time that the one loop radiative corrections give significantly 
contributions to the mass of the singlino. We will show the numerical analysis of the full one-loop 
level singlino mass in Section 4.5. Interestingly, we find that this radiative correction is roughly 
proportional to the soft SUSY breaking scale, oc MsusY/(4vr)^. Thus, this contrition 

can dominate the singlino mass in relatively large MgusY- Furthermore, thanks to the radiative 
corrections the singlino mass can reach 62 GeV, which is half of the mass of the SM Higgs boson. 
Therefore, the singlino annihilation cross section can be resonant via the s-channel SM Higgs 
boson exchange. This is a new scenario in the nMSSM. 

In next section, using effective Lagrangian, we will calculate the thermal relic abundance 
of the singlino dark matter and experimental constraints from Higgs invisible decay searches 
and from direct dark matter searches. In order to compare with literature, we will consider the 
resonant case with the Higgs boson exchange and with the Z boson exchange. In Section 4.5, 
we will evaluate the full one-loop singlino mass. 


4.2 Resonant singlino dark matter via SM Higgs boson 

In this section, using the low energy effective Lagrangian we calculate thermal relic abundance 
of the dark matter which annihilate via the SM Higgs boson or the Z boson. 

Let us consider the case where only the singlino s is light and other SUSY particles are 
relatively heavy, vew ^ AlsusY- Note that the masses of the singlet scalars are also heavy 
in the nMSSM with vew ^ AIsusY) see Eq. (3.61). In this case, the low energy effective 
Lagrangian includes singlino s, the SM Higgs boson h, Z boson, other fermions and other gauge 
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Figure 4.1: One of the Feynman diagrams for the resonant annihilation of the singlino dark 
matter via the SM Higgs boson exchange. 


bosons, and it can be written as 

-£es ^ ^ss + ^hss + , (4.9) 

where singlino s is written by Dirac 4-component spinor, Aeg is the SM Higgs-singlino effective 
coupling and gzss is the Z boson-singlino axial effective coupling^^. In this Lagrangian, the 
singlino dark matter can annihilate to SM particles, via the s-channel exchange of the SM Higgs 
boson or the Z boson. In the following, we estimate the thermal relic abundance of singlino dark 
matter with this effective model regarding Aefr, 9zx^x° parameters by solving 

Bolzmann equation [145]. 


The resonant case with the Higgs boson exchange 

First, we consider a dark matter annihilation via the SM Higgs boson exchange, ss —)• /i* —>• 
All. Figure 4.1 represents a relevant diagram for the resonant annihilation of the singlino dark 
matter. In order to obtain an annihilation cross section of the singlino, we use the optical 
theorem. 


a{ss —)• /i* —)■ All) 


ImA4(ss —^ /i* —^ ss) 
“^EcmPcm 


(4.13) 


Let us derive this neutralino-neutralino-Z boson axial vector coupling. Generally speaking, the neutralino- 
neutralino-Z boson coupling is given as 


C 


-.Ok+Ot 


_ QE _ O * 

0 ^IJ U ryU. I 

Xi - 1 - X- 


-\L,* 


u r7U 


(4.10) 


with 


O^. =_^_ 

2 cos 6w 


{-NisNjs + NiiN;R. 


(4.11) 


When the neutralino mass matrix includes no CP-violating phase, one gets Ofj = Ofj*. Therefore, we obtain 

.9zx^x9 


P = -x°- 


-7m75Xj 


0 T’M 

A Zj . 


(4.12) 


with gzx°x'^ = 
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with Ecm = \/sj PCM = Y 'S(l “ 4m?/s)/2. We assume the following full propagator of Higgs 
boson, 


2 I • ! —rn / \ ’ 

where rrih = 125.5 GeV and 


(4.14) 


rh(s) = r;,,o—, (4.15) 

ruh 

with F/j o = 4.07 X 10“^ GeV [146]^^. Then the pair annihilation cross section of the singlino 
can be obtained as follows, 


a{ss —>• /i* —)■ All) 


2 V s {s — + sr| 


(4.17) 


Thus, one can see a{ss h* ^ All) oc Vg. Since the freeze-out of the WIMP dark matter occurs 
when they are non relativistic (roughly estimation, v ~ 0.3), this dependence of the velocity of 
the singlino gives suppression to the annihilation cross section. In other words, this annihilation 
rate is p-wave suppressed. To distinguish whether the dark matter annihilation cross section is 
s-wave process (not suppresed) from the interaction structure of the dark matter is possible, and 
it is summarized in Ref [147]. For example, xxff is p-wave annihilation, xlbXflbf (GP-odd 
Higgs exchange) and </><(>// are s-wave annihilation, and so on. 

The thermal average of the annihilation cross section times the relative velocity of the anni¬ 
hilating particles can be obtained as follows. 


^ (T{s)vre\ 

> j (Pp^d?p^e-Ei/Te-E2lT 


(4.19) 


One can reduce this formula to the single-integration [145], 


{aVre\){T) 


1 

Sm^T{K2[m/T]Y 



AnY)Ki 



(4.20) 


where Ki[x\ are the modified Bessel functions of order i [148]. 


Although Eq. (4.15) is our assumption, this is satisfied in leading order. The dominant decay channel of the 
SM Higgs is bb (Br(h —>■ bb) ~ 0.5). This transition rate is 


r(s) 



(4.16) 


Therefore, the decay width is roughly proportional to Ecm around rrih = 125.5 GeV. 

^'^Strictly speaking, the relative velocity in the thermal average (...) is not the nonrelativistic relative velocity 
Vr = \vi — V 2 \ but the so-called Mpller velocity, 

V = Y (iTi — V2Y — (I’l X 1 ) 2 )^. (4.18) 
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Using thermal average of the annihilation cross section, the relic density of the singlino in 
the expanding universe can be evaluated by solving the following Bolzmann equation, 

^ + SHng = -{aVrei){T){nl - n|eq), (4.21) 

where H is the Hubble parameter, and ng is the number density of the singlino. The details of 
the calculation for the relic density of the dark matter are written in Ref. [149]. The Hubble 
parameter is related to the total energy density during the radiation-dominated era. 


i da ISttGp 

a dt V 3 


(4.22) 


where a is the cosmological scale factor, G is the gravitational constant G = 1/Mpj = 6.708 x 
10“^®GeV“^, and the total energy density p = (Tr^/SO)^* (T)r^. (7*(T) denotes the relativistic 
degrees of freedom {m <T). We have used the fitting formula of g*{T) of Ref. [150], where 
in the reference corresponds to 

In Figure 4.2, the black lines show the ratio of the thermal relic abundance Vigh? to the current 
dark matter density = 0.1199 [128] where we take the Higgs boson mass as mh = 125.5 
GeV. The horizontal axis is the mass of the singlino mg and the vertical axis is the singlino- 
singlino-Higgs coupling Aeff. Note that we have handled mg and Aeg as free parameters. The 
thick black line represents the appropriate parameter region in which the dark matter relic 
abundance is consistent with the current dark matter abundance. While, the singlino relic 
density overdose the universe at the dark-shaded region. As discussed previous section, since 
the couplings of the singlino with SM particles are too small and the annihilation cross section 
is also too small, the singlino dark matter could not be sufficiently diluted in the universe. On 
the other hand, around mg ~ mhj2, thanks to the s-channel exchange of the SM Higgs boson 
the annihilation cross section becomes large, and the singlino dark matter can be appropriately 
diluted (white region). 


This effective potential Eq. (4.9) receives two kinds of experimental constraints: constraint 
from the branching ratio of the Higgs to singlino pair (Higgs invisible decay searches) and the 
direct detection of the dark matter. The decay width of the Higgs to singlino pair is given by 


A^ 

r(/i —)• ss) = —^mh 

^ ^ 167r 



Thus the branching ratio of the Higgs to singlino pair is 


(4.23) 


Br(/i —)■ ss) 


T{h ss) 
r/i,o + r(h —)• ss) 


(4.24) 


’^^Although the value of is one of the SM, one can use it in the nMSSM. In our scenario the mass of the 
singlino dark matter is about 60 GeV, and so its freeze out temperature is 0(1) GeV. Therefore, there are no 
sparticle contributions to gr*. 
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Figure 4.2: The singlino thermal relic abundance and experimental constraints/future prospects 
in the case of the singlino resonant annihilation via the SM Higgs boson s-channel exchange. 
The black contour denotes the ratio of the thermal relic abundance Ctsh? to the current dark 
matter density Ctch? = 0.1199 [128]. The singlino thermal relic density overdose the universe 
at the dark-shaded region. The regions above the red solid lines are excluded by the Higgs 
invisible decay {h —)• ss) searches of CMS (Hr)]''"- < 58 %) [151] for upper line (yellow-shaded 
region) and by the global fit of the Higgs couplings (19 %) [152] for lower line. The dashed 
red lines correspond to the future sensitivity of high luminosity LHC (6.2 %) [153] and ILC 
with C = 1150fb“^ at y/s = 250GeV (0.4 %) [154]. The blue-shaded regions are excluded by 
XENONIOO [155] and LUX [156]. The regions above the blue and the green dashed lines can 
be probed by the future direct dark matter searches of LUX [157] and XENONIT [158]. As the 
future sensitivities of XENONIT, we have used 0.14 (top), 0.326 (middle) and 0.66 (bottom) to 
the parameter /^r. 



In collider, this search can be particularly performed on the Higgs boson produced via association 
with Z boson. 


pp{e~^e ) —)• Z* —Zh ii + invisible. (4.25) 

In the Figure 4.2, the regions above the red solid lines are excluded by the Higgs invisible 
decaysearches of CMS, Br(/i —)• ss) < 0.58 at 95 % C.L. (upper line) [151], and by the global 
fit of the Higgs couplings, Br(h —^ ss) < 0.19 at 95 % C.L. (lower line) [152]. The regions 
above the red dashed lines can be probed by the future Higgs invisible decay searches of high 
luminosity LHC, Br(/i —)• ss) < 0.062 at 95 % (upper line) [153] and ILC at y/s = 250 GeV, 
1 at“^, Br(/i — ss) < 0.004 at 95 % (lower line) [154]. 
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4.2 Resonant singlino dark matter via SM Higgs boson 


The direct dark matter searches can set limits on the spin-independent cross section of dark 
matter-nucleon elastic scattering, sN —>■ sN. If one integrate out the SM Higgs boson, effective 
Lagrangian becomes 





CXg 

4ttvew 


G 


fll2 



(4.26) 


where fi are the SM fermions, and is a field strength of SU(3) gauge. Note that, here 
we have not considered the Z boson exchange in the spin-independent cross section of dark 
matter-nucleon elastic scattering. It is because the Z boson-singlino coupling gzsS always has 
more suppression by Msusy than the Higgs-singlino coupling Aefr (see Eqs. (4.40), (5.72)). 
We have also not considered the squark exchange since we have assumed sparticles are heavy 
enough. Using this effective Lagrangian, we can obtain the cross section of dark matter-nucleon 
spin-independent elastic scattering as follows [138,159], 




^ — ''ett f2 '^ n'^s 

+ mNY' 


(4.27) 


with = 0.939 GeV and 


= {N\Y,mgqq - ^G^,G>^''\N). (4.28) 

<? 

We use /at = 0.326, which is the lattice result [160]. Then, 

2 

aCsN^sN) = 8.051 X 10-^3 ^ [cm^] (4.29) 

~ 0.8 X 10"^® (^) 

Note that, there are various estimation of the parameter in the literature. It is because a 
treatment of the heavy quarks, especially the contribution of the strange quark to the nucleon, 
is difficult. These results in the literature give a range to the value of fjy as follows [161], 


0.14 < fN < 0.66. 


(4.31) 


Including this uncertainty, the cross section of dark matter-nucleon spin-independent elastic 
scattering is gives as, 

a{sN sN) ~ (0.15 - 3.3) x 10"^® (^) 

In the Figure 4.2, the blue-shaded regions are excluded by the current direct dark matter 
searches of XENONIOO [155] and LUX [156]. The region above the blue dashed line can be 
probed by the future direct dark matter search of LUX [157]. For applying these constraints and 
future prospects, we assume = U.ch?'. Here we have used /at = 0.326. The regions above the 
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green dashed lines can be probed by the future direct dark matter searches of XENONIT [158]. 
As the future sensitivities of XENONIT, we have used 0.14 (top), 0.326 (middle) and 0.66 
(bottom) to the parameter /w. One can see that the region where s is consistent with the 
current dark matter relic abundance lies around Aefr ~ 0(0.01) and mg ~ 60 GeV. In this region, 
resonant pair-annihilation of s occurs via the Higgs boson with mg ~ mhj2. This allowed region 
can be covered by the future Higgs invisible decay searches and direct dark matter searches, 
especially by XENONIT. 


The resonant case with the Z boson exchange 

Next, let us consider a dark matter annihilation via the Z boson exchange, ss ^ Z* —)• //. 
Then the pair annihilation cross section can be obtained as follows [162], 

a{ss —>• Z* —)■ //) = (4.33) 

^s(s-4m|) 

The Lorentz-invariant function OJ^ (s) is defined as 

w^(s) = i^cf0{s - 4mj)l3f{s, mf)u}f (s)) , (4.34) 

/ 

where Nc^ is the color factor (Nc^. = 3 for quarks and = 1 for leptons), and a kinematic 
factor /?/ is given as 


with 


/3/(s,m/) 



(4.35) 


wf(s) = 


9Zs 


12|C; 


Mi, 


2\2 


(s - Mz) 


s - M| -I- iVzMz 

+ (\C(/^\‘^{s + 2m‘j) + - 4mj)^ (s - 4m|) 


(4.36) 


Here total with of Z boson is 2.4952 ± 0.0023 GeV, and Z boson-fermion couplings are 


C 


iffZ _ 

9 


'V 

2 cos 6w 

iffZ _ 

9 

^3,/- 

A 

2 cos 9w 


{Tsj - 2sin^ OwQf), 


(4.37) 

(4.38) 


Substituting the pair annihilation cross section via the Z boson exchange Eq. (4.33) into the 
thermal average of the annihilation cross section Eq. (4.20), and solving the Bolzmann equation 
Eq. (4.21), we have estimated the thermal relic abundance of the singlino dark matter. 
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Figure 4.3: The singlino dark matter thermal relic abundance in the case of the singlino 
resonant annihilation via the Z boson s-channel exchange. The black contour denotes the ratio 
of the thermal relic abundance to the current dark matter density = 0.1199 [128]. 
The singlino thermal relic density overdose the universe at the dark-shaded region. 


In Figure 4.3, we show the singlino dark matter thermal relic abundance in the case of the 
singlino resonant annihilation via the Z boson s-channel exchange^®. The black contour denotes 
the ratio of the thermal relic abundance to the current dark matter density = 

0.1199 [128]. Thus, the singlino thermal relic density overdose the universe at the dark-shaded 
region. Similar to the case of the resonance annihilation via the Higgs boson s-channel exchange, 
the greater part of the parameter region except mg ~ 44^/2 are suffering from the overabundance 
of the universe. We find that gzss can become small to O(10“^-10“^) at the case when one 
assume the singlino annihilation cross section is resonant by the Z boson exchange. While, one 
can estimate gzss at 


gzss = 




g 


2 cos 9-\/y 

Vg 


2 cos 9-\/y 


+ iV54A^54) 


'^EW 

Msusy^ 


COS 2/3. 


(4.39) 

(4.40) 


This estimation implies Msusy ^ 1-2 TeV at the resonant case with the Z boson exchange. 

However, as discuss next section, the Higgs-singlino coupling can also be estimated relating 
to Msusy- It implies that when Msusy ^ 1-2TeV, Agff > 0.1 (see Eq. (4.42)). Such a large 

’^®Similar to the resonant case with the Higgs boson, although there is the upper bound on the invisible Z width, 
r(Z — >■ XiXi) < 1.76 X 10~^ GeV [163] as a constraint from the experiment, we have not calculated this bound. 
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Agff and light singlino mass {rris ~ 45GeV) may be excluded by current experimental bounds^^ 
(see Figure 4.2). Therefore, in terms of the overabundance of the singlino dark matter in the 
universe, the resonant scenario with the SM Higgs boson is the last resort for the nMSSM. 


4.3 Radiative Singlino mass 


Now, we calculate the singlino mass rris and the Higgs-singlino coupling Agg in the nMSSM. 
From the tree-level calculations, these values are evaluated as 


+g' - g {Z^iN^^Nl, + 

~ V2A^sin2/3, 

Msusy 

where is a unitary matrix which can diagonalize the CP-even mass matrix as follows 


(4.41) 

(4.42) 
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-^Il2 

Ml,, 

-^Il2 

-^l?22 

Ml,, 


■^It23 

Ml,, 




diag ■ 


(4.43) 


where rrii is the lightest CP-even Higgs mass, that is the SM-like Higgs boson. On the other 
hand, the tree-level singlino mass is evaluated as Eq. (4.8). Obviously Agg ~ 0(0.01) and 
nis ~ 60 GeV can not be satished at the same time. 

However, one-loop corrections to the neutralino mass can raise the singlino mass with rel¬ 
atively large Msusy- The typical diagram, Higgs-Higgsino loop diagram, which contributes to 
the singlino mass is given in Figure 4.4. Note that a sum of the contributions of the neutral 
Higgs-Higgsino vanishes at the leading order. The vertex of the CP-even Higgs-Higgsino-singlino 
is different from the one of the CP-odd Higgs-Higgsino-singlino in only “i”. Thus, since the loop 
contribution is proportional to the square of the vertex, the CP-even Higgs loop contribution 
cancel out the CP-odd Higgs one up to their mass dependence. So, one can easily estimate the 
radiative corrections to the singlino mass by calculating the charged Higgs-Higgsino loop. This 
loop gives the following contribution to the singlino mass, 


m 


1-loop _ 


(47r)2 

A^ 

{AttY 


/ieffsin2/3 • F 

Msusy sin 2/3, 


2(mf2 -F4A/reff) \ 
sin 2/3 J 


where the loop function F (x) is dehned as 


(4.44) 


F(x) 


xlogx 

-I ? 

X — 1 


(4.45) 


’^^Strictly speaking, they are not excluded at the case when tan /3 is large and there is a extra contribution to 
the singlino mass which is beyond the nMSSM. It is because the Higgs-singlino coupling Aeff is suppressed at large 
tan/3 (see Eq. (4.42)), but the Z boson-singlino coupling gzss is not suppressed. 
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Figure 4.4: 


Typical one-loop diagram which contributes to the mass of the singlino. 


and it satisfies F (1) = 1. We hnd that radiative corrections is proportional to the soft SUSY 
breaking scale. It is because this loop diagram should have a chirality flip on the Higgsino 
propagator, and this flip gives Msusy to a numerator of the loop contribution. 

Next, we numerically calculate the neutralino 5x5 mass matrix including full one-loop 
corrections [164, 165]"^®. We assume that the all complex parameters to be real in neutralino 
mass matrix Eq. (4.2), then the neutralino mass term becomes 

^ = -^^X-^UeeV’x> (4-46) 

where V'x is the Dirac 4-component spinor = ((|J, (|j). At one-loop 

level, the radiative corrections to the neutralino sector are given by [164,166] 


A< = - 


i’x (i^ 




V’x + 


(47r) 


(V’xP^S^Pi'i/'x + V’xP^S'^^kV’x + + H-c.) 


(4.47) 


where the correction comes from a self-energy for neutralinos, and the corrections 

SVH(p2) 

come from the wave function constants for the neutralinos. The momentum of the 
external line is represented by p. A pole mass of the neutralino at one-loop level is obtain by 
the following equation. 




2 (47r) 


S^(p2) + (S^(p2))^ + ((S^(p2))^ + S^(p2)) 


+ MCe (u^(p^) + (SV))^ 


(4.48) 


The explicit formulae of the one-loop corrections to the neutralino S‘®^(p^) and S^/^(p^) are given 
in appendix B.2. Here, we have performed DR renormalization in which we have subtracted the 
In the limit of k = 0, one-loop corrections in the NMSSM reduce to the one in the nMSSM. 
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MsusyfTeVl 


Figure 4.5: The singlino mass, which is the lightest mass eigenvalues of the neutralino mass 
matrix, at the tree level (dashed) and the full one-loop level (solid) as a function of Msusy- We 
take A = 0.75, all dimensional parameters equal to Msusy and tan/3 = 2 (red), 5 (blue) and 10 
(green). 


1/e poles (see Eqs. (B.91, B.92)). Note that 55 component of the neutralino mass matrix does 
not diverge, which is consistent with the fact that the quadratic term of the singlet superfield 
is not included in the superpotential. The cancelation of the divergence of the neutral Higgs- 
Higgsino loop is trivial, since the vertices are different in only “i” as discussed before. While, 
the cancelation of the divergence of the charged Higgs-Higgsino loop is rather non-trivial. The 
charged Higgs-Higgsino-singlino coupling is proportional to sin /3 cos /3, on the other hand the 
charged NG boson-Higgsino-singlino coupling is proportional to — sin/3cos/3 (see Eq. (3.65)). 
Since the divergence is independent of the mass of the virtual particle, thus, all divergences 
are canceled out. This cancellation can be understood in terms of the supergraph. One-loop 
correction to the singlet superfiled is roughly given as 


T ~ 


de^ 


1 

(47r)2 Msusy^ 


( 5 ^ 5 ) 2 , 


(4.49) 


and thus the 55 component of the neutralino mass matrix is hnite. 

Figure 4.5 shows the dependence of the singlino mass to MgusY; which is the lightest mass 
eigenvalues of the neutralino mass matrix, at the tree level (dashed) and the full one-loop level 
(solid). In this figure, we take A = 0.75, all dimensional parameters equal to MgusY and 
tan/3 = 2 (red), 5 (blue) and 10 (green). We find that the singlino obtains sizable one-loop 
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corrections to the mass when Msusy ^ 1 TeV. Since this feature is due to the suppression of 
the singlino mass at the tree level, the two-loop level corrections to the singlino mass is estimated 
to be smaller than the one-loop one. We also find that both the tree level and the one-loop level 
mass have a tan /3 suppression. 

We have checked the validity of the our full one-loop calculation code by three ways. First, 
we have compared the full one-loop result with the estimation by Eq. (4.44), then we find that 
these results agree with very well. Second, we have checked that 55 component of the full one- 
loop correction matrix vanishs in the SUSY limit, Msusy 0. Third, we have chosen the Ax 
to vanish the combination m ^2 + varying the soft SUSY breaking scale. Then we find that 
even if Msusy is much higher scale, 55 component of the full one-loop correction matrix vanishs. 

Note that with Msusy ~ 0(10) TeV, tan/3 ~ 0(1) and A ~ 0(1), one can simply obtain 
Aoff ~ 0(0.01) and ~ 60 GeV^®. Moreover, the Higgs boson mass becomes around 125 GeV 
in such parameter sets with the help of the additional quartic coupling A. We will show these 
validity by using the numerical calculations in the next section. 

4.4 Numerical Results 

In this section, we numerically investigate the singlino resonant dark matter scenario and the 
Higgs boson mass in the nMSSM. We calculate the mass of the Higgs boson using the two-loop 
renormalization group equation including the matching condition Eq. (3.63) [57]. The detail of 
the calculation for the Higgs boson mass is given in Section 2.3.1.1. 

In Figure 4.6, we show the singlino mass mg (red lines), the effective Higgs-singlino dark 
matter coupling Aefr (blue lines) and the Higgs boson mass mh (black dashed lines) in Msusy- 
tan /3 plane. Here, the singlino mass mg is obtained as the lightest mass eigenvalue of the 
one-loop full neutralino mass matrix Eq. (4.48), and Aeff is given by the tree-level calculation 
Eq. (4.41). For simplicity, all parameters are chosen to be real. The trilinear coupling A is 
taken to be Amax which is a maximal value avoiding Landau singularities up to the GUT scale, 
2 X 10^® GeV (see Section 3.4.3). All soft SUSY breaking parameters except Ax are set to 
Msusy (Ac^ = eg = 1). In order to obtain a sizable contribution to the Higgs boson mass, we 
choose A| = I Msusy^- As one can see from Figure 4.2, the viable regions for the singlino dark 
matter are 55.5 GeV < mg < 62.7 GeV and 0.005 < Agff < 0.034. In Figure 4.6, these regions 
correspond to the red-shaded band and the region between the two blue lines respectively. One 
can see that the singlino resonant dark matter scenario is successful with tan/3 ~ 0(1) and 
Msusy ~ 0(10) TeV. On the other hand, the green band represents 125 GeV < mh < 126 GeV. 
We also find that the current dark matter relic abundance and the Higgs boson mass can be 
explained simultaneously with tan/3 ~ 0(1) and Msusy ~ 0(10) TeV. The blue (dark blue)- 

’^^The one-loop Aefr can be roughly estimated as Agg°°^ ~ m^usy which is negligible in comparison 

with Aeft®. 
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Figure 4.6: Contours of rris (red lines), Aeff (blue lines) and rrih (black dashed lines) in MgusY- 
tan/3 plane assuming A = Amax at each point. On the red-shaded region ( 55.5 GeV < nis < 
62.7 GeV ), the resonant annihilation via the Higgs boson can occur. The green-shaded region 
satisfies 125 GeV < nih < 126 GeV. The blue (dark blue)-shaded region is excluded by the 
current limits from LUX [156] (XENON [155]). The yellow-shaded region is excluded by the 
Higgs invisible decay search at the CMS [151] and the magenta dashed line is the current bound 
by the global fit of the Higgs coulings [152]. 


shaded region is excluded by the current limits from LUX [156] (XENON [155]). The yellow- 
shaded region is excluded by the Higgs invisible decay search at the CMS [151] and the magenta 
dashed line is the current bound by the global fit of the Higgs couplings [152]. Note that, in 
the calculation of these experimental bounds, we have assumed rrih = 125.5 GeV at each point 
in the plane. It means that these experimental abounds exclude the low-scale and tan (3 0(1) 

nMSSM. 

If we choose the lower value of the green line moves to left because the Higgs boson 
mass obtains more contribution from the quartic coupling (see Eq. (3.63)). On the other hand, 
with smaller value of mi 2 + XAx{S) the singlino mass becomes lighter and the red-shaded region 
moves to right (see Eq. (4.44)). The blue lines are not sensitive to the choice of mf 2 and Ax, 
because Agg is determined by the soft SUSY breaking scale and tan /3. The important point is 
that in any case with MgusY ~ 0(10) TeV and low tan/3 the current dark matter abundance 
and the measured Higgs boson mass can be realized simultaneously. This opens a window for 
the singlino dark matter in high-scale supersymmetry. 
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Figure 4.7: Contours of ms (red lines), Aefr (blue lines) in MsusY-tan/3 plane under mh = 
125.5 GeV by changing A, 0 < A < Amax- On the purple line, the singlino relic abundance Qsh"^ 
is consistent with the current value, = 0.1199 [128]. In the light blue region, 

The left side of the blue (green) dashed line can be probed by the future dark matter search 
LUX [157] (XENONIT [158]). As the future sensitivities of XENONIT, we have used 0.14 (left), 
0.326 (middle) and 0.66 (right) to the parameter /^r. ILC [154] can cover the left side of the 
magenta dashed line. Other lines are the same in Eigure 4.6. 


Finally, in Figure 4.7 we show these regions in detail. The input parameters are the same as 
Figure 4.6 except A. All soft SUSY breaking parameters except Ax are set to Msusy (Acj? = 
cs = 1), and we take = |Msusy^- The A is varied at each point to fix the Higgs mass 
boson to be 125.5 GeV. The range of varying A is 0 < A < Amax- Therefore, in this figure 
the Higgs boson mass is fixed to be 125.5 GeV except the dark-shaded regions, where one can 
not explain mh = 125.5 GeV. The singlino relic abundance is consistent with the current 
value on the purple line, Vlch? = 0.1199 [128]. In the light blue region, . The 

left side of the dashed lines can be covered by LUX (blue) [157], XENONIT (green) [158] and 
ILC (magenta) [154]. As the future sensitivities of XENONIT, we have used 0.14 (left), 0.326 
(middle) and 0.66 (right) to the parameter 

Again, we show that the current dark matter relic abundance (purple line and light blue 
region) and the observed Higgs boson mass (not dark-shaded regions) can be explained simul¬ 
taneously with tan/3 ~ 0(1) and MgusY ~ 0(10) TeV. In addition, we find that the future 
experiments, especially the direct dark matter search by the XENONIT, can probe a sign of the 
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singlino dark matter. 

There are other proposals of the direct dark matter search whose future sensitivities are 
comparable to or stronger than prospect of XENONIT. One of them is XMASS experiment, 
which is located at Kamioka Observatory. The XMASS experiment also use an ultra-pure liquid 
xenon and its scattering with dark matter. An expected sensitivity of the XMASS1.5 (XMASSII) 
project is that spin independent nucleon-dark matter cross section is 10“^® (10“^^) [cm^] at 

moM ~ 60GeV [167,168] (cf. [cm^j (XENONIOO), [cm^] (XENONIT), IQ-^^ [cm^] 

(LUX85days) and 10“^® [cm^] (LUXSOOdays)). Other ones are DarkSide-G2 experiment and 
LUX-ZEPLIN (LZ) experiment. Their expected sensitivities of the spin independent nucleon- 
dark matter cross section are 10“^^ [cm^j (DarkSide-G2 [169]) and 10“^® [cm^] (LZ [170]) at 
m-DM ~ 60 GeV. 


4.5 Discussions 

The NMSSM is another model of the singlet extension of MSSM [100]. As we have discussed 
previous chapter, this model is imposed the discrete Z 3 symmetry and the superpotential is 
given as 

W^NMSSM = \SH 2 H 1 + ^53 + IPYukawa. (4.50) 

In the NMSSM, the singlino can obtain a radiative correction to the mass in addition to the 
tree-level mass ~ 2k{S). The singlino resonant dark matter scenario may be successful 

with small tan /3 and small k in high-scale SUSY scenario. In the small k limit, a singlet-like CP- 
odd scalar boson Ai becomes a pseudo Nambu-Goldstone boson because of the existence of the 
global U(l) Peccei-Quinn symmetry. Therefore, one may be able to make a distinction between 
the singlino resonant scenario in the nMSSM and NMSSM by the search for h —>• AiAi [143]. 


Since there are some new sources of CP violating phases in the nMSSM, the EDM are gen¬ 
erally generated through relative phase between and Mgaugino at the one-loop level. We 
have estimate the one-loop electron EDM, which is generated from chargino-sneutrino loop and 
neutralino-selectron loop. The expected electron EDM is (see Eq. (2.106)) 
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' \ 2 
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(4.51) 


where cj) = arg (/iefrAIgaugino)- One can obtain jdej ~ 0(10“^®) e cm with tan /3 ~ 0(1), MgusY ~ 
0(10) TeV and sinep ~ 0(1). Interestingly, the electron EDM of this size does not conflict with 
the current bound [81] and can be probed by some future experiments [94-96]. 
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4.6 Conclusion of the Resonant Singlino Dark Matter 

In this chapter, we have studied the phenomenology of the singlino resonant dark matter scenario. 
We find that including one-loop corrections to the neutralino masses, the singlino can explain 
the current dark matter relic abundance through the resonant annihilation via the Higgs boson, 
if the soft SUSY breaking scale is high scale. We have shown that with high-scale SUSY breaking 
~ 10 TeV and low tan /3, the dark matter relic abundance and the SM Higgs boson mass can be 
explained simultaneously in this scenario. 

Even for the high-scale SUSY, we have also shown that the parameter region where the 
singlino dark matter is consistent with the current dark matter relic abundance can be probed 
by the future experiments (see Figure 4.2, 4.7). Therefore, the singlino dark matter signal can 
be “a first sign” of the high-scale supersymmetry. 
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Chapter 



Towards a Scale Free Electro weak Baryogenesis 

This chapter is based on the work by the author [18]. We propose a new 
electroweak baryogenesis scenario in high-scale SUSY models. We consider a 
singlet extension of the minimal SUSY standard model introducing additional 
vector-like multiplets. We show that the strongly first-order phase transition 
can occur at a high temperature comparable to the soft SUSY breaking scale. 

In addition, the proper amount of the baryon asymmetry of the universe can be 
generated via the lepton number violating process in the vector-like multiplet 
sector. The typical scale of our scenario, the soft SUSY breaking scale, can 
be any value. Thus our new electroweak baryogenesis scenario can be realized 
at arbitrary scales and we call this scenario as a scale free electroweak baryo¬ 
genesis. This soft SUSY breaking scale is determined by other requirements. 

If the soft SUSY breaking scale is 0(10) TeV, our scenario is compatible with 
the observed mass of the Higgs boson and the constraints by the electric dipole 
moments measurements and the flavor experiments. Furthermore, the singlino 
can be a good candidate of the dark matter. 

5.1 Electroweak Baryogenesis in the nMSSM 

In this chapter, we focus on the baryon asymmetry in the universe and the electroweak baryoge¬ 
nesis mechanism in the nMSSM. This chapter is based on the work by the author [18]. We will 
propose a new electroweak baryogenesis scenario in the nMSSM with high-scale SUSY breaking. 
First, we briefly summarize the our electroweak baryogenesis scenario in the nMSSM. 

The Electroweak baryogenesis (EWBG) [171-173] is one of the most promising mechanisms 
to generate the baryon asymmetry of the universe (BAU) r] = ns/s = (0.86 ± 0.01) x 10“^*^ ~ 
10“^^ [128]. In this mechanism, the first-order phase transition of the Higgs field occurs and 
the bubbles are nucleated initially. Then the CP asymmetric distributions of the particles are 
generated around the bubble walls if there is a source of CP asymmetry. Finally, these CP 
asymmetric distributions turn into the BAU due to the decoupling of the sphaleron process. 
This phase transition which associates with this sphaleron decoupling effect is called as the 
strongly first-order phase transition. 

Within the standard model, this EWBC mechanism can not be realized by two reasons. 
First, the strongly first-order phase transition can not occur while maintaining the Higgs boson 
mass 125 CeV [174,175]. Second, there is no CP-violating source enough to generate the proper 
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amount of the baryon asymmetry [176-178]. Thus, this mechanism requires new physics which 
can cause the strongly first-order phase transition with new CP-violating sources. The typical 
scale of this new physics seems to be comparable to the eleetroweak scale since this mechanism is 
supposed to occur around the eleetroweak scale. Now, the new physics models with such a rela¬ 
tively low scale suffer from severe constraints from the collider searches, the EDM measurements 
and the flavor experiments. 

In this chapter, we propose a new EWBG scenario in which EWBG occurs at arbitrary 
scales. As a new physics model, we consider SUSY models which have a new physical scale, 
the soft SUSY breaking scale, MgusY- In this new scenario, the particles with the masses of 
O(Msusy) play important roles. When the temperature of the universe drops across O(Msusy)) 
the appearance of the universe changes drastically. Eirst, the dominant terms of the potential 
for the scalar fields change from the thermal terms to the soft SUSY breaking terms. Second, the 
particles with masses O(Msusy) disappear due to the Boltzmann suppression. These changes 
deform the shape of the potential for the Higgs fields and they may cause the strongly first-order 
phase transition at the temperature O(MgusY)- In this mechanism, the value of MgusY is not 
constrained. Thus, EWBG can be realized at arbitrary scale MgusY if there is a proper amount 
of the GP-violating sources. 

We consider the nMSSM [13-15] specially. The potential of the nMSSM is suitable for 
the first-order phase transition. The ordinary EWBG scenarios in the nMSSM have been well 
studied in the literature [129,179]. In our new scenario, we add extra vector-like multiplets 
to the nMSSM which are coupled to the singlet superfield. In addition, we introduce a lepton 
number violating term in the vector-like multiplet sector. 

Here, let us see the outline of our scenario. In this scenario, the singlet scalar held obtains 
sizable thermal potential from the vector-like multiplets only at high temperatures. Then, the 
absolute held value of the singlet scalar held becomes smaller at high temperatures than at the 
zero temperature. As a result, the potential for the Higgs held gets deformed. Eurthermore, 
the global minimum of the potential for the Higgs held is generated far from the origin when 
the temperature is around MgusY- At this time, the strongly hrst-order phase transition occurs 
from the origin (symmetric vacuum) to this minimum (breaking vacuum). Subsequently, the 
baryon(H)-|-lepton(L) number is generated . After the strongly hrst-order phase transition, 
the Higgs held is trapped at the breaking vacuum. As the temperature decreases below MgusYj 
the breaking vacuum is lifted up and disappears. Then, the Higgs held returns to the symmetric 
vacuum. In this interval, non zero B — L number is generated from the B + L number by 
the lepton number violating term. As a result, the BAU is not washed out by the sphaleron 
process at the symmetric vacuum. The lepton number violating process is active only when 
T > Msusy since the number densities of the vector-like multiplets get Boltzmann-suppressed 

concrete estimation of the B + L number generated by the hrst-order phase transition is beyond the 
scope of this chapter and it is devoted to future work. 
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when T < MgusY- Thus, the BAU is generated and fixed at the temperatures smaller than 
MgusY- Finally, the Higgs field goes to the electroweak symmetry breaking vacuum when the 
temperature becomes the electroweak scale. 

In this scenario, the whole processes occur at T ~ MgusY- Surprisingly, the scale MgusY 
becomes a free parameter up to the small electroweak scale corrections which are needed to 
realize the electroweak symmetry breaking vacuum. Thus we call this scenario as a scale free 
electroweak baryogenesis. On the other hand, the favored value of the scale MgusY can be 
determined by other experiments. Considering the Higgs mass 125 GeV [134,135] and SUSY 
flavor/CP problem, MgusY ~ 0(10) TeV seems to be favored. Moreover, the singlino, the 
fermionic component of the singlet superfield, can be a good candidate of the dark matter. 
With Msusy ~ 0(10) TeV, the proper amount of the singlino dark matter can be obtained 
by resonant annihilation via the exchange of the standard model Higgs boson [17]. We show 
that the lifetime of the singlino dark matter is long enough even though there is the lepton 
number violating term which induces its decay. Therefore, this scenario can realize the proper 
Higgs boson mass, the right amount of the dark matter and the BAU without SUSY flavor/CP 
problem if MgusY ~ 0(10) TeV. 

5.2 The nMSSM with Vector-like Matters 

In this section, we briefly introduce our model, the nMSSM [13-15] with vector-like multiplets. 
We show the matter contents, the symmetries and the interactions in our model. 

The superpotential and the soft SUSY breaking terms of the nMSSM are given in Sec¬ 
tion 3.4.1. In addition, we add extra vector-like multiplets to the nMSSM. These vector-like 
multiplets play important roles. First, they give the sizable thermal corrections for S to cause 
the first-order phase transition. Second, they give the lepton number violation at high tempera¬ 
tures. As the vector-like multiplets, we add 16 (Q^ U', D' , L' , E' , N') -|- 16 {Q' , U', D' , L', E' , N') 
multiplets (with SO(IO) notation). We express the MSSM multiplets as Qi,Ui, Di, Li, Ei with 
i = 1, 2, 3 denoting the generation. In order to forbid unwanted terms of the vector-like multi¬ 
plets, we impose additional Z3 and Z2 discrete symmetries (see Table 5.1). Z3 symmetry forbids 
the terms like S^LH 2 which cause a rapid decay of the singlino, the dark matter candidate in 
our model (see Sec. 5.6 for details). Z2 symmetry is the vector-like multiplet parity where all 
vector-like multiplets are assigned as odd while the other multiples are assigned as even. We 
consider the situation where this vector-like multiplet parity Z2 is slightly broken and the small 
mixings between the vector-like multiplets and the MSSM multiplets exist. 

The allowed superpotential by the symmetries Z^, Z3 and Z2 in the vector-multiplet sector 
is 

Wsym = XiS (O'Q' + fj'u' + b'D' + L'L' + E'E' + V'V') 

+ kiL'HiE' + k2L'HiN' + ksQ'Hib' + k^Q'H 2 U ', (5.1) 
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Table 5.1: The charge assignment. 
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where we take a universal coupling Ai for SX'X' type terms for simplicity. There are corre¬ 
sponding soft SUSY breaking terms like A-terms Ax^SX'X', A^^L' HiE and soft mass terms 
|Y'p . As mentioned above, we assume that the vector-like multiplet parity Z 2 
is slightly broken The terms which appear after the broken of Z 2 are 

= e%S {Q'Qi + fjiU’ + bit)' + t’Li + biE’^ 

+ C (QiHib' + Q'Hibi + QiH2U' + Q'H2Ui + UHiE' + L'HiE^ 

+ €nN'^ . ( 5 . 2 ) 

We set partially universal couplings e^, e* and cat for simplicity. In this chapter, we consider 
the superpotential 

W = lUyukawa + WnMSSM + Wsyra + > (5-3) 

where WVukawa is the ordinary Yukawa terms in the MSSM superpotential. There are also the 
soft SUSY breaking terms for the MSSM multiplets like the soft masses for the stops m?. 

The lepton number (L) and the baryon number {B) of the vector-like multiplets are set as 
follows. Q', U', D', L', and E' have the same quantum numbers as the corresponding MSSM 
multiplets. X' has the opposite charge of X'. The lepton number of the N' is decided by the 
term k 2 L'HiN' to conserve the lepton number: N' has the same quantum number with E. Note 
that the term only ejyN'^ violates the lepton number explicitly. 

In this model, a singlino which is the fermionic component of the singlet superfield can be 
a good candidate of the dark matter [17]. However, the singlino has a finite lifetime in this 
model since the i?-parity is slightly broken due to the ejsfN'^ term. In Sec. 5 . 6 , we show that 
our eleetroweak baryogenesis scenario is compatible with the singlino dark matter scenario. 

5.3 Overview of our scenario 

In this section, we present the overview of our scenario. Since there are several steps in this 

scenario, we briefly outline the series of the thermal history below. The details of each step are 

7?-symmetry Z|* is also broken softly. Though, we assume that the terms introduced by the broken of 
Z 5 * are negligible except the tadpole terms of S. In addition, we assume that the size of these tadpole terms are 
still 0{Msvsy) with our setup. 
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5.3 Overview of our scenario 



Figure 5.1: The outline of the thermal history of our scenario. The details are given in the text. 


given in the subsequent sections. 

Figure 5.1 shows the rough sketch of the thermal history in our scenario. Each graph shows 
the potential for the Higgs field and the graphs are aligned from left (i) to right (v) as time 
goes. The shaded circle indicates the field value of the Higgs field. T denotes the temperature 
of the universe and B (L) denotes the baryon (lepton) number in the universe. H is the Hubble 
parameter at each time point. The Hubble parameter during the radiation dominated era is 
given as Eq. (4.22). F^p/j is the effective sphaleron rate where the sphaleron process changes the 
B + L number with conserving the B — L number only if F^p/i > H. The situation F^p/j > H is 
realized when the held value of the Higgs held is smaller than the temperature (see Eq. (5.56)). 
F^ is the effective lepton number decreasing rate coming from term. The lepton number 

violating process which changes the L number is active only > H. This condition > H 
corresponds to T > Msusy- If T < Msusy, the number densities of the vector-like multiplets 
are suppressed exponentially since their masses are O(Msusy)- As a result, this lepton number 
violating process would be decoupled since this process is caused by the scattering (or decay) 
processes of the vector-like multiplets (see Eq. (5.60)). 

Here, we briehy outline the thermal history (see Figure 5.1). 

(i) At enough high temperatures compared to O(Msusy)) the potential for the Higgs held is 
lifted and the Higgs held exists at the origin of the potential (symmetric vacuum). Both 
F^p/j and Tj/, are larger than H. At this time, B = L = 0 holds since there is no conserved 
number in the thermal equilibrium. 

(ii) As the temperature decreases, the global minimum(breaking vacuum) of the potential for 
the Higgs held appears far away from the origin. The hrst-order phase transition of the 
Higgs held occurs at T = Tist- Note that both the temperature Tut and the held value 
of the Higgs held at the breaking vacuum are O(Msusy)- At this time, EWGB occurs 
and the B + L number is generated in the interval of tewbg [171-173]. In the interval 
of Tewbg-, does not work {1/tewbg b^) and the B — L number is not generated. 
On the other hand, the held value of the Higgs held at the breaking vacuum is larger than 
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the temperature in this scenario. It makes the sphaleron rate smaller < H at the 
breaking vacuum. Thus the sphaleron process is decoupled and generated B + L number 
is not changed at the breaking vacuum. 

(hi) After EWBG, the Higgs field is trapped at the breaking vacuum. During this time, the 
sphaleron process is decoupled (r^p/j < H). On the other hand, the lepton number violating 
process is active (F^ > H) and the L number decreases gradually. Thus, the B number is 
conserved and the generated B L number is converted to the B — L number. 

(iv) At T = Troll ^ Msvsy, the breaking vacuum (the local minimum of the potential for 
the Higgs held) disappears. Then the Higgs held returns to the symmetric vacuum again 
through the second-order phase transition. The sphaleron process becomes active again 
(F^p/i > H) since the Higgs held exists at the symmetric vacuum. On the other hand, the 
lepton number violating process becomes decoupled due to the Boltzmann suppression of 
the vector-like multiplets at this time (F^ < H). As a result, the generated B — L number 
is conserved. Thus the B number and L number are hxed in the thermal equilibrium. 

(v) After the temperature becomes lower than the eleetroweak scale 0{vew), the Higgs held 
settles down at the eleetroweak symmetry breaking vacuum. At this time, both the 
sphaleron process and the lepton number violating process are decoupled. Thus, the 
generated B — L number is conserved and the BAFi exists until today. 

In this scenario, there are two nontrivial points. 

• The strongly hrst-order phase transition of the Higgs held occurs at Tist 0 {Ms\jsy)- 

• The lepton number violating process is active only when T > O(Msusy) 

The hrst point is discussed in Sec. 5.4. The second point is discussed in Sec. 5.5. In these 
sections, we show that these conditions are satished actually. The essential point is that the 
typical scales of the system such as the potential and the masses of the relevant particles are all 
O(Msusy)- On the other hand, the scale Msusy is not constrained by this scenario. Thus, we 
call this scenario as a scale free eleetroweak baryogenesis 

In addition, the singlino dark matter scenario [17] can be compatible with this scenario. This 
fact is nontrivial since the i2-parity is explicitly broken due to the lepton number violating term 
in our model. Fortunately, the lifetime of the singlino is long enough and the singlino can be a 
good candidate of the dark matter, as we show in Sec. 5.6. 

do not consider the CP-violation sources explicitly. The estimation including them is devoted to future 

work. 
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5.4 Strongly First-Order Phase Transition 

In this section, we show that the strongly first-order phase transition of the Higgs field occurs 
at T ~ O(Msusy)- In Sec. 5.4.1, we introduce the relevant potentials. Sec. 5.4.2 is devoted to 
the intuitive understanding of its behavior. In Sec. 5.4.3, we analyze the full potential defined 
in Sec. 5.4.1. 

5.4.1 Full Scalar Potential 

In this chapter, we consider the following potential 

T) = Vo{4>i) + Vcwi<P^) + Vricbi, T) , (5.4) 

where [i = l,2,s) are the field values of Vq, Few and Vt are the tree-level, the 

Coleman-Weinberg and the thermal potential respectively. 

Here we assume some conditions to make the potential simpler since the complete one- 
loop potential is highly complicated. First, only 0(1) couplings are taken into account. Thus, 
we neglect the MSSM Yukawa couplings except the top Yukawa coupling yt. We also do not 
consider e couplings which are introduced by the broken of the vector-like multiplet parity Z 2 
(see Eq. (5.2)). The couplings of the Higgs field with the vector-like multiplets are assumed as 
k = ki = k 2 = 0(1) and A: 3 ,fc 4 <C 1 to make the potential simple (see Eq. (5.1)). Then, the 
superpotential becomes 

2 

lEpot = ytQ3H2U3 + XSH 2 H 1 + ^5 

+ Ai5 (q'Q' + U'U' + b'b' + L'L' + E'E' + iV'iV') 

-F kL'Hit + kt’HiN' . (5.5) 

Second, we partially neglect the H 2 and S dependences of the one-loop potential. As we will 
see later, the strongly first-order phase transition occurs in tan /3 ~ 0 direction and these de¬ 
pendences are irrelevant. Third, we set all A-terms to be zero and some soft SUSY breaking 
masses to be the same values for simplicity. Fourth, we assume that the scalar components of 
the vector-like multiplets are heavy enough and their effects to the thermal self energy can be 
neglected. 

Here we show the each potential Vq, Vqw and Vt- 

Vq: We can write the tree-level potential from the superpotential Eq. (5.5) and the soft terms 
Eq. (3.28) as 

rnj. Q(f>l 2ts4>s + (f>^(pj. -I- X"^, (5-6) 

’^'^The CP-violating sources can enter in T-terms. However, we do not consider them since we show the possibility 
of the strongly first-order phase transition at high temperatures in this chapter. The study with explicit CP- 
violating sources can be found elsewhere. 


- 75 - 






Chapter 5. Towards a Scale Free Eleetroweak Baryogenesis 


where 


= —nil cos^ (3 — m\ sin^ /3 + 771^3 sin 2/3 , (5.7) 

\2 -2 

= — sin^ 2/3 + ^ cos^ 2/3 , (5.8) 

4 8 

and (j? = (l)\ + (j) 2 , tan/3 = </) 2 • 

Vcw: For the Coleman-Weinberg potential, we consider the terms from the top/stops 
from the vector-like multiplets 

= + (5.9) 


Each term has the form as 


Nc 

327r2 


E 

_ 2 =scalars 




(5.10) 


where Nc is the color factor. Mj’s are the masses of the corresponding particles, 
can be written as 


^cw = 


327r2 

Mt = yt4>2 , 


E^7 


± 






-2Mf 




(5.11) 

(5.12) 

(5.13) 


where Mt is the mass of the top quark and Mj ^ are the diagonalized masses of the stops 
with given 4‘,(j)s- Here, we assume the universal soft mass m? for the left- and the right- 
handed stops. Note that in this potential we have neglected the stop At term and D 
terms, and we have taken account of the filed dependence on the singlet scalar. These are 
different points from Eq. (2.51) in Section 2.2.2. 

Eor the vector-like multiplets, we can diagonalize the mass matrix analytically with the 
assumption written in Sec. 5.4.1. Thus can be written as 


T/vec _ 

'^CW — 


1 

327r2 


2 E 

±4=1,2 





where Msi±, Ms 2 ± and are the diagonalized masses of the vector-like particles. The 
mass matrices of the vector-like matters are given in Appendix D.1.1. The diagonalized 
massed of the vector-like particles are, 

Msi± — 2 (fi \J (ini' “ + k'^4>iY + A:2^2^2^ , 

(5.15) 
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^s2± = ^ i^^L' + '^N' + ± ^ (ml, - mj^, - k?4>lY + 4Af A:2(/.2,/,2^ , 

(5.16) 

i + k^4>l ± ^Jk^cj>i + 4Xjk^l<Pl^ . (5.17) 

Here, we assume = m|,, rnj^, = rnj^, = m\i = m|,,. 

Vt: For the thermal potential, we consider the improved one-loop thermal potential. It means 
that the thermal self energy for all scalars and the longitudinal components of the gauge 
bosons are taken into account. Thus we consider the following set of the thermal potential 


Ft(-/>*, T) = F^((/), T) + F/((^, 0., T) + F#(0.,T) + F)7“(0, T). (5.18) 


Each term have the form as X^j=particies '''^here Cj’s are the numerical 

constants and is defined as [180] 

^(x, T) = dz In (l ^ , 

+ 17> forboson(B), 

+ (ic < 1) for fermion(F). 

is the improved one-loop thermal potential for the Higgs field coming from the Z-boson, 
the W-boson and the top-quark. 


(5.19) 

(5.20) 


F/( 


h,T) = W,i{Mt/T,T) + 


‘^-rB 


iMwlT,T) + -F" (Mz/r,T) 


1 

+ 3 




[Mw/T,t)+^V,^[Mz/T, 


T 


(5.21) 

(5.22) 


where j2, M| = ^^0^/2, -|- 19g‘^T'^/6 and M| = M| -|- 19g'^r^/6-|- 

59g''^T^/18. defined as [180]. Note that if cp <T holds, the Higgs field (p obtains 

thermal mass terms: 

Vr^ip), T) ~ sin^ /3T2 + ^ {2f + g^^) cp"^ . (5.23) 


Vr^ comes from the thermal loops of the charged Higgs boson and the CP-odd Higgs 
boson. We have to take this effect into account since a relatively light charged/CP-odd 
Higgs boson is favored to induce the first-order phase transition. This term can be written 
as 

V^icP, (Ps,T) = Vi (Mcharged/T, t) + (iFodd/T, t) , (5.24) 

-^n h a r gfid = -|- 2A^(/)^ -|- ^cp'^ + Hyl , (5.25) 
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(5.26) 

(5.27) 


-^^odd ~ ''^1 “I" ''^2 + 2A + A 0 + Ha , 

7? r? 9/2 X2 l2 

n = —T2 + -I- ^rp2 _|_ ^ rp2 

4 2 4 3 6 


where Mcharged is the mass of the charged Higgs boson and Modd is the mass of the CP-odd 
Higgs boson. 

is the one-loop thermal potential for (pg coming from the colored vector-like fermions 
and the Higgsinos as 


vS{^s,T) = 24V(^ (Ai,^,/T,r) + AVl (A(/),/r,T) . 


(5.28) 


The second term comes from the Higgsinos and we neglect their small mixing to the singlino 
and the gauginos. Note that if < T holds, (pg obtains the thermal mass terms: 

vSiPg^T) ~ (^A?r2 + Pl . (5.29) 


is the one-loop thermal potential coming from the vector-like multiplets L', L', E', 
E',N',N'. This potential can be written as 

p-“ = 2 Vi{Mg,±/T,T)+4Y,Vi{Mj±/T,T). (5.30) 

i±,i=l,2 it 

Mgi± can be obtained by the replacement of m|, —t m|, -|- 3g‘^T‘^/8 + k'^T'^/6 and mf^, —t 
m\i + PT‘^/3 in (see Eq. (5.15, 5.16)). Here, we neglect the corrections of order 

0{g'‘^T‘^) in the thermal self energy. 


5.4.2 Tree-Level Analysis including Thermal Mass Terms 

In this section, we give the intuitive understanding of the potential. We consider the simplified 
potential which has only the tree-level terms and the thermal mass terms. As the thermal mass 
terms, we include the terms Analysis including the full terms is written in the next 

subsection. Here, we show that the potential is deformed due to the thermal mass terms for 
the singlet field pg. We also show that the global minimum of the potential for the Higgs field 
appears far away from the origin only at high temperatures. 

The potential with only the tree-level terms and the thermal mass terms 14r-i-th can be written 
as 


Vtr+t^iP, P, Pg,T) = - M^)P^ + ivlT^ + ml^)Pl + 2tsPg + , (5.31) 

where ^ sin^ /3 -|- |(2g^ + 5^) and Vs = Xi + ^ (see Eq. (5.23,5.29)). The field value of the 
singlet scalar field can be driven from the minimization condition dVtr+th{Pi,T)/dpg = 0. It is 
derived as 


_ts_ 

0 + X'^p'^ + y|T2 


~ 0{Msusy) ■ 


(5.32) 
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Since ts ~ ^(-^susy) the absolute field value of the singlet scalar field becomes 

0{Msusy)- Note that it decreases when the field value of the Higgs held </> or the tempera¬ 
ture T increases. This is one of the key features of our model. After substituting the held value 
of the singlet scalar held, the potential becomes 


/3, T) = 



For convenience, we rewrite the potential as the following form 



s 


= a(/3, TfX^ + (-6(/3, Tf + c(/3, Tf) X - 


(5.34) 


where 



(5.35) 



XHg 



Note that f{T) > 1 holds. In addition, a{l3,T), b{l3,T) and c{j3,T) are increasing functions with 
respect to T. 

From here, we consider the following conditions: 

(i) Only the electroweak symmetry breaking vacuum is realized at the zero temperature. 

(ii) The global minimum of the potential for the Higgs held appears far away from the origin 
at high temperatures. 

For simplicity, we mainly consider two directions. One is the direction with /lyac being the angle 
of the vacuum at the zero temperature. The other is the direction with /3tr being the typical 
angle of the hrst-order phase transition. As we will see later, /3tr ~ 0 is favored to realize the 
hrst-order phase transition. 

Zero temperature conditions 

First, let us consider the conditions to have only the electroweak symmetry breaking vacuum at 
the zero temperature. 

For the /3vac direction, in order to realize the electroweak symmetry breaking vacuum prop¬ 
erly, we need 


■9Ftr+th(<(>,/3vac,H = 0) I 

d(j) 


(5.37) 
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ayto+th(^,/3vac,T-0) . ^ 

d(p l<t>=rEW ' ^ 

where z;ew — 174.1 GeV is the vacuum expectation value of the Higgs field at the zero temper¬ 
ature. ~ 0 indicates that (j) is at the vicinity of the origin. These conditions can be rewritten 
as 


5(/3vac,0) > 1, (5.39) 

&(/3vac, 0)^ = XewY 

~ l + 2XEw^(a(/3,ac, 0)2-1), (5.40) 

where Xew = I'^S- Not® ^EW ^ 1 holds since we assume the soft SUSY breaking 

scale is much larger than the eleetroweak scale. In order to satisfy these conditions, we need 


a(/3vac j 0) > 1 , 

^(/3vac; 0) = 1+0 


9 

m% 


In addition, we impose the condition not to generate the minimum at /3 

= 0 ). 


d(j) 


I 


>0 for /3 9^ /3v 


This condition can be rewritten as 


(5.41) 

(5.42) 


(5.43) 


6(/3,0) < 1 for /3 /3vac • (5.44) 

For the /3tr direction, there should be no global minimum at the zero temperature. Thus, 
the condition Utr-i-th(^</’) Ar, 0 ) — Utr+th(0) Ar, 0) > 0 is imposed and can be rewritten as 

(a(Ar,0)-l)VAAr,0)'<l- (5.45) 


High temperatures conditions 

Next, let us consider the conditions to have the global minimum far away from the origin at 
high temperatures. 

Suppose that at the critical temperature Tc, two minima of the potential appear at the 
origin and at Y = Xc,l3 = Ar- The condition becomes 

v{Xc, Ar, Tc) = u(0, Ar, Tc), (5.46) 

A(Yc,Ar,Tc) = 0, (5.47) 

where the prime means the partial derivative by Y. To have the positive solutions of Xq and 
Tc, the necessary and sufficient conditions are 

o(Ar,Tc) < 1, 

(a(Ar,Tc)-l)V6(Ar,rc)2 >1. 
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Solutions 

Let us see that the conditions Eqs. (5.41), (5.42), (5.44), (5.45), (5.48), (5.49) can be satisfied 
simultaneously. We divide these conditions to the pairs of Eqs. (5.42, 5.44), Eqs. (5.41, 5.48) 
and Eqs. (5.45, 5.49). 

First, we see the conditions Eqs. (5.42, 5.44). To satisfy these conditions simultaneously, let 
us parameterize b{/3, 0) as the following form 

6(/3, 0) = 6i + 62 cos(2/3 — 263) . (5.50) 

Note that 61 , 62 and 63 are the function of mf, m^, mi 2 , ts and m^Q (see Eqs. (5.7, 5.36)). 
If we take these values to satisfy 61 + 62 — Ij ^2 > 0 and 63 ~ Pvnc, these conditions can be 
satisfied easily. 

Second, we consider the conditions Eqs. (5.41, 5.48). Note that the conditions Eqs. (5.41, 5.48) 
are the opposite conditions. In addition, a(/3, T) is an increasing function of T. Thus, the two 
conditions Eqs. (5.41, 5.48) can not be satisfied with only one direction. However, these con¬ 
ditions can be satisfied with the two directions /3vac and Ptr are needed. Next, let us see that 
Ptr ~ 0 is favored. Note that if the ratio a(/3vac) 0)/a(/3tr) 0) is larger, it is easier to satisfy the 
two conditions Eqs. (5.41, 5.48) at the same time. On the other hand, if > g‘^/2 holds, a(/3, 0 ) 
can be parameterized as 


a(/3, 0) = ai — 02 cos(4/3), (5.51) 

with 01,02 > 0. Thus, if /3vac is near 7 r/ 4 , /3tr ~ 0 is favored to give the ratio o(/3vac) 0)/o(/3tr, 0) 
larger and satisfy these two conditions. 

Finally, let us consider the conditions Eqs. (5.45, 5.49). The discrepancy between the con¬ 
ditions Eqs. (5.45, 5.49) can be reconciled by f{Tc). In other words, the thermal mass of 
4>s can work to generate the global minimum of the potential for the Higgs field only at high 
temperatures. Actually, if we find the values of o(/3tr,0), 6(/3tr,0) and f{Tc) which satisfy 

(o(/3tr,0)-l)2 + 6(/3tr,0)2 <1, (5.52) 

[f{Tcf/^a{f3tr,0) - 1 )' + f{Tc)^b{f3tr,0f > 1, (5.53) 

the conditions Eq. (5.45, 5.49) can be satisfied. 

The above solutions can be achieved simultaneously with appropriate parameters. Thus the 
global minimum far away from the origin can be generated only at high temperatures due to 
the thermal mass for the singlet field (ps- Note that small value of a(/3tr,0) and large value of 
6(/3tr,0) are favored in order to satisfy the above conditions. Small a(/3tr,0) is satisfied easily 
with tan /3vac ~ 1- On the other hand, large value of 6(/3tr, 0) corresponds to small 171^2 compared 
to |mf -|- m||. This situation makes the charged Higgs boson light. As we will see in the full 
potential analysis of the next subsection, the strongly first-order phase transition can actually 
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Table 5.2: The parameters at the benchmark point. 


tan /3vac 

A^ 


k 

is/mio 

m'ilmio 


Q/msfi 

2.0 

0.50 

0.50 

1.0 

0.58 

-0.1657 

-0.1675 0.001226 

1 


occur at the high temperature. In addition, the region with tan/3vac ~ CI(1) and the light 
charged Higgs boson is favored. 

5.4.3 Numerical Analysis with Full Potential 

In this section, we analyze the full potential introduced in Sec. 5.4.1. We show that the strongly 
first-order phase transition can actually occur at the temperature comparable to Msusy- At 
first, we show the thermal history at a benchmark point. Next, the conditions for the strongly 
first-order phase transition are discussed. Then, we present a scatter plot and show that the 
region with low tan/3vac and the light charged Higgs boson is favored in our scenario. 

Thermal history 

Now, let us see the typical thermal history of our scenario. Table 5.2 shows the benchmark 
parameters. We take tan/3vac = 2.0, = 0.50 and k = 1.0. The standard model coupling 

constants have scale dependence. We take the values at the scale 10 TeV: = 0.753, = 0.528 

and = 0.394. For simplicity, we assume that all of the soft SUSY breaking masses are same 
m? = q. In order to realize the eleetroweak symmetry breaking vacuum at 

the zero temperature, q) corrections are needed. However, such small corrections are 

negligible for the high temperature dynamics. Thus we do not consider the corrections Note 
that we have checked that at the benchmark point there is no Landau pole of couplings k up to 
the GUT scale using one-loop RGEs [181]. 

Figure 5.2 shows the potential for the Higgs field as a function of 4> with varying 

temperatures T. cpg and tan/3 are calculated to minimize the potential for each given ej) and T. 

Typically, 4>s/msfi j'NJ - 0.5 and tan/3 ~ 0.01 — 0.1 hold. At the high temperature (the red line 

T/nisfi = 0.4), the origin is the only minimum of the potential. As the temperature decreases, 
a global minimum appears at (p/rngfi ~ 0.4 (see the orange line T/rrisfi = 0.37). Then, after 
T/nisfi = 0.31 (the cyan line), the potential is lifted up and the local minimum disappears at 
T/msfi = 0.17 (the black line). 

Note that nisfi can be any value in this analysis. If nisfl is varied, the size of the corrections 
q) changes. In addition, the values of the standard model couplings change since 

impose the zero temperature conditions as a(/3va,c,0) > 1 and &(/lvac) = 1 (see Eq. (5.41, 5.42)). In order 
to impose these conditions easily, we absorb the tadpole and quadratic terms of the Coleman-Weinberg potential 
into the tree parameter. The explicit formulae of the tadpole and quadratic terms are given in Appendix D.l 
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Figure 5.2: The potential for the Higgs field T) as a function of 4> with varying temper¬ 

atures T. (ps and tan/3 are calculated to minimize the potential for each given cp and T. Here, 
we subtract the constant term from the potential to set V^\^{(p = 0,T) = 0. The region with 
small (p is enlarged in the right figure. 
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their values depend on the scale to calculate. However, up to these small corrections, the results 
do not depend on the value of rrisfi. Thus we call this scenario as a scale free eleetroweak 
baryogenesis. 

Strongly first-order phase transition 

Here, we show the conditions for the strongly first-order phase transition. 

First, let us see the condition for the first-order phase transition to occur. If the global 
minimum of the potential exists except the origin, the vacuum tunneling from the origin to the 
minimum can occur. We call this global minimum as the breaking vacuum and the origin as the 
symmetric vacuum. The finite temperature vacuum tunneling rate Ttran. per unit space-time 
volume V is given as the following form: 

El^ ^ ^ (5.54) 

where S{T) = S 3 /T and S 3 is the three-dimensional Euclidean action which is evaluated on the 
bounce solution [182,183]. The condition for the first-order phase transition to occur is given 
by 

J = 1. (5.55) 

For T ~ O(TeV), the first-order phase transition occurs at S{T) < 130 [184]. Here, we adopt 
the condition S{T) = 130 for the first-order phase transition to occur. Since this condition has 
only a logarithmic dependence on the temperature, we ignore this dependence. We show this 
condition of the first-order phase transition explicitly in Appendix C.1.2. 

Second, we show the condition for the strongly first-order phase transition. After the vacuum 
tunneling occurs, the Higgs field is trapped at the breaking vacuum. To cause EWBG, the 
sphaleron process have to be decoupled at the breaking vacuum since the B + L number should 
not be washed out. The sphaleron rate is evaluated as [185] 

r\42/27T A4> 

TsphOcTe 3 T ^ (5.56) 

with A(/) = + 4>2 at the breaking vacuum. In order to decouple the sphaleron process, 

Tsp/i <C is required. This condition is equivalent to Ac/i/T > 0.9 [129], which is derived at 
T ~ 0(100) GeV. Since this condition has only a logarithmic dependence on the temperature, 
we adopt the condition A(f>/T > 0.9 as the strongly first-order phase transition 
Eigure 5.3 shows S{T) and Aeff/T as a function of T. The three-dimensional 
bounce solution S{T) is analyzed numerically by CosmoTransitions software package [187]. 
The thick lines correspond to the benchmark point. The condition for the strongly first-order 

’^®With higher temperature, the condition value 0.9 becomes smaller [186]. Thus, the condition A(^/T > 0.9 is 
conservative. 
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Figure 5.3: The classical action S{T) for the three-dimensional bonnce solution 

{Blue) and ^(p/T (Red) as a function of T. We take Af = 0.49 {dashed), 0.50 {thick) and 0.55 
{dotted). The gray line represents S{T) = 130 and A(/)/T = 0.9. 


phase transition is lA.(j)lT > 0.9 when S{T) decreases to 130 at the first time. Note that the 
temperature T decreases as the time goes. From the Figure 5.3, we can see that the action 
S{T) becomes smaller than 130 at the first time with ~ 1.1 when T/rngfi ~ 0.34 at the 

benchmark point. Therefore, the strongly first-order phase transition occurs at this time. Then 
the B + L number is generated by the EWBG process and the BAU is generated thanks to 
the lepton number violating process (see Sec. 5.5). The other lines are drawn with the same 
parameters at the benchmark point except Ai. Note that the action value S{T) is sensitive 
to the parameter Ai. With larger Ai, the thermal effects on the (ps become stronger. Then, 
the potential gets more deformed. As a result, the action value S{T) and A(p/T\s=i 3 o become 
smaller. With A^ = 0.55, Acp/T is not larger than 0.9 when S{T) becomes 130 at the first time. 
Thus, the phase transition is not strong. On the other hand, with smaller Ai, S{T) does not 
decrease to 130. The strongly first-order phase transition occurs with 0.50 < Af < 0.55 for 
the benchmark point. Figure 5.4 indicates the profile of the bounce solution at the benchmark 
point. From this figure, we find that the typical wall width is L^jT ~ 30, and A/3 ~ 0.1. 


Scatter plot 

In order to show the favored region in our scenario, we present a scatter plot in the plane of 
tan/3vac and the charged Higgs boson mass Mcharged (Figure 5.5). We have scanned the following 
parameter ranges. 


200 GeV < Mcharged < 2 TeV , 
1.5 < tan/3vac <10, 
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TjTYis o~ 0.34 
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Figure 5.4: The bounce solution profile for the first-order phase transition at the benchmark 
point. The horizontal axis is the space coordinate r normalized by rngfi- r = 0 corresponds to 
the center of the bubble. The blue lines indicate the field values of 4>i, 02 and —0^. The red 
line corresponds to tan/3. 

0.3 < Af <1.0, 

0.5 < ts/ml^ <0.7, (5.57) 

with fixed values k = 1.0, = 0.5, m? = = mj^, = m^Q. We have estimated the bounce 

action by a simplified way in which we use the one-dimensional scalar potential 14nin(0; T). The 
fitting formula of the Euclidean action is given in Appendix C.2. We have checked the error 
of the fitting formula of the Euclidean action is at most ~ 20% compared with the results of 
by the three-dimensional scalar potential. Here we do not consider the mass of the standard 
model Higgs boson. It depends on Ax (for low tan /3 region) and m? (for large tan /3 region) 
which do not change our result so much. Thus when MgusY = 0(10) TeV, we can obtain the 
standard model Higgs boson mass 125 GeV easily with varying Ax or m? [57]. Therefore, we 
set rrisfi = 10 TeV here and do not consider their effects. At all points in Eigure 5.5 the first- 
order phase transition {^S{T) < 130) occurs. At the green points the strongly first-order phase 
transition {A4>/T > 0.9) occurs. We find that the region with low tan/3 and the light charged 
Higgs boson is favored in our scenario. This is consistent with the intuitive understanding in 
the previous subsection. 

5.5 Baryon Asymmetry of the Universe 

So far, we have seen that the strongly first-order phase transition can occur at a high temperature 
in our scenario. In this section, we show that the proper amount of the BAU can be generated. 
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Figure 5.5: The scatter plot in tan /3vac — Tfcharged plane. The red shaded region shows the 
exclusion region by B Xgj search. At all points the first-order phase transition (^S(T) < 130) 
occurs. At the green points the strongly first-order phase transition (Acp/T > 0.9) occurs. The 
star corresponds to the benchmark point. 

Here, we consider the lepton number violating process caused by ejvN'^. This process is needed 
for the BAU to exist until today since the generated B+L number by EWBG should be converted 
to the B — L number. 

The lepton number (L) and the baryon number (B) of all multiplets are defined in Sec. 5.2. 
It is important that only the term violates the lepton number. To make the discussion 

clear, we define N' number by the approximate U(l) symmetry N' : 1,N' : —1. This N' number 
is contained in the L number via the mixing e (see Eq. (5.2)). Note that the masses of the 
fermion and the scalar components of N', N' are O(MgusY)- 

Here, we see the L number decreasing process and the thermal history of our model to show 
the generation of the BAU. 

Lepton number decreasing process 

Let us see the details of the L number decreasing process initially. There are two steps in this 
process. At first, the L number in the standard model sector is converted to the N' number via 
the mixing terms e. Then, the N' number decreases due to the term ejyN'^. For simplicity, we 
consider the situation that the rate of the former process is larger than that of the latter one 
with assuming e > cat. Thus, the bottleneck process of the L number violation is the process 
caused by Therefore, we consider this process only and denote the rate of this process 

as T pi;/. 
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At first, let us see F^/ to estimate the effective L number decreasing rate We assume 
that the scalar component of N' does not decay to two fermion components of N' kinematically 
for simplicity Then, the N' number is violated only by the scattering processes. For 
T < 0(Msusy)j Fcan be estimated as 


F^'(r) ~ 


Cat (-^susyF")^^^ 

levr M|usy 


and the equilibrium N' number density njv' obeys 


MsusY \ 

T ) ’ 


(5.58) 


dnp^i 

dt 


—3Hn]\f' — F+ (lepton number conserving processes) . 


(5.59) 


Then the effective L number decreasing rate F^ can be written as 


TpU) = 

'^L 


MsusY 
IGvr Ni 



2Msusy \ 

T J • 


(5.60) 


The is the equilibrium L number density and W ~ 10 is the number of the light components 
of the leptons^^. Here we use njy/ = (l/6)(MsusYT)^^^exp(- MsusY/T)(fiL/T) and np = 
(N,/6)T3 (ml/T) [188]. 


Thermal history 

Now, let us consider the thermal history of our scenario (for the overview see Sec. 5.3). At 
first, the strongly first-order phase transition of the Higgs field occurs at Tist. The Higgs field 
is trapped at the breaking vacuum of the potential until T = Troii. Then the Higgs held returns 
to the origin again. At the benchmark point, Tist — 0.34 rus^o and Tmii — 0.15 rngfi hold. If 
> H holds during T > Troii and < H holds during Troii > T, the BAU exists as explained 
below. 

At the time T = Tist, the strongly hrst-order phase transition occurs. In our model, we 
assume that the B + L number is generated at this time. This EWBG process is supposed 
to occur in the time span tewbg and typically tewbg ^ holds. Since we consider the 
situation T^{Tist) is the same scale of H(Tist), the effects of F^; can be negligible. Then the 
B + L number is generated with the B — L number unchanged 

yB{Tist) + YL{Tist) > 0 , (5.61) 

YB{Tist)-YL{Tut) = ^, (5.62) 

the benchmark point, the scalar component of N' can decay to two fermion components of N'. This is 
avoided by choosing the coupling XiSN'N' larger and the SUSY breaking mass term for the scalar component of 
N' smaller. This choice does not change the result of the previous section. 

’^^Strictly speaking, a linear combination of the B number and the L number decreases by the N' number 
violating process. However, the amount of the BAU is changed only by a factor of a few with this effect. Thus, 
we do not consider this effect for simplicity. 
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where Yb/l is defined as the baryon/lepton number density divided by the entropy density. 

After EWBG, the Higgs field is trapped at the breaking vacuum during Tist > T > Trou- 
The sphaleron process decouples because the field value of the Higgs field is larger than the 
temperature. As a result, the B number conserves. On the other hand, the L number gradually 
decreases due to the L number violating proccess. The L number decreasing factor A"dec can be 
estimated as 


AUtoU) 

hfdec = / f'^dt 



HT 


dT. 


Thus, just before T = Tfoii, the L number and the B number become 


(5.63) 


Tl(T roll) ^ e-^^^^YL{Tut), (5.64) 

Tb(Tj.o1i) = YBiTist ). (5.65) 


After the Higgs field returns to the origin at Troii > T, the sphaleron process becomes active 
again. Note that the sphaleron process makes the B + L number wash-out towards the thermal 
equilibrium with conserving the B — L number. On the other hand, the B — L number decreases 
by the L number violation process e^- The decreasing factor can be estimated as 





HT 


dT. 


(5.66) 


Then the B number and the L number follow 


YB(Tf) + YUTi) oc YB(Tf) - YL{Tf) , (5.67) 

Yb^) - YUTi) ~ (YBiTrou) - ^^(Troii)) , (5.68) 

where Tf is the temperature at the sufficiently late time MgusY S> Tf and c = {ul — tib)/ riL is 
an 0 (1) factor. 

At the end, T = T^, the B number and the L number are estimated as 

YB(Tf) ~ d-^ • (1 - YB(Tut ), (5.69) 

YUTi) ~ -c-i • (1 - e-^dec) YB(Tut ), (5.70) 

with d = (ub — nL)/nB. If all particles except the standard model particles are heavy enough, 
c = 79/51 and d = 79/28 hold [188]. In order to obtain the sizable BAU, A"dec 1 and <C 1 

are favored (see Eqs. (5.69, 5.70)). This corresponds to T^ > 77 during T > Troii and < H 
during Trou > T. Note that both A^ec and N^j are proportional to and the quantity Adec/Au, 
is a function of Ti^f, Troii and MgusY- Thus, if Adec/A^, 3> 1 holds, the sizable BAU can 
exist until today since we can find the suitable value of ej^ which makes large N^ec (Afdec 1) 
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and small ^ 1). At the benchmark point, we obtain N^ec/^w r\j 30 To ensure 

-^dec ^ 1) we can choose ew ~ 10“^. With this choice of ew, the generated baryon asymmetry 
at the EWGB exists until today. In general, Adec/-^to 1 can hold since there is a hierarchy 
Msvsy > Tist > Troll- Thus, the BAU can exist by this mechanism in our scenario. 


5.6 Singlino Dark Matter 


In this section, we show that the singlino dark matter scenario is compatible with our new 
baryogenesis scenario. At first, we briefly review the properties of the singlino dark matter 
as we have shown thoroughly in the previous chapter. Then, we estimate the lifetime of the 
singlino dark matter with the lepton number violating term. We show that it does not suffer 
from experimental constraints. 

Let us review the singlino dark matter scenario (the detail is written in the previous chapter). 
In our model, after integrating out the particles with masses above the eleetroweak scale, the 
low energy effective Lagrangian becomes 

CeS = ^SM - ^ss - ^hss , (5.71) 

where h is the standard model Higgs boson and Tsm is the standard model Lagrangian. Here, 
s is the singlino, the lightest neutralino mainly composed by the fermionic component of the 
singlet superfield S. We denote the singlino as the Majorana spinor. The effective coupling Aeff 
can be estimated as 


Aeff sin 2/3. (5.72) 

JWSUSY 

The singlino mass nig is dominated by the one-loop corrections when Msusy is large. In our 
model, the singlino can get sizable corrections from vector-like multiplets sector. The singlino 
mass can be evaluated as 


mg ~ tt^ALsusy • (5.73) 

(dvr)^ 

In this model, the singlino can be a good candidate of the dark matter. If mg ~ 60 GeV 
and Aeff ~ 0(0.01), the singlino dark matter scenario is successful with resonant annihilation 
via the exchange of the standard model Higgs boson. Such a situation can be realized when 
AfsuSY ~ 0(10) TeV, tan/3 ~ 0(1) and A, Ai ~ 0(1). Note that the low tan/3 and 0(1) 

’^^This value depends on the value of the exponential factor in Eq. (5.60). Here, we set this exponential factor 
as the typical masses of the vector-like fermions Msusy = ts/rn^fi — |</)s| . However this exponential factor 
also depends on the masses of the vector-like scalar bosons since at least one boson particle participates in the 
scattering process. Typically these masses are heavier than |(/)s| and this exponential factor becomes larger. Thus, 
we have chosen the conservative value here since the ratio Ndec/N^ becomes larger with larger exponential factor. 
’^^’^Strictly speaking, the singlino mass is promotional to the T-terms {Ax-^SX'X') which are dropped off in the 
previous discussions. However, the effects of such A-terms ~ Msusy to the thermal dynamics are supposed to 
be small and do not change the previous results. 
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couplings are realized with our baryogenesis scenario. The soft SUSY breaking scale Msusy is 
determined by the requirement of the singlino dark matter scenario, especially by the effective 
coupling Aeff Eq. (5.72). 

In our model, there are the lepton number violating term (cat) and the SM-extraparticles 
mixing terms (e). Thus, this model does not conserve the i2-parity and the singlino can decay 
to the standard model particles. So, let us estimate the decay rate of the singlino. Note that the 
term breaks the lepton number by three AL = 3. In addition, the decay process breaks 

the vector-like multiplet parity Z 2 at least three times. 

Let us consider the dominant decay channel s —)• zzz/i/. The other channels are more sup¬ 
pressed since the number of final state particles increases if the decay products include charged 
leptons. To see the coupling of the sfui/, we consider the following fermion four point operator 
which arises from integrating out the particles whose masses are O(Msusy) 






i\D 

-''■'SUSY 


(5.74) 


Here, fg^uu is a numerical factor and e denotes e* or defined in Eq. (5.2). We denote as 
the Weyl spinors. The decay rate of the singlino due to this operator can be evaluated as 


r(5 




p2 ,6 


^ J S, 


2 


m~ 


30727r3 M|usy 


(5.75) 


The mass of the singlino is favored to be mg ~ 60 GeV in order to realize resonant annihilation 
via the exchange of the standard model Higgs boson. On the other hand, the typical value of 
the Z 2 breaking couplings e is 0(10“^) (see Sec. 5.5). Thus the lifetime of the singlino Tg can 
be estimated as 


Tg ~ 0.8 X 10^® 



10"^ 

e 


6 


fsvuv 


2 


/ Msusy 
VlO TeV 


/ 60 GeV 

V mg 


5 

[sec] . 


(5.76) 


Now, we estimate the upper bound on the factor fgyuv by a diagrammatic way. Let us 
consider the diagrams for the operator Ogyyy. To draw the diagram, we need the vertex . 

Thus, each diagram includes the vertex ew and three propagators of N'. Since the final state 
contains three neutrinos, these three propagators of N' should be converted to them. Therefore, 
there are three lines which start from N' to the neutrino. We call these lines as lepton lines. 
For each lepton line, at least one propagator of a Higgs multiplet or one vacuum expectation 
value of the Higgs field vew should be attached If vew'^ are attached to all three lepton 
lines, the diagram may have no loops and fgyuu is suppressed by {vew /■ If vew’’^ are 
attached to two lepton lines of three, the diagram has at least one loop and fgyuv is suppressed 
by (l/167r^)(uEvy/-^SUSY)^- If vew is attached to one lepton line of three, the diagram has at 

There are also the diagrams in which some lepton lines have no propagator of a Higgs multiplet and no 
vacuum expectation value of the Higgs field. However, such a diagram is highly suppressed since a lot of vertices 
are needed. Therefore, we ignore such diagrams here. 
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least one loop and fsuvi/ is suppressed by {1 /\Q'k‘^){vew/^S\JSy)- If ^evk’s are not attached to 
any lepton lines, the diagram has at least two loops and fswu is suppressed by (l/lGvr^)^. In 
any cases, the following inequality holds 

huuu < 10-^ , (5.77) 

if Msusy = 0(10) TeV. Note that this estimate of the upper bound on fsuvu is conservative. 

From Eq. (5.76) and Eq. (5.77), the lifetime of the singlino becomes long rg > 10^® sec. . On 
the other hand, there are experimental bounds on the lifetime of the dark matter. First, the 
lifetime of the dark matter should be much longer than the lifetime of the universe 10^^ sec. . 
Second, there are constraints from the cosmic ray searches, tdm ^ 10^® sec. [189]. Obviously, 
the lifetime of the singlino is much longer than the experimental bounds Thus there is no 
problem in the decay of the singlino and the singlino can be a good candidate of the dark matter 
in our scenario. 

5.7 Discussions 

We comment on the experimental constraints for the light charged Higgs boson and the SM- 
extraparticles mixings. First, let us consider the constraints for the light charged Higgs boson. 
The relatively light charged Higgs boson and heavy SUSY particles are favored in our scenario. 
It means that this model can be regarded as the two-Higgs doublet model at low energy regions. 
Even if SUSY particles are heavy, the existence of the light extra scalars is constrained by the 
flavor and the CP violation physics. In the viable parameter region of our scenario, the process 
H —>■ Xs 7 is the only relevant constraint from the flavor physics. The red shaded region in 
Figure 5.5 is excluded at 95 % C.L. by a current bound [190]. On the other hand, the electron 
EDM is one of the severe constraints on a new CP-violating phase [81]. In our scenario, a new 
CP-violating phase may enter into the potential for the Higgs field through only Ax term. The 
electron EDM is induced by the mixing between the CP-even and the CP-odd Higgs bosons 
which is estimated as ~ \^{ts/m^sQ){^Ax/msfi)- If Ax/rrisfl < 0.1, our scenario is compatible 
with the current bound of the electron EDM experiments [191]. 

Second, the flavor changing neutral current appears through the SM-extraparticles mixings. 
One of the severe constraints comes from the branching ratio of /i —)• ey, Br(/i —>• ey) < 5.7 x 
10“^^ (90 % CL) [79]. We have estimated this value at our benchmark point, 

Br(/r ^ ey) ~ x 10“®. (5.78) 

Thus if we take e < 10“^, the bound from Br(/r —)• ey) can be escaped easily. 

The experimental bounds by the cosmic ray searches come from the various decay channels of the dark 
matter. Especially, the decays to the charged leptons are important. However, in our model, the decays of the 
singlino to the charged leptons are more suppressed than the decay to three neutrinos since the number of the 
final states increases. Therefore, the bounds can be evaded more easily. 
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Finally, we comment on the neutrino masses. Although this model includes the matters 
which couple to the neutrinos, neutrino masses are protected to be zero. In order to generate 
the nonzero neutrino masses, we have to extend our model or change the imposed symmetry. 
For example, let us introduce three right-handed neutrino superhelds N" which have R- 
symmetry charge 3, Z 3 symmetry charge 1 and Z 2 parity even. We also introduce an extra 
local U(l) symmetry, and N" has charge 2 and other particles do not have charge under this 
symmetry. If there is an additional Planck suppressed term W = H 2 LN” where (j) is the 

singlet under and Z 3 , an even parity under Z 2 , and has charge —1 under the extra symmetry, 
the appropriate Dirac neutrino masses are generated when the extra symmetry breaking scale 
is ((/>) ~ 10^3 QgY_ 

As this chapter is a hrst study of a scale free electroweak baryogenesis scenario, much work 
is left to be done. First, we have to check whether the proper amount of the baryon number 
can be generated within our scenario including the explicit CP-violating phases. Second, the 
vacuum stability against the charged Higgs held direction has to be checked in detail. 

We have to comment on the stability for the charged Higgs held direction of the potential. At 
the zero temperature, there is a charge breaking global minimum in the charged Higgs direction if 
we consider the tree-level potential at the benchmark point (see Appendix D.2). This is because 
the charged Higgs boson mass can become negative in the relatively large (p region since the held 
value \(ps\ becomes small (see Eq. (3.44)). To see that there is no problem with this minimum, 
we have checked two conditions. First, we have checked that the charged Higgs boson mass 
including the thermal self energy is positive for all time of the universe"^^^. Second, we have 
calculated a tunneling rate from the electroweak breaking vacuum to the charge breaking global 
minimum at the zero temperature with tree-level potential. Then it turned out that the lifetime 
of the electroweak breaking vacuum is much longer than the one of the universe (5 ~ 0(1000)). 
Thus, we consider that this minimum gives no problem. The full analysis of the stability against 
the charged Higgs field direction is complicated and will be done in the future. 

5.8 Conclusion of Scale Free Electroweak Baryogenesis 

In this chapter, we proposed a new electroweak baryogenesis scenario with the high-scale nMSSM 
including vector-like multiplets. We have shown that the strongly first-order phase transition 
can occur in a high temperature comparable to Mgusv- The proper amount of the BAU can 
be generated via the lepton number violating process. Furthermore, the singlino dark matter 
scenario [17] is also compatible with our scenario. The key points are as follows: (i) the thermal 
mass term for the singlet scalar field generates the global minimum of the potential for the Higgs 
held far from the origin, (ii) the lepton number violating process converts the B + L number to 

simplicity, we do not include the mass corrections from the Coleman-Weinberg potential. We have numer¬ 
ically checked that the Coleman-Weinberg potential for charged Higgs boson gives typically positive contributions 
to the mass of the charged Higgs boson.). 
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the B — L number. Even though there is the lepton number violating process, the lifetime of the 
singlino is long enough. In this baryogenesis process, Mgusv can be an arbitrary value and it is 
almost a free parameter. Thus, we call this scenario as a scale free eleetroweak baryogenesis. The 
scale MgusY will be determined by other requirements. If Mgusv ~ 0(10) TeV, this scenario is 
compatible with the proper Higgs boson mass and the right amount of the singlino dark matter 
without SUSY flavor/CP problem [17]. In addition, this singlino dark matter scenario can be 
testable by future experiments of the search of the Higgs invisible decay and the direct direction 
of the dark matter (see previous section). 

Consequently, we have shown the possibility of the high scale baryogenesis scenario. We 
hope that this study becomes a first step of scale free eleetroweak baryogenesis scenarios. 
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6 

Conclusion 


The standard model of the particle physics has worked very well for a long time. Nevertheless, 
there are many unsolved problems within the SM, for example, the observed dark matter particles 
and baryon asymmetry of the universe. From theoretical viewpoint, the gauge hierarchy problem 
is still in question. The supersymmetric models are good candidates of the physics beyond the 
standard model since they can solve the hierarchy problem naturally. 

The minimal SUSY model contains a supersymmetric dimensionful parameter /z, and this 
parameter causes problem, which is also one of the hierarchy problem. In order to realize 
nature, namely the Z boson mass is to be at the electroweak scale, the /r parameter has to know 
the soft SUSY breaking scale. 

The nearly Minimal Supersymmetric Standard Model (nMSSM) is one of the promising 
models of the new physics: this model can avoid the /r problem, the domain wall problem, and 
the tadpole problem simultaneously. In addition, this model has natural candidate of the dark 
matter, namely singlino, and can generate the baryon asymmetry of the universe. 

In this thesis, we consider the phenomenology of the nMSSM. Especially, we focus on the 
phenomenology of the dark matter and the baryon asymmetry in the universe by the electroweak 
baryogenesis mechanism. 

First, we have studied the phenomenology of the singlino resonant dark matter scenario. We 
find that with high-scale supersymmetry breaking the singlino can obtain a sizable radiative 
correction to the singlino mass, which opens a window for the singlet dark matter scenario with 
resonant annihilation via the exchange of the Higgs boson. We have also shown that with high- 
scale SUSY breaking ~ 10 TeV and low tan/3, the dark matter relic abundance and the SM 
Higgs boson mass can be explained simultaneously in this scenario. 

Next, we have also proposed a new electroweak baryogenesis scenario with the high-scale 
nMSSM including vector-like multiplets. We have shown that the strongly first-order phase 
transition can occur in a high temperature comparable to the soft SUSY breaking scale. The 
proper amount of the baryon asymmetry in the universe can be generated via the lepton number 
violating process. Furthermore, we have calculated explicitly the lifetime of the singlino and we 


Chapter 6. Conclusion 


find that the singlino dark matter scenario is also compatible with our scenario. The key points 
are as follows: (i) the thermal mass term for the singlet scalar field generates the global minimum 
of the potential for the Higgs field far from the origin, (ii) the lepton number violating process 
converts the B + L number to the B — L number. Even though there is the lepton number 
violating process, the lifetime of the singlino is long enough. In this baryogenesis process, the 
soft SUSY breaking scale can be an arbitrary value and it is almost a free parameter. Thus, we 
call this scenario as a scale free electroweak baryogenesis. The soft SUSY breaking scale will 
be determined by other requirements. If it is 0(10) TeV, this scenario is compatible with the 
proper Higgs boson mass without SUSY flavor/CP problem. 

Therefore, we hnd that when the soft SUSY breaking scale is 0(10) TeV, this electroweak 
baryogenesis scenario is compatible with the singlino resonant scenario. In addition, these 
scenarios are also compatible with the observed mass of the Higgs boson and the constraints by 
the electric dipole moments measurements and the flavor experiments. 

Even for the high-scale SUSY, we have also shown that the parameter region where the 
singlino dark matter is consistent with the current dark matter relic abundance can be probed 
by the future experiments. Hence, the singlino dark matter signal can be a hrst sign of the 
high-scale supersymmetry. 

As a result of these two studies, we conclude that the nMSSM with a high-scale SUSY 
breaking is valid and can be probed by the direct direction of the singlino dark matter. 
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A 

Notations and Conventions 


A.l Notations 


We use the following metric tensor 




/I 0 0 0 \ 

0-10 0 
0 0-10 
Vo 0 0 - 1 / 


with the Greek indices are // = 0, Ij 2,3,4, and the totally antisymmetric tensor is 


,0123 _ , _ 1 

e — —£0123 — 


The Pauli matrix (t“ is defined as, 


v = 


0 1 


ct 2 = 


0 -i 
i 0 


1x3 = 


1 0 
0 - 1 / ’ 


a 0/ ’ 

and and are defined as follows, 

C7^ = (l, ^), d'^ = (l,-^). 

The antisymmetric tensor for two components is 

£12 = —£21 = 1 - 

The gamma matrices satisfy the following anti commutation relations. 

We use the chiral basis gamma matrices. 


= 


0 

0 / ’ 


= i7°7'7V = 


-1 0 

0 1 


(A.l) 

(A.2) 

(A.3) 

(A.4) 

(A.5) 

(A.6) 

(A.7) 

(A.8) 
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The Gell-Mann matrix A“ is defined as, 



A'* = 


A' = 



.(A.9) 


A.2 Group theoretical constants 


The generators of a simple Lie group SU(A^) are represented by x Hermitian matrix 
which satisfies Tr[t“] = 0, and = cr“/2 (a = 1,2,3) for SU(2), = A“/2 (a = 1,2,... ,8) for 

SU(3). While, in the following the t for U(l)y means the hyper charge Yi of operated field <l>j. 
The generators are normalized by 


Tr 



r(i) 


forU(l)y, 

i for SU(2),SU(3). 


(A.IO) 


The structure constants for a Lie group is defined by 






(A.ll) 


The quadratic Casimir operator for the fundamental representation is defined as, 


=C{i)6), (A.12) 

wehre 


Ci{i) = for U(l) and <l>j, 

(A.13) 

C2{i) = \ 

f I for SU(2) and^i = Q, L, 

[ 0 for SU(2) and$i = U, D, E, S, 

(A.14) 

C3{i) = \ 

f 1 for SU(3) and^i = Q, U, D, 

[ 0 for SU(3) and$i = L, E, S. 

(A.15) 

The quadratic Casimir operator for the adjoint representation is 



E/^e/cde = C{G)d,,, 

a 

(A.16) 

with 

( 0 forU(l), 



C{G) = 2 for SU(2), 

(A.17) 


[ 3 for SU(3). 
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B 

Quantum Corrections 


In this appendix, we collect the functions of the quantum correction, which 
are needed in this thesis. We first present the full set of two-loop RGBs for the 
coupling constants of the SM and the singlet extension SUSY model. Second 
the one-loop corrections to the mass of the neutralino are exhibited. Finally, 
we present the loop functions which we use in the text. 


B.l Renormalization Group Equations 

In this appendix, we assume the high-scale SUSY mass spectrum, vew ^ -^SUSY and Mgaugino ~ 
H ~ Y^rug = ©(Msusy)) where rug represents dimension two soft SUSY breaking mass term of 
Higgs and sfermion. Note that in the split case, Mgaugino ~ /.t MsusYj we should consider 
the RGEs of the Yukawa-like gaugino couplings Eq. (2.97) because the Higgs-Higgsino-gaugino 
coupling is still active at vew < Q < Msusy- Therefore, we should take into account the RGE 
of not only the SM couplings but also the Yukawa-like gaugino couplings at vew < Q < Msusy- 
The RGEs of the split mass spectrum case is written in Refs. [45,57]. 

We write the RGEs in the following notation, 

dgi ^ I3i,i 132,i ^ 

dlnQ (dvr)^ (dvr)^ ’ 

where Q is the renormalization scale. 

B.1.1 RGEs below SUSY breaking scale 

First, we present the RGEs up to two-loop order for the SM couplings, g', g, gs,yt,yb,yT and 
Aquartic in the MS scheme [45,192-194]. In the SM, one-loop level (3 function for the gauge 
couplings is given as 

= -g^^ , (B.2) 

where Uf {ni,) is a number of the gauge multiplet of the Weyl spinor (complex scalar). 


The one-loop /3 functions for the SM couplings are 
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^l,g' 

A,9 
/^l,9s 

f^^,yb 

h,yT 

/3i,A5 


/3 


41 
19 


-g , 


-79I 


f 

17 ,2 

9 9 

yt 1 

1 

1 

1 

1 

1 

4"^ 

f 

5 12 

9 9 

yb 1 

- q — 

12 ^ 

4^' 

/ 

' 15 ,2 

9 , 

5 r ( 

-^y - 

4" 


3 I ^ 2 


.2 , o„.2 


2"® 

^ 2 
2^- 


; , 


'4 I > 


2Aquartic ^SAquartic + ^Vt 4" ^Vb 4“ “^Vt 2 ^ 2^ 

“4(3?/^ + 3y^ + y^) + -g'^ + 


(B.3) 

(B.4) 

(B.5) 

(B.6) 

(B.7) 

(B.8) 


(B.9) 


The two-loop /3 function for SM couplings 

_ ,^3/199 2 9 2 , 44^2 

/32,9' - + 2 ^ + 


are 


P 2 , 


,yb 


z 17 2 5 

- yl/t 


1^2,g 

1^2, gs 

^2,yt = Vt 

+2/b 




s“ + f 9^ + 1292 - I 


5 2 ^2 J- 

-jVt - -jVb- 2 


2 ^2 

■Vb - Y' 

1 


-Xr 


3/11/2,92 oe 2 

+ 2 ^ -265s - 


Xt - 2y! 


Vt 


"131 


/9 255 9 

■5 + ■ " 


7 


255 2 9 2 11 9 9 2 , A 

Yo'd 4- 36y^ — 12y^ — —y^ — OAquarticJ 

/2 I 99 2 . . 2 1 2 I 6 2\ I 2 ^6 ,2 I 16 2 9 2 

. 48 ^ + 16 ^ +X--^yb + -^yr)+yr[Y9 -^5. 

1187 ,4 23 4 4 3 ,2 2 19 ,2 2 22 

Y^g'^ - ys - 1085s - 45''5' + ys 9s + Xgl 


3 2 

4~ 2 Aquartic 4“ 


56 

+yl 


2 /^91 ,2 I 99 2 . 2 1 2 112,62^ 

(,48^ +16^ + ~ 4 ^* " + 4^- j 

2/4 79 ,2 225 2 09 . 2 2 9 2 A 2/^25 

'6 y-^y ^ i^y ^ ^^y^ ~ ~ iy^ ~ + 5 ? y 


(B.IO) 

(B.ll) 

(B.12) 


(B.13) 


/2 2 

■g + -^9 


9 2 
4 "- 


I 8^2 

' 2 Aquartic 


1^2,yr = 5r 5? 


9 / 85 *9 

Vt (+ 


127 
45 


^/4 _ 23 ^4 


5 ^ - 108y: 


9/22, /2 2 I n 2 2 

|5 5 + y 5 5s + 95 5s 


31 


(B.14) 


45 2 00 2 27 2 32 27 

-^5 4- 20yg - —y^ -|- -yb 


yl 


2 /^25 ,2 45 

+56 (^^5 + 

+ 0.2 , 4W ,4 

2 Aquartic 2 ^ ^ 


27 2 27 


9^ + 20 y 2 - ^y^ 

23 4,9 ,22 

Y9 +19 9 


, /179 


^' 2 \, 2 /'‘-'^/ 2 , 165 2 
4 5r ) + 5r ( 5 + ig 5 


- 35 ; 


6Aquartic 
(B.15) 


- 102 - 











B.l Renormalization Group Equations 




Aquartic (-78Aquartic " 727/2 _ ^yl - 2%^ + IS^'^ + 54(/2) 

+ Aquarticy? “ 42y^ + gQ^2^ 

+Aquarticyb + Aquartic^? (^“2/? + 

+A,„„.,C (^9'“ - f 9- + f 9“9") + %? (i59? - 'ivl - 59“ - 169?) 


+yt 


-f9'‘-59‘ + 2l9“9= 


+ 49?(-39? + 159i+j9'"-169? 


+yl (fs" - 59 * + 99'"9") + 49? ( 59 ? - 29 “) + 9? (-f 9« - |/ + ll9'"9") 


379 


-^ 5 '® + -^9^ - 1^9'^9^ - -rrr 


305 


559 


24 


24 


9'V 


289 


24 


9'^9^- 


(B.16) 


B.1.2 RGBs above SUSY breaking scale 

Next, we have calculated the RGBs up to two-loop order for couplings of the singlet extension 
of MSSM (NMSSM), g', g, gs,yt,yb,yT, ^ and n in the DR scheme. The following results are 
consistent with Ref. [100]. Note that in the limit of tc = 0, the following RGBs reproduce the 
one in the nMSSM. R is becanse the tadpole term does not affect the RGBs of the dimensionless 
coupling constants. Burthermore, in the limit of k = 0 and A = 0, the following RGBs reproduce 
the one in the MSSM. 

Before we go forward the derivation of the SUSY RGBs, we comment on the matching condi¬ 
tion on the Yukawa couplings. At SUSY breaking scale, the following relationships are imposed, 

yu,SM{MsuSY) = yu,SUSY{Ms\JSY)smP, 

2/d,5M( A/susy) = yd,SUSY{MsvsY)cos/3, (B.17) 

where yi^susY is the Yukawa coupling in the superpotential, and yi^sM is the one in the SM. 
The higher-order corrections to this matching condition (conversion factor from MS to DR 
regularization scheme and threshold corrections of the heavy sparticles) are given in Ref. [45]. 
For simplicity, we have omitted its subscript {SUSY/SM) in the text. For example, in previous 
section we use the SM Yukawa couplings yi^sM, while in this section we use the SUSY Yukawa 
couplings yi,susY- 


Let us deviate the SUSY RGB of the Yukawa couplings. First we define the superpotential 
as follows. 
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W = (B.18) 

where <hQ are the bare Chiral superfields. Thanks to the non-renormalization theorem, a loga¬ 
rithmic dependence of the coupling constant can be related to a logarithmic dependence of the 
wave function renormalization constant [97,195]. 

vijkii +yiikij +yijnk > (b.19) 

(B.20) 

Lat\ (B.21) 




Vijk 


(4yr)2n 

Pn,Mij 

(4yr)2n 

l^n,Li 

(4yr)2n 


where subscript n represents the loop order, and 7 *- is an anomalous dimension matrix. 


77 




2dlnQ 


(B.22) 

(B.23) 


where are the renormalized Chiral superfields. 

As using the superpotential Eq. (B.18), the anomalous dimension matrix 7 ^ is given as 
follows [196-198], 




1 

(47r)2 

1 

(47r)2 


2l3'rSaCJi)S‘-~,C’ + 


(B.24) 

(B.25) 


where we assume the Chiral superfield 4>* is belong to the fundamental representation of the 
gauge group Ga- Here, we use the well-known one-loop level (3 function for the gauge couplings 
in the SUSY model. 


^ .,3 

(47r)2 (47r)2^“ 



(B.26) 


The two-loop level (3 function for the gauge couplings is also related to these functions [196-198], 


/^2,ga 

(47r)'^ 


(4vr) 




(B.27) 


where r is the number of the generator of the gauge group Ga- 
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B.l Renormalization Group Equations 


In the NMSSM, neglecting all Yukawa coupling except the third generation, the explicit 
formulae of the one-loop anomalous dimension (i = j) Eq. (B.24) are 


(4vr)^7Q^ 

= Vt+yl- :j^9“ - 

8 2 

- 3 "- 

(B.28) 

(4vr)^W 



(B.29) 

(4vr)^7B^ 

9 2/9 89 

= ^yb-gg -^9s, 


(B.30) 

(4vr)^7L^ 



(B.31) 

(4vr)^7£^ 

= 2 y 2 - 2 y' 2 , 


(B.32) 


= 3y^ + y? - ^g'^ - \g^ 

+ a 2 , 

(B.33) 

(4vr)^7S 

= 3yl-\g'^-\g^ + X^ 

•> 

(B.34) 

(4vr)^7s^ 

= 2X^ + 2k^, 


(B.35) 


and off-diagonal one-loop anomalous dimensions {i ^ j) are zero. Using the one-loop anomalous 
dimensions Eqs. (B.28-B.35), one can derive the two-loop level RGEs of the NMSSM in the DR 
scheme. 


The one-loop /3 functions for the NMSSM couplings are 

/3i,5' = 115 '^, 

A ,9 = 

+ + ^y'i+ yi^ ^ 

= yb (^- - 35 ^ - y -h Vt + Qyl + , 

1 ^ 1 ,Vt — yr {—^g'^ — 3 ( 7 ^ -|- -|- 3 yl + 4 y^) , 

Pi,\ = A {—g'‘^ — 3g^ + 4A^ -|- 2k^ + 3yf + 3yl + , 

(3i,k = K (6A^ -h 6 av^) . 


(B.36) 

(B.37) 

(B.38) 

(B.39) 

(B.40) 

(B.41) 

(B.42) 

(B.43) 


The two-loop /3 functions for the NMSSM couplings are 


^2,g' 

1^2,gs 

k,yt 




/199 '2 I Q 2 I 38 2 26 2 14 2 „2 

5 + 9 y —gs - -:ryt - -:ryb “ 62 /r - 2A 




V 9 " ’ ' 3 3 3 

( 3 y '2 + 25 y 2 + 24 y 2 _ _ Qy^ _ 2 y 2 - 2X^) , 


/II 


gt + 95 ^ + I4ys - 4yt - 4y^ 

yt (- 22 yt - 5 yt - 3A^ - 5 yh! - Hyt^^ - yWr - - 2/?A' - 2A2 ac 2 
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+ 2 s“!/? + + 69 ^!/? + 169^9? + + ys'* “ ysJ 

+ 59 'y + ^«+ 89 ^ 91 ), 

fem = m(- 22 yt - 5 yf- 3 yi- 3 X‘- 5 ylyf- 3 yly^- 3 ylyl-iyf\^- 2 \^k^ 
+ l9'M + la'M + 2y’Vr + 69"» + 169=9? + ^ 9 “ + ^9* - f 9? 
+ ^9'V + y9“9? + 89=9?„ 


(B.47) 


(B.48) 


/52,?7r - 


Vr { -lOyt - 9yt - 3A" - Qy^yl - _ 3y2y2 _ ^y 2^2 _ ^^2^2 ^ 25,2^2 


-‘^a'^Vb + ^g^vl + '^^alvl + + y / + , (B.49) 

/32,a = A (^-lOA^ - 9y\ - 9yl - 3y^ - - 9\^yl - 9\^yl - 2>\^yl - 12 \^k^ 

-^aivl + 25 '^A^ + - ^^'^6 + 2 y'^y^ + Gy^A^ + IGy^y^ + IGy^y^ 

+f5'^ + y/ + 3y'V), (B.50) 

P2,K = K (^-24k^ - 12 A'^ - 24^2A^ - l^yJX^ - 18y^A^ - Gy^A^ + Gy'^A^ 

+18y2A2^. (B.51) 

These formulae are consistent with the RGEs of Ref. [100]. 

B.2 One-loop Corrections to the Mass of the Neutralino 

In this section, we collect the one-loop radiative corrections to the mass of the neutralino [1G4]. 

Note that, we have found that Ref. [1G4] includes some typos in the equations of the one-loop 

corrections, and were provided with fixed one-loop corrections by the author [1G5]. 

The self-energy matrix for neutralinos is given as follows, 


_r^ 


a=l 

2 

+2 B„(p. ,y-, Mwjryyg 

a=l 
5 3 

+ E E ^o(P, rn^l , 

a=l b=l 
5 3 


a=l b=2 
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+ ^ Bo(p, Mz)m^org_co,^orfo,GO,xO 

a=l 

2 

+6 ^ Bo(p, mt, 

a=l 

2 

- 8 y^-So(p,m -T^o M/ -^^-ow - 

Xa’ ' Xa Xj,W,Xa Xi ,W,Xa 

a=l 

5 

-4 ^ Bo{p, m^o,Mz)m^or^Sz^^or^o^z,x°^, 


(B.52) 


a=l 


where we have neglected the terms which are proportional to the SM fermion mass except top 
quark. Xi represents the neutralino of not a mass eigenstate but a gauge eigenstate, p is the 
momentum of the external line. 

The one-loop corrections from the redefinition of the neutralino held via the wave function 
renormalization are given as follows, 

e5(p) = 


a=l 
3 3 


a=l 
5 3 


Xa Xi,G+,Xa 




Bi(n, TTii,., mij )r^o-* T^o-*-, m/i, iT^o ob^o o 

XpV*,Vl, Xi,^i.,^b 2^^ Zj,h-b,Xa XiM^a 

a=l 6=1 a=l 6=1 

^53 ^5 
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6 3 2 
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E "^xg’ (b.53) 
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y^Bi(p,m-,mj:^-)r^o ^ _r^o f7+ - - y^Bi(p,m-,Miy)r-^:o „+ _r-^:o^+ - 
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6 3 


■^R 

xlW,Xa 


Xj,K^'^b Xi^K^Ub Z^ Xa ’ x?,^,Xa 

a=l 6=1 a=\ 

- ^ B,(p.»xS>^rfr?«,z,,srg,z,,.. (b.54) 


a=l 


The vertices with the gauge eigenstate neutralino Xi; ^ and Y (rlp^ ) are defined in the 
Ref. [164], 


B.3 Loop Functions 


In the calculation of the radiative corrections to the Higgs boson mass from the Coleman- 
Weinberg potential, we have used the following loop functions, 

1 


f{Q ,x,y) = 


gix,y) = 


x-y 

1 1 


-ifi{Q^x)-h{Q^y)), 


2x-y 


(B.55) 


{x - y) 

1 

4(x — yY 

1 

(x - yf 


{x - yf 


(/2(<3^a:) + /2(<3^y)) 
■[h{Q\x)-h{Q\y)), 


with 


fi{Q^,x) = 2 x(ln^-l 


/ 2 (( 3 ^x) = 41n- 


X 


(B.56) 

(B.57) 

(B.58) 


In the one-loop threshold corrections of the Higgs quartic coupling at high scale, we have used 
the following loop functions, 

2x In X 
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(B.59) 

(B.60) 

(B.61) 

(B.62) 

(B.63) 
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(x 2 - l)(y 2 _ 1 ) 

1 + xy 
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and these functions are normalized such that -F(l) = G(l) = 77(1) = 77i(l) = 772 ( 1 ) 

5(1) = /l/2/3/4(l) = /5/6/7/8(l5 1) = 1- 


(B.64) 

(B.65) 

(B. 66 ) 

(B.67) 

(B. 68 ) 

(B.69) 

(B.70) 

(B.71) 

(B.72) 

(B.73) 

/(I) = 


On the other hand, in the one-loop threshold corrections of the Higgs quartic coupling at weak 
scale, we have used the following loop functions, 


n 2 Q 1 

Fi(g) = 61n^ + -ln^--Z 
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(B.74) 
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where ^ = mj^/M^ cw = cos^vvj sw = sin^vv and 


'W 


(B.75) 


(B.76) 


Z[x] = 


2((arctan(l/(() for x > 1/4 
Cln[(l-7C)/(1-C)] for x< 1/4, 


(B.77) 
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C(x) = V|1 - 4x1- 


(B.78) 


In the calculation of the branching ratio of the pL 67 , we have used the following loop 
functions, 


fc{x,y) = xy 
fN{x,y) = xy 


5 - 3 ( 3 : + y) +xy 
(x - iy(y - 1)2 
-3 + X + y + xy 


21ogx 
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(x - y)(x - 
2x log X 


(x-l)2(y-l)2 [x-y){x- 

and these functions are normalized such that /c(l, 1) = 


2 logy 

- 1)3 (x - y){y - 1)3_ 
2y log y 

1)3 ^x-y){y-ir\^ 

1/2 and /Ar(l, 1) = 1/6. 


(B.79) 

(B.80) 


In the calculation of the one-loop self energy, we have used the following functions^^, which 
are called Passarino-Veltman function [199-201], 
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(B.81) 

(B.82) 

(B.83) 


(B.84) 


where we use the dimensional regulation, n = 4 — 2e, and p is a momentum of the external line. 
After the integration of the loop momentum, Aq function becomes 


Ao(m) = m 2 Q -h 1 - In . 

Here l/e = 1/e — 7 ^; +In dvr and is Euler’s constant (0.57721. ..). 

Bq function becomes 

- fB{x+) - fsix-), 

where 

s{p, mi, m2) ± \/s{p, mi, m2)2 — 4^2(777,2 _ jg) 

= p2 _ j.^2 ^2^ 

= ln(l — x) — xln ^1 — —^ — 1 . 

’^^Our notation of A and B functions are the same as Ref. [199]. 
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B.3 Loop Functions 


The function Bi can be expressed by Aq and Bq as follows, 

1 

The zero momentum of the external line limit, Bq and Bi functions can be expressed as 
follows. 


Bi{p,mi,Tn2) = [Ao{m2) - Ao{mi) + s(p,mi,m2)So(p,mi,m 2 )] • (B.90) 
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Bo(0, mi,m 2 ) = --In 
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1 1 

+ X + - + 


m2 ' 

m? I ’ 
Inx 


2 ' l-x ' (1 -a:)2 
where M = max(mi,m 2 ), m = min(mi,m 2 ) and x = 


— 9{1 — x)\nx 


(B.91) 
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Chapter 


C 

Vacuum Transition 


In this appendix, we briefly review the vacuum decay and the vacuum tran¬ 
sition rate per unit space-time volume at zero/finlte temperature case. Then, 
we give the httlng formula for Euclidean action in four and three dimensions. 

C.l Vacuum Decay 

When the scalar expectation values on the vacuum are at a global minimum of the scalar 
potential, this vacuum is stable. On the other hand, when its scalar expectation values are at a 
local minimum of one, this vacuum becomes unstable. Eventually, the unstable vacuum (false 
vacuum) decay into the global minimum (true vacuum) by a quantum fluctuation (quantum 
tunneling) and by a thermal fluctuation (thermal tunneling) for a scalar field located at the 
local minimum. When the vacuum transition rate per unit space-time volume from the false 
vacuum to the true vacuum is larger than the Hubble parameter, a bubble nucleate and all false 
vacuum decay quickly. Namely the universe is filled with the true vacuum. On the other hand, 
if the vacuum transition rate per unit space-time volume is smaller than the Hubble parameter, 
the lifetime of the false vacuum is longer than the age of the universe and it becomes meta-stable 
vacuum. 

The scalar potential of the SUSY models often has a global minimum which is not an ordi¬ 
nary electroweak symmetry breaking vacuum, and the electroweak symmetry breaking vacuum 
becomes unstable [202-205]. One of the reasons is that the large ^ term and large tan/3 lead to 
the large scalar trilinear couplings like y^/ztan 13 H 2 LE, which can generate the charged breaking 
global minimum. Then, the condition that the electroweak breaking vacuum is not unstable 
gives the upper bound on /r and tan j3. These vacuum meta-stability conditions give an allowed 
region of the deviation from the standard model prediction in SUSY models (e.g. [65,206-209]). 

C.1.1 Quantum Tunneling at Zero Temperature 

A possibility of the quantum tunneling of the false vacuum had been first suggested by Kobzarev, 
Okun, and Voloshin [210] . Then, Callan and Coleman had established a calculation method [211- 
213]. The vacuum transition rate from the false vacuum to the true vacuum can be evaluated 
by semiclassical technique. At this time, the imaginary part of the energy of the false vacuum 
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corresponds to the vacuum transition rate to the true vacuum at zero temperature. In the 
semiclassical technique, one can evaluate the energy of the false vacuum state using the path 
integral method in Euclidean space-time. The vacuum transition rate per unit volume at zero 
temperature is evaluated as, 

= (C.l) 

A precise value of the prefactor A is difficult to evaluate. However, it does not depend dramat¬ 
ically on the parameters of the theory, and one can roughly estimate it at the fourth power of 
the typical electroweak scale in the potential, 

A ~ (100 GeV)^ (C.2) 


In contrast, the power index ^4 is a sensitive parameter of the vacuum transition rate. It can 
be evaluated by an 0(4) symmetric solution as follows, 

S 4 = SeS(p)] - Se41^^], (C.3) 

where p is a radial coordinate in four-dimensional space-time, 

P = \/(i-io)^ + (x-xo)2, (C.4) 


here the bubble nucleate on (to, xq). The Euclidean action in four dimensions SEiifp] is as 
follows. 


Se 4 [ 4 >{p)] = 2vr^ 



+ V{cP^) , 


(C.5) 


where cjii are the real scalar held which construct the scalar potential. Note that if is complex 
scalar held, the factor 1/2 is removed, (j)^ represents the value of the helds at false vacuum. The 
bounce conhguration is a stationary point of the action, namely (p satishes the held equations, 

dV{j)) (fp 3 dp 

dp dp"^ p dp 

In addition, the bounce conhguration also satishes the following boundary condition. 


lim p{p) = , 

p^OO 


dPip) 


dp 


= 0 . 


p=0 


(C.7) 

(C.8) 


On the other hand, the present Hubble parameter is observed as Hq ~ 1.5 x 10-^2 QgY 
When the vacuum transition rate per unit volume Ttrans./G is larger than the fourth power of 
the current Hubble parameter. 


Ae-^* > 


(C.9) 
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the lifetime of false vacuum is shorter than the age of the universe and the false vacuum decay 
quickly. This inequality leads to 


S '4 < 4 In 



~ 400. 


(C.IO) 


Therefore, when the S '4 < 400 the false vacuum must decay into the true vacuum. This is the 
vacuum meta-stability condition at the zero tempreture. 


C.1.2 Thermal Tunneling at Finite Temperature 

Linde had pointed out that the argument of Coleman is valid only at zero temperature [182]. 
It is because that when the temperature is as large as the typically particle scale, the potential 
changes drastically due to the high-temperature effects. 

At the finite temperature, vacuum decay (thermal tunneling) rate can be evaluated by [183, 
184], 

= , (C.ll) 

where the prefactor A^ depends on the temperature. However, similarly to the zero temperature 
case, one can estimate 

At ~ (C.12) 

The power index S 3 is the Euclidean action in three dimensions which is evaluated by an 0(3) 
symmetric solution, 

S3 = SE3mP)]-SE3[^^], (C.13) 


with 


poo 

Se 3 [ 4 >{p)] = dvr / r'^dr 

Jo 


EUf I 


(C.14) 


At the finite temperature, the bounce configuration cf) satisfies the held equations, 

dV ( 0 ) d?^ 2 d^ 

dd) dr"^ r dr ’ 


(C.15) 


and 


lim (/>(r) = (pJ, 

r—)-oo 


d(j){r) 

dr 


= 0 . 

r=0 


(C. 16 ) 

(C.17) 
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The thermal tunneling condition is given as follows, 



dt 




Ttrans. 

V 


> 1 . 


(C.18) 


The Hubble parameter during the radiation-dominated era is (the same as Eq. (4.22)), 





45M2j 




(C.19) 


The equation of motion Eq. (C.19) and the adiabatic expansion condition. 


lead to 


d{aT) 

dt 


dT + aT = 0, 


dt 

Hr 


(T)- = (-It) 

2il4p[^ 

T^' 


(C.20) 


(C.21) 


Substituting Eqs. (C.19), (C.21) into Eq. (C.18), the thermal tunneling condition becomes 


/ 


today 


dT 


16 C^-^pi -g^/T 


j't) 


> 1, 


(C.22) 


namely 





(C.23) 


where Tc is the critical temperature that the false vacuum and the true vacuum degenerate. 
If we take ~ 3 x 10“^ that is the typical value at T > 1 GeV [150], the thermal tunneling 
condition leads to 


53 ^ f 140 for Tc = 1 TeV, 
^'^\l30 forrc = 10TeV. 


(C.24) 


Therefore, we have used S 3 /T < 130 as the thermal tunneling condition in the text. 


C.2 Fitting Formula for Euclidean Action 

According to Ref. [214], we give the fitting formula for Euclidean action in four and three 
dimensions. Note that in Ref. [214], ip denotes real scalar field. On the other hand, we use (j) as 
a complex scalar field. Thus, a factor y/2 or 2 is different from the literature. 
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Let us consider the following potential for complex scalar (j), 

VicP) = \4>\^-cPo,i)i^-M, (C.25) 


where A > 0 and the dimensional coefficients (po,ij(po ,2 > 0. Obviously, this scalar potential has 
two minimum: the origin (p = 0 and another point cp = (pm (~ + '^o, 2 )/ 2 )- Namely, the 

origin is a local minimum (false vacuum), and (p = (pm is a global minimum (true vacuum). If 
the scalar field expectation value is at a the false vacuum, then this vacuum will be metastable 
and will decay into the stable true vacuum. 

For such the one-dimensional scalar potential, the Euclidean action in four dimensions from 
the false vacuum to the true vacuum can be obtained by the following fitting formula [214], 

£'4 = —L_ [13.832<5 - 10.8195^ 2.0765(5^] , (C.26) 

3A (2 — 0 )'* 

with 


^ _ 80o,1(/>O,2 


f2- (C-27) 

( 00,1 + 00 , 2 ) 

Similarly, the Euclidean action in three dimensions from the false vacuum to the true vacuum 
can be obtained by the following fitting formula [214], 


^3 = 


:(0o,i + 00,2) ^y^2 ' [8-2938(5 - 5.53305^ 0.8180(5^] 


81VA 


(2 - 5)^ 


(C.28) 


In our numerical analysis Figure 5.5, we have reduced the three-dimensional to the one¬ 
dimensional scalar potential in the direction of the phase transition. Then we use the following 
polynomial fitting, 


V (0) = (a0^ -I- 60 -t- c)0^, 


(C.29) 


with 


(— 3 ^ 4 ^ -|- 2A)Vm + A^Vm 
^{l-2A^)Vm+A^VM 

{—AA F 3)Vm + A^Vm 
[0m(l-^))]' ’ 


(C.30) 

(C.31) 

(C.32) 


A = ^. (C.33) 

9m 

This scalar potential also has two minimum. The the origin (0 = 0, V{(p) = 0) is a local 
minimum (false vacuum) and {(pm, 10n) is a global minimum (true vacuum). Note that we have 
used this htting function only when Vm < 0 , since we need the fitting formula of the scalar 
potential that the vacuum on the origin can decay into the vacuum on the global minimum. 
This polynomial satisfy the following conditions, 


- 117 - 











Chapter C. Vacuum Transition 


• The origin is a multiple root. 

• The polynomial passes through {4>m, Vm) and {(pm, Kn)- 

• The derivative of V with respect (p at (p = (pm vanishes: V'{<pm) = 0. 

Here we have taken Vm as a local maximum value of the one-dimensional scalar potential. 
Finally, A, (/>o,i and (/>o ,2 can be written by the parameter a, b, c as follows, 


(p0,l, <Po,2 = 


A = a, 

—b ± — 4ac 


2a 


= <pM 


1-2A^ + A^B ±{A- 1)^{2A-IY - A^B 

A^B -3A + 2 ’ 


(C.34) 


(C.35) 


where 


Vm 

Vm' 


(C.36) 
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Detail Calculations for Chapter 5 

In this appendix, we present detail calculations for Chapter 5: the masses 
of the vector-like matters, the Coleman-Weinberg potential for the vector-like 
matters and for the top/stop multiplets, and the conditions of existence of a 
charge breaking minimum in the charged Higgs direction. 

D.l Coleman-Weinberg Potential 

In this section, we collect the Coleman-Weinberg potentials which are used in Section 5. Here, 
we show explicitly tadpole and quadratic terms of the Coleman-Weinberg potential. It is because 
that these terms can be absorbed by the redefinition of the tree parameters (see Section 5.4.3). 


D.1.1 Masses of Vector-like Matters 


First, the vector-like multiplet superpotential is (see Eq. (5.5)), 

Wvec. = \iS{L'L' + E'& + N'N') + kHi{L'E' + L'iV'), (D.l) 


where we assume that Ai and k are real for simplicity. Note that the strongly first-order phase 
transition occurs in tan (3^0 direction in our model. So, we neglect the colored vector-like 
mattes since we assume that they do not have the dependence. In the following, the doublet 
matters L' and L' are denoted by 



The vector-like fermion mass terms are given as, 


-C 


vec.ferm. 


(L';,Z',iv'*,iv') 


-XiS 

0 

V 0 

/ 0 

AiF 



-XiS* 

0 

0 

\ 

0 

kH^ 

0 

k{Hly 

0 

AiS 

* 

0 

AiF 

0 

/ 

XiS* 

0 k{Hyy 

\ 

0 

0 

0 


0 

0 

Ai5* 


0 

XiS 

0 

/ 


/L'A 

L'2* 

N' 

\N'*J 

(L'A 

E' 

\e'*I 


(D.2) 

(D.3) 


(D.4) 


+ {L^,L'^,E'\E') 
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The eight mass eigenvalues are obtained by diagonalizing the vector-like fermion mass matrix. 
The lighter four mass eigenvalues are 

M‘j_ = ^ (^2A?|5|2 + A:2|//0|2 _ ^+ 4AfA;2|S|2|F0|2^ , (D.5) 

and the heavier four mass eigenvalues are 

^ (^2A?|S|2 + A:2|^0|2^y'^4|^0|4 + 4^2^2|5|2|^0|2^ _ (D_6) 


The vector-like scalar mass terms are given as, 


-^vec.scal. — {^i L 2 , N , )A1vec.seal.neutral 


+ (L;*,F,Z;,E'*)AfL.scal.charged 


/ 

L' 2 * 

N' 

\R'*/ 
(L'2\ 

L'f 

& 

\E'*J 


(D.7) 


with 


M 


vec. seal. neutral 


M 


vec. seal. charged 


^ml,+Xl\S\^ 
0 

-\ik{H?yS 
0 


0 

ml, + + Xi\Sf 

0 


Xik{HfyS 

ml, + k^\H?\^ + Xi\S\^ 


XikH^S* 


-XikH^S* 

0 

m%, + k^\H?\^ + Xi\Sf 
0 


(D.8) 


XikH^S* 


m%,+Xi\S\^ 

(D.9) 


0 

ml,+Xl\S\^ 

0 

XikS*Hf 


Xik{H': 

0 


?)*S 


0 

XikS{Hf)* 

m%,+Xl\S\^ 0 

0 m%, + k'^\H?\^ + Xl\Sf 


where we neglect the D term contributions and the i/® dependence because the strongly first- 
order phase transition occurs in tan /3 ~ 0 direction in our model. 

The eight mass eigenvalues are obtained by diagonalizing the vector-like scalar mass matrix. 
When we take the ml, = ml ,, and ml, = ml,, the eight mass eigenvalues are given 

as follows, 


M: 


s, neul=p 


neu2=f 


i (^ml + m%, + 2A?|5|2 + ^ ^{ml - ml, + 

i (^ml + ml, + 2A?|5|2 + klHff ^ ^{ml - ml, - + 4AfA:2|s|2|^0|2^ 
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M. 


s, chal=p 


M. 


s, cha2=f 


i (^ml, + m|, + 2A2|5|2 + ^ ^{ml, - m|, + ^ 21 ^^- 012)2 + 4A2A:2|5|2|/?0|2^ , 

i + m|, + 2A?|5|2 + ^ ^{ml, - + 4,Xjk^\S\^\H0\^^ . 

(D.IO) 


D.1.2 For Vector-like Matters 


When one takes ml, = ml, = mh, = = m?^, = m|.,, the Coleman-Weinberg potential for 

the vector-like matters is given as 

1 


V^^{HlS) = 


327r2 

1 


327r2 




4A/A (in 


0+4Mi(ln(^ 




,(D.ll) 


with 


iV4% = ^ (^2ml, + 2A?|5|2 + Ai2|i/0|2 ^ ^, (D.12) 

where Q is the renormalization scale. 

If this Coleman-Weinberg potential is expanded around the zero temperature vacuum: 

ts 


m = 0 , 


S = so = - 


m 


2 ’ 

s,0 


(D.13) 


this potential becomes as follows, 

T/vec / TjO c^ _ T/vec i T/vec i /T/vec \* , x/vec , T/vec 

'^CWV-^1 ? — '^CW, 0 + '^CW, si + I'^CW, slJ + '^CW, s2 + KcW, h2 


+K 


CW, h4 


+ 


(D.14) 


with 


T/vec 

'^CW, 0 


T/vec 
'^CW, si 


T/vec 
'^CW, s2 


327r2 

Afgp 

27r2 


A 

27r2 


(mi + Af4)= (in ( iA.. 


2„2 

0 




xh 




-ml - Xllln ( ) + {ml + X^) In ( ) (5 - sp), 


+ xl 




2 042 2i / AlSo^ I z' 2 I o\2 2^ 1 + A^Sq 

- 3AiSo In ( -h -i 3AiSoj In 






(D.15) 


(D.16) 

l^-^pp, 

(D.17) 


T/vec 

'^CW, h2 


A 

47r2 


—ml + ml In 


m 


L> 


+ Afs 




+ 2AJsq In 


+AiSc 

Afsg 


l^iT, 
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(D.18) 


T/vec 
'^CW, h4 


327r2 3(m|, + Af Sq)^ 


6Afmf,/So + SA^Sg + In 


4„4 


+ AiSq 




+I2)^m\isl In 


—4 In 


m 


L' 


+ xh 




+ GA^SqI n 


m- 


-j + 


Q' 



10' 

V ) 

3 


l^iT- 


(D.19) 


The generated tadpole term and the generated quadratic mass terms ^2 ^12 

can be absorbed by the redefinition of the tree parameters, ts, fnifl, nrf. 

D.1.3 For Top/stop 

The Coleman-Weinberg potential for the top/stop multiplet is given as 




32ir'-! 


L ± 


± 


In ■ 


M? 






- - - 2Mf In 


Q2 


with 


(D.20) 


Ml^ = ml + y^\Hlf±ytX\S\\H% 

/\ /r2 _ 21 rjO 12 

JXlt - % 1-^21 > 


(D.21) 

(D.22) 


where we take mq = mfj = m?, and neglect the A term and D term contributions. If this 
Coleman-Weinberg potential is expanded around the zero temperature vacuum, this potential 
becomes as follows. 


with 


^cw, 0 
^cw, 2 


^cw(-^l) ) “S') — Few, 0 + ^CW, 2 + 


(D.23) 


167r2"'* 



3 

IGvr^ 



l^iT + 22/?m? 





(D.24) 

(D.25) 


The generated quadratic mass terms 2 can be absorbed by the redefinition of the tree 

parameters, ml and m|. 


D.2 Charge Breaking Minimum 

In the nMSSM, there is often a charge breaking minimum in the charged Higgs direction. Es¬ 
pecially, the benchmark point in the text has a charged breaking global minimum. Thus the 
electroweak symmetry breaking vacuum becomes metastable vacuum. It is because that when 
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the (j) is relatively large, |(/>s| becomes small (see Eq. (5.32)). At this time, small \(ps\ leads to 
the light charged Higgs boson mass, 

^^charged = ml + ml + + ^(/)^ . (D.26) 

Therefore, the nMSSM tends to have the charge breaking minimum in the charged Higgs direc¬ 
tion. Note that, when the coupling A is 0(1), the charged Higgs boson is always lighter than 
the typical mass of the Heavy Higgs boson (see Eq. (3.66)), 

Miarged = . (D.27) 


In the following, we consider the condition that the charge breaking minimum occurs. We can 
take (H^) = 0 at elsewhere using the SU(2) rotation. Thus, the Higgs potential of the nMSSM 
can be expanded as follows (the same equation as Eq. (3.38)), 


Ho 


+ + \H^\^ + \H+\^ - 2|H'?p|i70|2 _ 2\H^\‘^\H+\‘^ + 

8 

+ ^\H^f\H+f + {-XAxSH^H^ + tsS-ml2H^H^ + H.c.). (D.28) 


Now, let us consider the minimization condition for and “77^”. These conditions are 

given as follows. 


QHq 

dH^ 


dHl 


dS 


vev 


<9Ho 

dH+ 



„+|2 



+ |2 
2 I 


\ A ^2 

AA>.s-m 

Vl 



0 , 

V 2 (ml-h X^(vl + s^) + ^(vl -vl + \v^\^) - XAxs— - mfa—) 
V 4 V2 V2 J 

0 , 

s (^m^j + }?{v\ -\- vl -{- I'l’jl^) + {ts ~ 

0 , 

, ( 2 , \2 2 5^/2 2 I -|-|2\ , 5^ 2\ 

(vj) I m 2 + A s - —(ui -V 2 - \v^\ ) + yUi 1 

0 , 


(D.29) 


(D.30) 


(D.31) 


(D.32) 


namely, 

m\ = {ml 2 + XAxs)—-^{vl-vl-\v^\^)-^\v^\^-X^{vl + s^), (D.33) 

Ul 4 / 

fnl = {ml2 + XAxs)—-X^{vl + s^) + ^{vl-vl-\v^\^), (D.34) 

V2 4 
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_ -ts + >^AxViV 2 

+^2 + 1^2 \^) ' 

and when we assume ^ 0, then 

n2 

2 \2 2 I i/ / 2 2 I +|2\ y 2 

m2 = -X s + —{vi -V2- |?;2 I ) - 

where represents the vev of the scalar field. 

Using Eq. (D.34) and Eq. (D.36), one find the following relationship, 


(mi2 + XAxs)— - 
V2 



then, we can get a sin 2/3 as a function of cf, 


sin 2/3 


2(mf2 + A^as) 
(A2 - ^)(/-2 


where we use 


(D.35) 


(D.36) 


(D.37) 


(D.38) 


vi = (/icos/3, (D.39) 

V 2 = (/isin/3. (D.40) 


Next, using Eq. (D.33) and Eq. (D.36), we can get 




9^ 


-1 


ml-ml + ^{vl - vl) - ^vl + X^vl - (m^a + A^as) — 

2 2 Vl 


(D.41) 


It is non-trivial to show a existence of the solution of these conditions Eqs. (D.35), (D.36), (D.42) 
and (D.41). However, when Ax ~ 0, these conditions become a bit simple. Then, Eq. (D.42) 
becomes 


sin 2/3 


2mf2 

(A2-f)02- 


(D.42) 


As one can see, the angle /3 can be expressed as a function of (f for given input parameters, 
A, mf, m|, rrig, ts and mf 2 - In addition, ?;i, V 2 and can also be expressed as a function of 
(j), using the fact that /3 can be expressed as a function of (f). Finally, s can also be expressed as 
a function of 4> (see Eq. (D.35)). 

Thus, we can obtain the vevs (ui, V 2 , s, V 2 ) that may be the charged breaking vacuum 
solution as a function of cj) for given input parameters. At this time, if these vevs can satisfy 
the following three necessary conditions, the solution surely exists, and the charged breaking 
minimum exists for given input parameters. 

Since sin 2/3 < 1, 


2mf2 

A2-4 


<4>^. 


(D.43) 
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Since an absolute value of V 2 is real, 

> 0. (D.44) 

Since we do not use Eq. (D.36) itself, we should impose it, 

ml + \^s^{(t)) - ^ (uf((/>) - vl{(t)) - \v^\‘^{(t))) + = 0. (D.45) 

In fact, the sample point in the text satisfies the above conditions, and there is a charge 
breaking global minimum in the charged Higgs direction at the zero temperature. However, we 
show that the electroweak symmetry breaking vacuum is actually meta-stable vacuum, and its 
lifetime is much longer than the one of the universe (see Section 5.7). 
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